В настоящее время я просматриваю учебник R по анализу медиации. После репликации набора данных и моделирования каждого посредника. Есть этап, на котором нужно вычислить веса регрессии для получения плотностей. В руководстве авторы показывают, как вычислить это для переменной numeri c (которая была смоделирована с использованием линейной регрессии) с помощью функции dnorm
.
fitM1<-glm(M1 ~ A + C, family = gaussian("identity"), data =dat)
num1<- with(extdat1, dnorm(M1, mean = predict(fitM1, newdata = within(extdat, A<- a1),
type = "response"),
sd = sqrt(summary(fitM1)$dispersion)))
denom1 <- with(expdat1,
dnorm(M1, mean = predict(fitM1, newdata = within(expData, A <- a2),
type = "response"),
sd = sqrt(summary(fitM1)$dispersion)))
expdat1$W1 <- num1/denom1
Однако пытаюсь применить это к моим данным , мой посредник - двоичная переменная. Как рассчитать вес регрессии logisti c?