在常规数据集上,调用将如下所示:
library(ROCR)
mymodelFit1 <- glm(data = myData, Outcome~ predictor1+
predictor2+ predictor3,
family = binomial(link = "logit")))
prob <- predict(mymodelFit1, newdata=myData, type="response")
pred <- prediction(prob, myData$Outcome)
perf <- performance(pred, measure = "tpr", x.measure = "fpr")
plot(perf)
auc <- performance(pred, measure = "auc")
auc <- auc@y.values[[1]]
auc
现在,我已经使用鼠标包来填补缺失的数据:
library(mice)
impData <- mice(analysis_set,m=5,maxit=50,meth='pmm',
seed=500, printFlag = FALSE)
使用插补数据集生成逻辑模型:
mymodelFit2 <- with(data = impData, exp = glm(Outcome~ predictor1+
predictor2 + predictor3,
family = binomial(link = "logit")))
问题是:如何计算新模型(mymodelFit2)的AUC?
&#34小鼠&#34;为缺失值生成x个集合(在这种情况下为5,在m =参数mouse()中指定)。该模型具有x数(在这种情况下为5)的估计值。
是否有功能或简单的方法来计算这些模型的AUC?我怎样才能获得汇集模型的AUC?