Question

我进行了三种不同疗法的实验，希望了解每种疗法中want的存活率是否存在显着差异，以及是否有显着的阻滞作用。我使用“汽车”包装进行了Anova试验，发现没有阻滞作用，但是治疗之间的生存率存在显着差异。当我尝试执行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）中：产生了NaNs

如果我放弃方块效应或只用治疗和t的存活率制作一个新的数据文件，就会出现相同的错误。

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

如何在REGW.test

0 个答案: