具有分类/因子变量的多级多重插补 (MICE)？

Question

我有一个数据集，我尝试在包中使用多重插补mice，miceadds并micemd在多级设置中使用分类/因子变量。我可以将该方法2l.2stage.pois用于连续变量，效果非常好。分类变量的插补也适用于polyreg，但这并没有利用多级数据。

这是说明我的问题的可复制代码。

dat <- data.frame(x=rep(c("A", "B"), 50),
                  y=rep(c(10:14), 20),
                  r=rep(c(1:50), 2),
                  z=rep(c(1:4), 25))

dat[40:44, c("x")] <- NA
dat[60:64, c("y")] <- NA
dat[80:84, c("r")] <- NA

predm <- mice::make.predictorMatrix(data=dat)

predm[, c("z")] <- -2
predm[, c("r")] <- 3
predm[row(predm) == col(predm)] <- 0
meth <- c("x"="polyreg", "y"="2l.2stage.pois", "r"="", "z"="")

imputed = mice(dat, method=meth, predictorMatrix=predm, m=1, maxit=1)
imp <- complete(imputed)

table(imp$x, dat$x, useNA = "always")
table(imp$y, dat$y, useNA = "always")

本质上，我试图polyreg用可能应该以 开头的东西来代替2l.。我已经尝试过了2l.2stage.bin，但那是针对二进制变量的。谢谢！

Answer 1

mice目前，在包中或相邻包中没有实现多级多项式回归插补R。解决方法是使用pmm。根据经验，对于非聚类数据中的分类变量效果很好，另请参阅FIMD书中pmm的建议。mice

多级方法在as和aspmm中实现。miceadds2l.pmmmicemd2l.2stage.pmm