没有给予LSD在R中使用农业包装的价值

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

Answer 2

在建立数据框架之后，在建立模型之前，尝试将因子向量作为因子。以下是完整的代码：

library(agricolae)
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) 
df$factor<-as.factor(df$factor) #make the factor "factor"
model<-aov(var~factor, data=df)
(LSD.test(model, "factor"))$groups