在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