Good morning,
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)
Thank you