Question

我正在尝试分析分割图设计中的数据。我认为我没有足够的经验自己编写Tukey的HSD测试，尽管我已经能够为其他实验设计做到这一点。我已经在堆栈溢出和Google上进行了搜索，但是我无法获得一组有效的代码。

这是我的数据集的头。如果需要，我可以分享更多内容：

 dput(head(rawData))
structure(list(ï..Plot = 2111:2116, Variety = structure(c(5L, 
4L, 3L, 6L, 1L, 2L), .Label = c("Burbank", "Hodag", "Lamoka", 
"Norkotah", "Silverton", "Snowden"), class = "factor"), Rate = c(4L, 
4L, 4L, 4L, 4L, 4L), Rep = c(1L, 1L, 1L, 1L, 1L, 1L), totalTubers = c(594L, 
605L, 656L, 729L, 694L, 548L), totalOzNoCulls = c(2544.18, 2382.07, 
2140.69, 2401.56, 2440.56, 2503.5), totalCWTacNoCulls = c(461.76867, 
432.345705, 388.535235, 435.88314, 442.96164, 454.38525), avgLWratio = c(1.260615419, 
1.287949374, 1.111981583, 1.08647584, 1.350686661, 1.107173509
), Hollow = c(14L, 15L, 22L, 25L, 14L, 13L), Double = c(10L, 
13L, 15L, 22L, 11L, 9L), Knob = c(86L, 80L, 139L, 156L, 77L, 
126L), Researcher = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = "Wang", class = "factor"), 
    CullsPounds = c(1.75, 1.15, 4.7, 1.85, 0.8, 5.55), CullsOz = c(28, 
    18.4, 75.2, 29.6, 12.8, 88.8), totalOz = c(2572.18, 2400.47, 
    2215.89, 2431.16, 2453.36, 2592.3), totalCWTacCulls = c(466.85067, 
    435.685305, 402.184035, 441.25554, 445.28484, 470.50245)), row.names = c(NA, 
6L), class = "data.frame")

我对查看不同费率的不同字母类别感兴趣。我想把它们放在一个图表上，该图表具有每个Variant（总共六个）的条形图，y轴上的totalCWTacCulls，x轴每个上的5 Rates

这是我用于此数据的方差分析：

require(agricole)
model <- with(rawData, sp.plot(Rep, Rate, Variety, totalCWTacCulls))

我正在尝试使用HSD.test函数，但是对于使它起作用的工作我没有足够的了解。这是我目前的尝试：

HSD <- HSD.test(model, "Rate", alpha = 0.05, group=TRUE)
print(HSD$groups)

但是，我得到的答复是我的MSerror丢失了。以前使用此方法分析实验时，我没有遇到这个问题，但这是我对分割图设计进行的第一个分析。

任何帮助或指导将不胜感激。谢谢！

如何在分割图设计上进行Tukey的HSD测试？

0 个答案: