首页 > 解决方案 > '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”命令确实提供了我有输出,但没有包括模型信息(例如观察和组的数量)。

谢谢!

标签: rlogistic-regressionimputationmulti-levelmultilevel-analysis

解决方案


推荐阅读