Невозможно получить уточненные средства для glmer с помощью lsmeans

голоса
0

У меня есть GLM, что я хотел бы получить скорректированные средства для использования lsmeans. Следующий код делает модель (и, кажется, делает это правильно):

library(lmerTest)
data$group <- as.factor(data$grp)
data$site <- as.factor(data$site)
data$stimulus <- as.factor(data$stimulus)

data.acc1 = glmer(accuracy ~ site + grp*stimulus + (1|ID), data=data, family=binomial)

Однако, используя при попытке использовать любой из кода ниже, чтобы получить уточненные средства для модели, я получаю ошибку

Ошибка в lsmeansLT (модели, test.effs = test.effs, DDF = DDF):
Модель не линейная модель смешанных эффектов.

lsmeans(data.acc1, stimulus)

или

data.lsm <- lsmeans(data.acc1, accuracy ~ stimulus ~ grp)
pairs(data.lsm)

Любые suggestiongs?

Задан 20/03/2017 в 13:10
источник пользователем
На других языках...                            


1 ответов

голоса
1

Проблема заключается в том, что вы создали обобщенную линейную смешанную модель с использованием glmer()(в данном случае смешанная модель логистической регрессии) не является линейной смешанной моделью с использованием lmer(). lsmeans()Функция не принимает объекты , созданные , glmer()потому что они не являются линейными смешанными моделями.

Ответы на этот пост может помочь: я не могу получить выход lsmeans в glmer

И этот пост может быть полезно , если вы хотите , чтобы понять / вычислить предельные эффекты для смешанного GLMS: Есть ли способ получить «маргинальные эффекты» от `glmer` объекта

Ответил 20/03/2017 в 16:01
источник пользователем

Cookies help us deliver our services. By using our services, you agree to our use of cookies. Learn more