如何解决MICE多次插补后比较多级模型时出错
我使用MICE软件包安装了两个嵌套模型,以进行多次插补。我现在想比较模型,但是函数D3产生以下错误:
as.data.frame.default(data)中的错误: 无法将类“ c(“ lmeStructInt”,“ lmeStruct”,“ modelStruct”)'强制转换为data.frame
数据集包含来自家庭的问卷调查数据。当父母双方都填写有关孩子的调查表时,使用lme拟合了一个多层次模型。我有以下代码:
MLM_analyses <- data.frame(werkdoc$Total_Problems_TScore,werkdoc$Internalizing_Problems_TScore,werkdoc$Externalizing_Problems_TScore,werkdoc$Groep,werkdoc$totaalzichtbaarheid,werkdoc$totaaloverig,werkdoc$PANAS.pos,werkdoc$PANAS.neg,werkdoc$age.child,werkdoc$Geslacht,werkdoc$age.parent,werkdoc$Geslacht.kd,werkdoc$GezinTotaal,werkdoc$opleiding.laag.midden.hoog)
qpredR <- quickpred(MLM_analyses,exclude = c("werkdoc.Total_Problems_TScore","werkdoc.Internalizing_Problems_TScore","werkdoc.Externalizing_Problems_TScore","werkdoc.GezinTotaal","werkdoc.age.child"))
meth <- make.method(MLM_analyses)
meth[c("werkdoc.Total_Problems_TScore","werkdoc.age.child")] <- ""
qpredR[,c("werkdoc.GezinTotaal")] <- 0
diag(qpredR) <- 0; qpredR[which(meth == ""),] <- 0
rowSums(qpredR); mean(rowSums(qpredR[-which(rowSums(qpredR)== 0),]))
ini <- mice(MLM_analyses,predictorMatrix = qpredR,maxit = 0,m = 1,printFlag = F,method = meth,seed = 2020)
ini$loggedEvents
dsImp <- mice(MLM_analyses,maxit = 5,m = 5,seed = 2020)
plot(ini)
#Create the models
MLM.ext.1 <- with(dsImp,expr = lme(werkdoc.Externalizing_Problems_TScore ~ werkdoc.Groep + ns(werkdoc.totaalzichtbaarheid,2) + werkdoc.totaaloverig + werkdoc.PANAS.pos + werkdoc.PANAS.neg + werkdoc.age.child + werkdoc.Geslacht.kd + werkdoc.age.parent + werkdoc.Geslacht + werkdoc.opleiding.laag.midden.hoog,random = ~ 1 | werkdoc.GezinTotaal,na.action = na.omit))
MLM.ext.2 <- with(dsImp,expr = lme(werkdoc.Externalizing_Problems_TScore ~ werkdoc.totaaloverig + werkdoc.Geslacht,na.action = na.omit))
D3(MLM.ext.1,MLM.ext.2)
有人知道这里出了什么问题吗?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。