我正在尝试使用 emmean/lsmeans 进行事后测试,以进行组的多重比较。事实上,我正在尝试重新创建一个在这里找到的示例:https ://rcompanion.org/handbook/G_11.html 。
因此,我有一个看起来像这样的模型:
library(ordinal)
model = clm(Likert.f ~ Instructor + Question + Instructor:Question,
data = Data,
threshold="symmetric")
library(car)
library(RVAideMemoire)
Anova(model,
type = "II")
然后我正在尝试探索交互效果的组分离:
library(lsmeans)
marginal = lsmeans(model,
pairwise ~ Instructor + Question,
adjust="tukey") ### Tukey-adjusted comparisons
marginal
cld(marginal,
alpha = 0.05,
Letters = letters, ### Use lower-case letters for .group
adjust = "tukey") ### Tukey-adjusted comparisons
但我收到以下错误:
Error in UseMethod("vcoc"): no applicable methos for 'vcov' applied to an object of class "clm"
我也对 emmeans 进行了同样的尝试,但出现了同样的错误。
有人可以指导我出了什么问题吗?:)