r - 使用 Agricolae 包的 R 中没有给出 LSD 的值

Question

您好:我正在尝试使用 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

Answer 1

输出中没有给出 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