Построение эффекта взаимодействия в R в многоуровневом анализе

Question

Я хочу построить эффект взаимодействия с помощью R. Ниже приведена картинка, на которой он должен выглядеть. Я использую многоуровневый анализ с пакетом lme4.

Моя модель выглядит следующим образом:

h3a_c_3 <- lmer (PSS_mean ~ WLoad_mean_gmc * PCT_mean_gmc + sex_m + age_m + startup_m + (1 | код команды), данные = данные) </p>

PSS_mean обозначает Perceived Stress (зависимая переменная), PCT_mean_gmc - это команда PsyCap (Модератор), а WLoad_mean_gmc (Predictor) - рабочая нагрузка.

Я уже пробовал несколько вещей, но пока ничего не получалось.

Большое спасибо заранее!

Сара

https://drive.google.com/open?id=15JjC-s2OqDMkfl9BHMNgXnpkpUlazpfO

Answer 1

Я не могу сказать наверняка без каких-либо данных или примеров для тестирования, но я думаю, что это должно работать:

library(ggeffects)
x <- ggpredict(h3a_c_3, c("WLoad_mean_gmc", "PCT_mean_gmc"))

x
plot(x)

Вы можете найти исчерпывающую документацию для пакета ggeffects здесь: http://strengejacke.github.io/ggeffects

Существует даже пример для смешанных моделей с условиями взаимодействия здесь .