r - 'mitml'-package:是否有可能在具有二元结果变量的多级逻辑回归输出中获得 ICC?
问题描述
我正在尝试使用 lme4 包运行一系列多级逻辑回归分析。分析运行良好,没有插补,还为我提供了 ICC,使用以下代码:
main_effects <- glmer(SecurityCC ~ Sens_centr + NE_centr + (Sens_centr + NE_centr | Study_ID), data = data_sens_temp, family = binomial)
summ(main_effects)
但是,在使用“mitml”包进行插补后,我似乎无法再获得 ICC 值,我确实想检查它,因为我想评估执行多级分析是否有用。我的代码如下:
## Impute the whole dataset once
type <- c(-2, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1)
names(type) <- colnames(data_sens_temp)
type
imp <- jomoImpute(data_sens_temp, type = type, n.burn = 5000, n.iter = 100, m = 10, seed = 752)
# Extract list of completed data sets
imp.list <- mitmlComplete(imp, print = "all")
# Model with predictors
# Hypothesis 1
fit.imp.1.A1 <- with(imp.list, glmer(SecurityCC ~ Sens_centr + NE_centr + (Sens_centr + NE_centr | Study_ID), family = binomial))
testEstimates(model = fit.imp.1.A1, extra.pars = TRUE)
当我有一个分类结果变量时,“testEstimates”命令似乎没有输出 ICC,因为我也用一个连续的结果变量对此进行了测试,然后我得到了一个 ICC 值。有谁知道如何通过“mitml”包的估算获得多级逻辑回归的 ICC 值?
相关地,我还想知道是否有可能使我使用“mitml”包得到的输出看起来更友好(例如,看起来类似于带有 glmer 输出的“summ()”命令。“testEstimates”命令确实提供了我有输出,但没有包括模型信息(例如观察和组的数量)。
谢谢!
解决方案
推荐阅读
- airflow - 来自 XCom 值的 Airflow BashOperator 参数
- python - 以交互方式突出显示 igraph 中的节点和边
- node.js - stream-chat-js 包的 Angular 编译错误。没有导出成员“客户”
- react-native - 如何在 React Native 的 TextInput 中禁用退格键?
- javascript - 无法在 Fullcalendar Django 中添加事件
- firebase-realtime-database - recyclerView GridLayoutManager 中的图像显示不正确
- python - 需要拆分字典
- reactjs - 关于 useReducer + Context 的具体问题
- r - 刻度填充手册 - 映射到不同绘图上不同颜色的色标值
- python - 在 C++ 中存储数据,如 python 中的字典