1

我用三种不同的处理方法进行了一个实验，想看看每种处理方法中蝌蚪的存活率是否有显着差异，是否有显着的阻滞效应。我使用“汽车”包进行了方差分析，发现没有块效应，但是治疗之间的存活率存在显着差异。当我尝试执行 Ryan-Einot-Gabriel-Welsch Post Hoc 测试时，我收到以下错误代码：

if (Tprob[ntr - 3] > Tprob[ntr - 2]) Tprob[ntr - 2] <- Tprob[ntr - : 参数长度为零另外：警告消息：在 qtukey(p = (1 - alpha)^(i/ntr), i, df = DFerror) ：产生的 NaN

如果我放弃块效应或制作一个仅包含治疗和蝌蚪存活百分比的新数据文件，我会得到同样的错误。

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
4

0 回答 0