Question

dredge results snapshot我在MuMin包中使用挖泥功能，使用逻辑回归得出多因子组合。

library(MuMIn)
AA <- dredge(Forward.model,evaluate = TRUE,rank = "AIC",m.lim = c(8,13),trace=2)
write.xlsx(AA,File location)

这导致大约10000种不同的模型组合。现在我想计算所有这些组合的AUC。 R中有包可以做到吗？我不太熟悉在R中编写循环。

Answer 1

我会按照以下步骤操作：

1） 对于每个模型，获取每个样本的预测值。设pred.mat矩阵和预测，其中每行代表一个人，每列代表不同的模型，y带有二分法响应变量的矢量。

2） 使用以下方法为每个模型计算AUC ：

    # Load pROC library (install it if not)
      library(pROC)
    # Compute the AUC values
      auc.val <- apply(pred.mat,2,function(x)roc(pred=x,resp=y)$auc)

因此，auc.val将显示每个模型的AUC值。

R

1 个答案: