R,mitools::MIcombine,没有p值的原因是什么?
R, mitools::MIcombine, what is the reason for no p-values?
我目前 运行 一个简单的线性回归模型,在 R 中有 5 个多重估算数据集。
例如model <- with(imp, lm(outcome ~ exposure))
为了汇集汇总估计,我可以使用 mitools 包中的命令 summary(mitools::MIcombine(model))
。但是,这不会给出 p 值的结果。我也可以使用 mice 包中的命令 summary(pool(model))
,这确实给出了 p 值的结果。
因此,我想知道 MIcombine 不产生 p 值是否有特定原因?
查看文档后,mitools
库不提供 p-values 似乎没有什么特别的原因。虽然,该软件包的重点是插补,而不是模型结果。
但是,您不需要这些包中的任何一个来查看您的结果——连同每个模型 p-values。我开始将其写为评论,但决定包含代码。如果您不知道...您可以使用基数 R 的 summary
。我意识到 mice
的输出是比较性的,mitools
也是。我认为这也很重要。
如果您的调用输出是 model
,那么这将起作用。
library(tidyverse)
map(1:length(model), ~summary(model[.x]))
我目前 运行 一个简单的线性回归模型,在 R 中有 5 个多重估算数据集。
例如model <- with(imp, lm(outcome ~ exposure))
为了汇集汇总估计,我可以使用 mitools 包中的命令 summary(mitools::MIcombine(model))
。但是,这不会给出 p 值的结果。我也可以使用 mice 包中的命令 summary(pool(model))
,这确实给出了 p 值的结果。
因此,我想知道 MIcombine 不产生 p 值是否有特定原因?
查看文档后,mitools
库不提供 p-values 似乎没有什么特别的原因。虽然,该软件包的重点是插补,而不是模型结果。
但是,您不需要这些包中的任何一个来查看您的结果——连同每个模型 p-values。我开始将其写为评论,但决定包含代码。如果您不知道...您可以使用基数 R 的 summary
。我意识到 mice
的输出是比较性的,mitools
也是。我认为这也很重要。
如果您的调用输出是 model
,那么这将起作用。
library(tidyverse)
map(1:length(model), ~summary(model[.x]))