R中带有ggeffect()的固定效果

问题描述 投票:1回答:1

我想知道是否有人有使用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)

Plot with expected confidence intervals

这为我提供了很好的围绕自变量的干净置信区间。但是,一旦添加状态固定效应,即使仍然存在超强关系,置信区间也将从0扩展到1。

m2 <- glm(depvar ~ indvar + state, data = data, family = "binomial")

margin <- ggeffect(m2, "indvar")

plot(margin)

Plot with super wide confidence intervals

非常感谢!

我想知道是否有人有使用ggeffect()并具有许多固定效果的经验。尽管我的自变量与因变量密切相关,但我得到了这些巨大的信心...

r
1个回答
0
投票

感谢发布。

© www.soinside.com 2019 - 2024. All rights reserved.