问题描述
我想将 densityplot()
包中的 mice
输出(densityplot
对象的 mids
方法)打印到 PDF 上的单独页面。对于下面的 MWE,我想在每一页上打印一个图(即第 1 页上的 bmi,然后第 2 页上的 hyp,然后第 3 页上的 chi)。
library(mice)
imp <- mice(nhanes)
densityplot(imp)
编辑:
我稍微修改的解决方案(以编程方式执行此操作):
# Which variables had missing observations
imputed_vars <- names(d_imp$method)[d_imp$method != ""]
pdf(file = "./MICE Density Plots.pdf")
# Plot each imputed variable separately
for(var in imputed_vars){
densityplot(d_imp,formula(paste("~",var)))
}
dev.off()
解决方法
你的意思是这样的:
imp <- mice(nhanes)
densityplot(imp)
densityplot(imp,~bmi)
densityplot(imp,~hyp)
densityplot(imp,~chl)
BR