Почему мой тест Тьюки (от Agricolae) не имеет значения? - PullRequest
Новая
       95

Почему мой тест Тьюки (от Agricolae) не имеет значения?

0 голосов
/ 21 июня 2019

Я пытаюсь использовать HSD.test из пакета Agricolae, чтобы выполнить тест Тьюки на моих данных.Тем не менее, мой результат всегда просто слово NULL.

Мои данные: 

structure(list(Sample = c(1011L, 1012L, 1014L, 1024L, 1025L, 
    1026L), Collection = c(1L, 1L, 1L, 1L, 1L, 1L), Irrigation = structure(c(3L, 
    3L, 3L, 5L, 5L, 5L), .Label = c("Rate1", "Rate2", "Rate3", "Rate4", 
    "Rate5"), class = "factor"), Variety = structure(c(2L, 1L, 3L, 
    3L, 2L, 1L), .Label = c("Hodag", "Lamoka", "Snowden"), class = "factor"), 
    Suc = c(0.7333, 0.4717, 0.5883, 0.6783, 0.8283, 0.6833), 
    Gluc = c(0.03, 0.04, 0.043, 0.075, 0.057, 0.087), L = c(59.48, 
    57.59, 59.25, 66.45, 68.29, 65.65), a = c(4.36, 6.85, 3.43, 
    1.7, 0.78, 2.84), b = c(26.82, 27.6, 26.2, 26.14, 25.37, 
    27.19), NoDefect = c(100L, 100L, 100L, 92L, 100L, 100L), 
    Defect = c(0L, 0L, 0L, 8L, 0L, 0L)), row.names = c(NA, 6L
), class = "data.frame")

Объединение факторов сбора и орошения (я хочу сделать Тьюки для группы из этих 10 всегообработка 2 типов сбора и 5 типов орошения): 

Tukey_data <- data %>%
      mutate(Tukey_ID = paste(Collection, Irrigation, sep="_"))

Anova: 

Tukey_data$Tukey_ID <- as.factor(Tukey_data$Tukey_ID) 
TukeyAov <- aov(Suc ~ Tukey_data$Tukey_ID, data = Tukey_data)

Tukey's: 

HSD.test(TukeyAov, "Tukey_ID", group=FALSE)

Ответ: NULL

Есть идеи, почему мой ответ нулевой?Большое спасибо!

...