相关问题:link
我有如下数据:
data <- structure(list(Choice = c(0,1,0),gender = structure(c(1,2,1),label = NA_character_,class = c("labelled","numeric")),uniqueID = structure(c(1,3,4,5,6,7,8,9,9),"numeric"))),row.names = c(NA,-50L),class = c("data.table","data.frame"))
我想用这些数据运行逻辑回归 (glm),但我想控制 uniqueID
(面板数据样式)。
最简单的方法是手动执行此操作,只需运行以下回归:
result <- glm(Choice ~ as.factor(gender)+ as.factor(uniqueID),family = binomial(link=logit),data=data).
问题是完整的数据有太多的 uniqueID
,因此 summary
变得不可读。寻找答案,我找到了this post,建议pglm
。
但是,pglm
显然不支持 within
。我链接的帖子解释了原因。
然后我尝试了 lme4::glmer
,但我无法从文档中真正弄清楚语法应该是什么。
谁能帮助我使手册可读(通过从输出中排除 as.factor(uniqueID)),或者帮助我了解 glmer
的语法?