1

我正在使用mice包来获取完整的数据。我认为问题在于我没有对所有数据进行插补,因此其中一些数据具有 NA。(一些数据缺失的变量仅用于预测其他变量的缺失,所以我不想估算这些变量。

我可以用这段代码复制问题:

require(mice)
impute <- mice(
    nhanes, 
    imputationMethod = c(
        "",        # age
        "pmm",     # bmi
        "pmm",  # hyp
        ""         # chl
    ),
    seed = 101)
x11()
densityplot(impute)

Error in density.default(x = c(NA_real_, NA_real_, NA_real_, NA_real_,  : 
  need at least 2 points to select a bandwidth automatically

我怎样才能得到密度图?如果我"""pmm"for替换chl或者只是运行,impute <- mice(nhanes)那么它将与这个例子一起工作,产生这个: 在此处输入图像描述

但我不能用我自己的数据做到这一点,所以我正在寻找另一种方法....只是为了得到 和 的密度图bmihyp在上面运行mice我的代码之后,它不会估算值chl

编辑: 我知道我可以使用我之前的问题的答案中的方法使用ggplot,但在这种情况下,我真的需要使用densityplot

4

1 回答 1

2

'chl' 列中仍然会有大量缺失数据。您可以使用鼠标功能complete将原始非缺失值和估算值组合在一起。但是,如果您将数据和公式的角色颠倒过来(最终)在 hte 帮助页面的详细信息中描述,则 densitplot.mids 函数将执行此操作。

densityplot( x=impute , data= ~ bmi+hyp)

在此处输入图像描述

于 2012-09-12T17:55:34.180 回答