패키지의 비 네트를 사용하여 R의 중개 패키지를 둘러 보려고합니다.
mediate()
함수 의 출력을 이해하는 데 어려움을 겪고 있습니다.
require("mediation")
require("sandwich")
data("framing")
med.fit <- lm(emo ~ treat + age + educ + gender + income, data = framing)
out.fit <- glm(cong_mesg ~ emo + treat + age + educ + gender + income,
data = framing, family = binomial("probit"))
summary(out.fit)
# OR for sending a message to a senator for treated compared to untreated.
exp(as.numeric(coef(out.fit)[3]))
# mediation
med.out <- mediate(med.fit, out.fit, treat = "treat", mediator = "emo",
robustSE = TRUE, sims = 100)
summary(med.out)
...
Estimate 95% CI Lower 95% CI Upper p-value
ACME (control) 0.0802 0.0335 0.1300 0.00
ACME (treated) 0.0808 0.0365 0.1311 0.00
ADE (control) 0.0142 -0.1030 0.1325 0.78
ADE (treated) 0.0147 -0.1137 0.1403 0.78
Total Effect 0.0949 -0.0316 0.2129 0.14
Prop. Mediated (control) 0.7621 -2.0926 4.9490 0.14
Prop. Mediated (treated) 0.7842 -1.9272 4.6916 0.14
ACME (average) 0.0805 0.0350 0.1304 0.00
ADE (average) 0.0145 -0.1087 0.1364 0.78
Prop. Mediated (average) 0.7731 -2.0099 4.8203 0.14
...
이것은 치료 효과 중 8.08 %가 치료되는 것들 중에서 정서적 상태를 통해 매개됨을 의미합니까? 아니면 계수의 변화 treat
입니까?
누군가 출력을 설명 할 수 있다면 대단히 감사하겠습니다.