R中的多個模型組合的ROC曲線和AUC

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值對每個模型的。