I am trying to run my data to make multiple comparisons. The way I was told originally was that RStudio gives results with Fisher's protected LSD to show means separation. I ran my numbers without changing anything and have the Fisher's LSD results, however some people are wanting me to try to run the same statistics with a Student-Newman-Keuls test. How do I change my code to give me SNK results?
I do not really think this matters much with this question but this is the code I am running:
setwd("/Users/mc1499/Documents/Thesis Measurements/2018 & 2019 Spring") dat<-read.csv("TQR pooled.csv", na.strings = ".") head(dat) block<-as.factor(dat$Rep) trt<-as.factor(dat$Treat) trtdate<-as.factor(dat$DateTreatment) yr<-as.factor(dat$Year) str(dat) RCI0WAT<-lme(X0WAT ~ Treat*DateTreatment,random = ~1|Year/Rep, na.action=na.omit, data=dat) anova(RCI0WAT) predictmeans(RCI0WAT,"Treat:DateTreatment",pairwise=TRUE)