Как включить только один контраст категориального предиктора уровня 3+ в структуру случайных эффектов модели Лмера? - PullRequest
Новая
       108

Как включить только один контраст категориального предиктора уровня 3+ в структуру случайных эффектов модели Лмера?

0 голосов
/ 19 сентября 2019

Некоторые из смешанных моделей, которые я запускал в прошлом, имели категориальные предикторы с более чем двумя уровнями.Когда я включаю такие переменные в структуру случайных эффектов модели lmer(), я получаю результаты, которые предполагают, что существует случайный эффект одного контраста (или сравнения), но нет никакого случайного эффекта остальных контрастов.

Мой главный вопрос : Как, если возможно, я могу указать только один контраст для включения в качестве случайного эффекта с категориальным предиктором, имеющим 3 или более уровней?

Пример : скажем, я собрал данные о 20 лицах, каждый из которых прошел 3 условия (A, B, C).Из предыдущей работы я знаю, что эффект перехода от условия A к условию B является фиксированным (без случайности), но я хочу включить случайный эффект для перехода от условия A к C. 

df <- data.frame(PID = rep(1:20, 10), group = rep(c("A", "B", "C"),200), y = rnorm(200))

model <- lmer(y ~ group + (1 + group | PID),
              data = df,
              REML = TRUE)

#Partial summary
summary(model)
Random effects:
 Groups   Name        Variance  Std.Dev.  Corr     
 PID      (Intercept) 6.298e-02 0.2509490          
          groupB      1.907e-08 0.0001381 0.91     
          groupC      4.423e-08 0.0002103 0.75 0.96
 Residual             9.007e-01 0.9490689

Как можнобыть замеченным выше, включая group, поскольку случайный эффект приводит к двум случайным эффектам.Один для сравнения между группой A и группой B, а другой - для сравнения между группой B и группой C. Существует ли способ указать только сравнение между группой A и группой C в синтаксисе случайного эффекта?

...