具有类别/因子变量的多级多重插补（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，但这是针对二进制变量的。谢谢！