Нет значения для LSD в R с использованием пакета Agricolae

Question

Привет: Я пытаюсь рассчитать LSD для некоторых данных при использовании пакета agricolae в R. Это кажется достаточно простым из документации, но часть ЛСД отсутствует в результаты.

Я обнаружил, что некоторые другие имели эту проблему, но не смогли найти ответ на свои вопросы. 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