我有一个在R中运行鼠标算法得到的中间对象。我正在尝试使用plot.mids来评估收敛。我希望能够为两个指定的变量绘制收敛诊断。这是一个例子:
library(mice)
nhanes1 <- nhanes
nhanes1[2,"age"] <- NA
mi1 <- mice(nhanes1,m=3,maxit=8,seed=2)
plot(mi1)
这里,只显示变量chl的图表,虽然您可以点击后退箭头查看其他图表,但我想在报告的同一图表中一起显示chl和bmi。我尝试通过y =公式指定变量,但不认为我理解如何使用它。
plot(mi1, chl~bmi)
这似乎试图将chl的手段与bmi(我认为)的手段联系起来,这不是我所追求的。