Многократное сравнение групп с использованием clm: Ошибка в UseMethod ("vco c"): нет применимых методов для 'vcov', примененного к объекту класса "clm"

Question

Я пытаюсь сделать пост-хо c тесты с 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, но с та же ошибка.

Может кто-нибудь подсказать мне, что не так? :)