Как исправить код ошибки «аргумент нулевой длины» в REGW.test - PullRequest
Новая
       94

Как исправить код ошибки «аргумент нулевой длины» в REGW.test

1 голос
/ 03 мая 2019

Я провел эксперимент с тремя различными видами лечения и хочу посмотреть, есть ли существенные различия в выживаемости головастиков при каждом лечении, и были ли существенные эффекты блока. Я выполнил Anova с использованием пакета «машина» и обнаружил, что не было никаких эффектов блока, однако существовали значительные различия в выживаемости между процедурами. Когда я пытаюсь выполнить пост-специальный тест Райана-Эйно-Габриэля-Вельша, я получаю следующий код ошибки:

Ошибка в if (Tprob [ntr - 3]> Tprob [ntr - 2]) Tprob [ntr - 2] <- Tprob [ntr -: аргумент имеет нулевую длину Дополнительно: предупреждающее сообщение: В qtukey (p = (1 - alpha) ^ (i / ntr), i, df = DFerror): произведено NaN </p>

Я получаю ту же ошибку, если отбрасываю эффекты блока или создаю новый файл данных только с процентами выживаемости лечения и головастика. 

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
...