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