运行 乘法插补数据的卡方/从列表列表中提取估计值
Running chi-square for multiply-imputed data / extracting estimates from a list of lists
我 运行 mice
可以创建 50 个估算数据集。现在,我正在尝试计算一些变量的卡方统计量。
可以获得每个数据集的单独估计值:
list <- with(imp, chisq.test(var1, var2, correct = FALSE))
但是无法合并结果。可以合并 miceadds
中的结果,但该函数需要一个卡方统计向量。
如何从list$analyses[[i]]$statistic
中提取相关信息?我假设 lapply
是合适的,但我不确定如何为列表列表指定它。
编辑:提供的答案 here 似乎是合适的:
x2 <- lapply(list$analyses, `[`, 'statistic')
x2 <- unlist(x2, use.names = FALSE)
我 运行 mice
可以创建 50 个估算数据集。现在,我正在尝试计算一些变量的卡方统计量。
可以获得每个数据集的单独估计值:
list <- with(imp, chisq.test(var1, var2, correct = FALSE))
但是无法合并结果。可以合并 miceadds
中的结果,但该函数需要一个卡方统计向量。
如何从list$analyses[[i]]$statistic
中提取相关信息?我假设 lapply
是合适的,但我不确定如何为列表列表指定它。
编辑:提供的答案 here 似乎是合适的:
x2 <- lapply(list$analyses, `[`, 'statistic')
x2 <- unlist(x2, use.names = FALSE)