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

问题描述

我有一个 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]]))

解决方法

您需要 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]))))