我有一个数据集，我试图在多水平设置中使用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")

本质上，我试图用可能以2l.开头的东西来代替polyreg。我试过2l.2stage.bin，但那是二进制变量。谢谢你，谢谢