Student-Newman-Keuls Test

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

There's an implementation of SNK in the agricolae package


https://cran.r-project.org/web/packages/agricolae/
I'm guessing there are other packages that have implementations, as well.

2 Likes

This topic was automatically closed 21 days after the last reply. New replies are no longer allowed.