0

我正在尝试使用 Agricolae 包中的 HSD.test 对我的数据进行 Tukey 测试。但是，我的结果始终只是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")

结合收集和灌溉因素（我想对 2 种收集类型和 5 种灌溉类型的这 10 种总处理方法中的一组进行 Tukey 的处理）：

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

方差分析：

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

图基的：

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

回复： NULL

知道为什么我的回答是空的吗？非常感谢！

4

1 回答 1

0

我遇到了同样的问题，发现治疗变量需要完全按照模型中的方式进行标记。它将其读取为字符串，而不是标签。在您的情况下，您需要trt = "Tukey_data$Tukey_ID". 在另一个例子中，当我定义我的模型时，我使用as.factor(Tukey_ID)并且需要将治疗定义为trt = "as.factor(Tukey_ID)"

于 2021-05-10T18:05:37.300 回答