在R中使用Agricolae软件包没有给出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

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

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