运行卡方以获得多重插补数据/从列表列表中提取估计值

Question

我已经运行mice来创建50个插补数据集。现在，我正在尝试计算一些变量的卡方统计量。

可以获得每个数据集的单独估算值：

list <- with(imp, chisq.test(var1, var2, correct = FALSE))

但结果不能汇总。可以将结果组合在miceadds中，但该函数需要卡方统计量的矢量。

如何从list$analyses[[i]]$statistic中提取相关信息？我认为lapply是合适的，但我不确定如何为列表列表指定它。

Answer 1

编辑：提供的答案here似乎是合适的：

x2 <- lapply(list$analyses, `[`, 'statistic')
x2 <- unlist(x2, use.names = FALSE)