r - 可视化多重插补数据的多重密度图

Question

我有一个类的多重插补数据集，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]]))

标签： rmachine-learningplotdata-visualization

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]))))