r - 警告信息:模型并非都适合相同数量的观测值
问题描述
我正在安装几个具有不同结构的模型。每个结构的示例:
> m1 <- glmer(X~Y1+Y2+Y3+(1|B/C),control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), data = dataset, family=binomial(link="logit"))
>
> m2 <- glmer(X~Y1*A+Y2*A+Y3*A+(1|B/C),control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), data = dataset, family=binomial(link="logit"))
>
> m3 <- glmer(X~(1+Y1|B)+Y2+Y3,control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), data = dataset, family=binomial(link="logit"))
我对所有这些都使用相同的数据集,并且我希望使用包中的AIC()
公式来比较 AIC 值stats
。运行后AIC(m1, m2, m3)
,我收到以下警告消息:
Warning message:
In AIC.default(m1, m2, m3) :
models are not all fitted to the same number of observations
我检查了所有这些观察的数量,它们确实不同:
> nobs(m1)
[1] 157278
> nobs(m2)
[1] 150164
> nobs(m3)
[1] 157278
但我不明白为什么这个值不同,因为我使用的是完全相同的数据集(我从一开始就运行了两次,以及模型,以确保所有模型都相同)。
会不会和模型的公式有关?
解决方案
推荐阅读
- ios - 无法存储 Firebase 数据类型“AuthService”没有成员“setUserinformation”
- google-cloud-platform - GKE - 到 VPN 网络的 Kube-DNS 解析
- python - 如何使我的多个 if 语句起作用?或者最好如何在我的代码中使用嵌套循环
- php - Imagen、HEX 到 Byte 和 Byte 到 HEX 并显示 img。Visual Basic 到 PHP
- r - 如何在每个元素前面附加列号?
- haskell - 尝试使用 if 跟踪函数
- c++ - 如何简化 std::variant 类类型
- sql - 为具有两个字段的不同组合的记录选择行 ID
- istio - 如何获得 Istio 策略正在运行的反馈
- ios - SceneView 相机预览作为背景