可视化多重估算数据的多重密度图

Question

我有一个 class mids 的多重估算数据集，我想为每个估算变量创建密度图。出于某种原因，我无法使代码正常工作。我不断收到一条错误消息：

Error in `[.data.frame`(r, , xvar) : undefined columns selected

这是我试过的：

library(mice)
imp <- mice(airquality, seed = 123, 3, print = F)
vars <- c('Ozone', 'Solar.R')
map(1:length(vars), function(x) densityplot(imp, ~vars[[x]]))

Answer 1

您需要 densityplot 中的公式。尝试-

library(mice)

imp <- mice(airquality, seed = 123, 3, print = F)
vars <- c('Ozone', 'Solar.R')

list_plots <- purrr::map(seq_along(vars), function(x) 
                 densityplot(imp, as.formula(paste0('~', vars[x]))))

可视化多重估算数据的多重密度图

Visualize multiple density plots for multiply imputed data

plot

r

data-visualization

machine-learning