Я использовал примерный кодер, найденный в https://rdrr.io/cran/agricolae/man/waller.test.html, чтобы запустить waller.test из пакета Agricola R. Если я запускаю код со встроенным набором данных, он работает нормально. Когда я запускаю его с моим набором данных, настроенным на тот же самый фрейм данных с двумя столбцами в первом столбце Фактор с 4 коэффициентами и во втором столбце числовых значений для измеренных значений, я получаю следующую ошибку:
Error in if ((K - IN0/ID0) * (K - IN1/ID1) <= 0) b0 <- t :
missing value where TRUE/FALSE needed
Я пыталсяубедитесь, что набор данных отформатирован аналогично. фрейм данных с Фактором (lifeStage) в первом столбце и Числовым значением (Дыхание) во втором столбце. Когда я запускаю встроенный небольшой файл данных (sweetpotatoe), он работает нормально, мой набор данных большой, 8000 наблюдений останавливается в расчете waller.test с ошибкой выше.
mydata3col <-gather(data.df,Stage,Respiration,worker,drone, pupaD,pupaW,factor_key=TRUE)
str(mydata3col)
mydata2col <- mydata3col[-c(1) ]
model_mydata<-aov(Respiration~Stage, data=mydata2col)
out_mydata <- waller.test(model_mydata,"Stage", group=TRUE)