嗨:我正在尝试使用 R 中的 agricolae 包计算某些数据的 LSD。从文档中看似乎很简单，但结果中缺少部分 LSD。
我发现其他一些人也遇到了这个问题，但无法找到他们问题的答案。 https://stat.ethz.ch/pipermail/r-help/2013-December/364391.html
这是我正在使用的代码，以及我生成的一些数据，这些数据复制了我在真实数据上得到的结果。
library(agricolae)
## Create Data Frame
factor <- c(1,1,1,1,2,2,2,2,3,3,3)
var <- c(22,23,24,27,13,17,18,10,31,29,33)
df <- data.frame(factor,var)
## Run LSD test
model <- aov(var ~ factor, df)
out <- LSD.test(model, "factor")
## Results
out$statistics$LSD
NULL
最佳答案
输出中没有给出 LSD 值，因为看起来只有当您拥有完美平衡的数据并且您提供的数据是不平衡的(即 4 组“1”复制，4 组“2”复制，但只有3 组“3”重复)。当数据平衡时，给出LSD值，例如:
library(agricolae)
## Create Data Frame
factor <- c(1,1,1,1,2,2,2,2,3,3,3,3)
var <- c(22,23,24,27,13,17,18,10,31,29,33,27)
df <- data.frame(factor,var)
## Run LSD test
model <- aov(var ~ factor, df)
out <- LSD.test(model, "factor")
## Results
out$statistics$LSD
11.04747
关于r - 使用 Agricolae 包在 R 中没有给出 LSD 值，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/39861633/