【发布时间】:2017-11-23 09:31:27
【问题描述】:
我尝试使用 Influence.ME 库在 glmer 模型 (lme4) 中检测有影响的案例。
我收到以下错误:
data.adapted[, group.var] 中的错误:“闭包”类型的对象是 不可子集
我知道这个错误意味着什么,我正在尝试索引一个函数。问题是我不知道如何解决它。它是由我的数据矩阵的格式引起的吗?因为我估计了一个逻辑回归模型?
感谢任何帮助。
===== 示例代码有效:
data(school23)
model.a <- lmer(math ~ structure + SES + (1 | school.ID), data=school23)
alt.est.a <- influence(model=model.a, group="school.ID")
===== 我的代码产生了错误:
m1 <-glmer(outcome ~ predictor + (1 | subjects), family = binomial)
estex.m1 <- influence(m1,"subjects")
其中结果、预测变量和受试者三者都包含一个包含 21,007 个观察值的向量。 (注意。在 glmer-function 中包含 data 参数没有帮助)
【问题讨论】:
-
只是为了确定您是否尝试过
data <- cbind(outcome, predictor, subjects)和m1 <-glmer(outcome ~ predictor + (1 | subjects), family = binomial, data=data)? -
@abhiieor,感谢您的评论。是的,我做到了。