У меня есть вопрос относительно интерпретации сюжета, созданного с помощью функции межсюжета в R. Цель состоит в том, чтобы выяснить, какое влияние оказывает участие на то, мигрирует ли человек или нет, и как богатство человека (HH_type
) меняет этот эффект.
participation
— непрерывная переменная, HH_type
принимает значение 1 для бедных, 2 для средних и 3 для богатых.
Я оценил эту модель:
library(lme4)
Model_int1 <- glmer(migration ~ participation:HH_type + participation + HH_type + sex + age + education + hh_migrated + (1 | country/district),
data = migration_sudden,
family = binomial(link = "logit"))
и я получил негативное и значительное взаимодействие.
Для создания сюжета я использовал функцию interplot
:
library(interplot)
interplot(m = Model_int1, var1 = "participation", var2 = "HH_type") +
xlab("Wealth") +
ylab("Estimated Coefficient for Participation") +
theme_bw() +
ggtitle("Estimated Coefficient of Participation on migration by household wealth") +
theme(plot.title = element_text(face="bold")) +
geom_hline(yintercept = 0, linetype = "dashed")
и я получаю этот сюжет:
Теперь я не уверен, как интерпретировать этот сюжет. Правильно ли говорить, что бедные люди с большей вероятностью мигрируют, если у них есть широкие возможности для участия; для людей со средним достатком участие не влияет на то, мигрируют они или нет, а богатые люди с меньшей вероятностью мигрируют, если у них больше возможностей для участия?