我想知道是否有人有使用ggeffect()并具有许多固定效果的经验。尽管我的自变量与因变量密切相关,但是当添加状态固定效应时,我会得到这些巨大的置信区间。]
一个玩具示例:
set.seed(200)
indvar <- runif(500, min = 0, max = .5)
state <- as.factor(rep(c(1:50), 10))
statev <- as.integer(state) * runif(500, 0, 0.02)
depvar <- round(indvar + statev)
data <- data.frame(indvar, state, depvar)
m1 <- glm(depvar ~ indvar, data = data, family = "binomial")
margin <- ggeffect(m1, "indvar")
plot(margin)
这为我提供了很好的围绕自变量的干净置信区间。但是,一旦添加状态固定效应,即使仍然存在超强关系,置信区间也将从0扩展到1。
m2 <- glm(depvar ~ indvar + state, data = data, family = "binomial")
margin <- ggeffect(m2, "indvar")
plot(margin)
非常感谢!
我想知道是否有人有使用ggeffect()并具有许多固定效果的经验。尽管我的自变量与因变量密切相关,但我得到了这些巨大的信心...
感谢发布。