我用三种不同的处理进行了实验,想看看每种处理下蝌蚪的存活率是否有显着差异,是否存在显着的块效应。我使用“car”包进行了方差分析,发现没有块效应,但治疗之间的生存率存在显着差异。当我尝试执行 Ryan-Einot-Gabriel-Welsch 事后测试时,我收到以下错误代码:
if (Tprob[ntr - 3] > Tprob[ntr - 2]) Tprob[ntr - 2] 中的错误 <- Tprob[ntr - : argument is of length zero In addition: Warning message: In qtukey(p = (1 - alpha)^(i/ntr), i, df = DFerror) : NaNs produced
如果我删除块效应或创建仅包含处理和蝌蚪存活百分比的新数据文件,我会得到相同的错误。
setwd("C:/Users/garig/Desktop/TubDataAnalysis")
HylaSur<-read.csv(file="Tadpole_Survival_Nocray.csv", row.names='Tank_ID', sep=",") #data file
mod2<-lm(Mean_Hyla_Survival ~ Block + Treatment, data=HylaSur, na.action = na.omit)
Results_mod2<-Anova(mod2, type=3)
Results_mod2
shapiro.test(resid(mod2))
Posthoc<-REGW.test(mod2, "Treatment" ) #Ryan-Einot-Gabriel-Welsch
Error in if (Tprob[ntr - 3] > Tprob[ntr - 2]) Tprob[ntr - 2] <- Tprob[ntr - :
argument is of length zero
In addition: Warning message:
In qtukey(p = (1 - alpha)^(i/ntr), i, df = DFerror) : NaNs produced
我检查了代码,看起来你至少需要有四个级别。这是你的情况吗?