Попытка запустить код
cld(posthoc, alpha=.05, Letters=letters)
для получения буквенных значений для моих данных, но выдает эту ошибку:
Ошибка в UseMethod ("cld"): нет подходящего метода для 'cld 'применяется к объекту класса "c (' emm_list ',' list ')"
arcdata
Код следующий
install.packages("car")
arcdata<- read.csv(file.choose())
arcdata
dput(arcdata)
str(arcdata)
library(ggpubr)
plot<- ggline(arcdata, x="preg", y="arc", color = "herd",
add =c("mean_se", "dotplot"), palette= c("#00AFBB", "#E7B800"))
plot
model1<-lm(arc~herd*preg, data = arcdata)
anova(model1)
plot(model1, 2)
aov_residuals<-residuals(object = model1)
shapiro.test(aov_residuals)
plot(model1, 1)
leveneTest(arc~herd*preg, data= arcdata)
install.packages("lsmeans")
install.packages("multcompView")
library("lsmeans")
library("multcompView")
posthoc<-lsmeans(model1, pairwise~herd*preg, adjust="tukey")
posthoc
cld(posthoc,
alpha=.05,
Letters=letters)
Последний (cld) раздел кода не производит ожидаемые группы букв, возможно, это как-то связано с альфа = .05?