我正在为一组物种生成生态位模型,我想使用AUC作为生态位生态质量的指标。开发Maxent的Steven Phillips在他的Maxent手册中提供了用于计算R中AUC的代码。但是,我正在阅读报告部分AUC比率作为更健壮和概念上合理的度量的论文。我想我理解如何使用ROCR R软件包计算部分AUC,但是如何计算AUC比率?
以下是Phillips的教程脚本:
presence<-read.csv("bradypus_variegatus_samplePredictions.csv")
background<-read.csv("bradypus_variegatus_backgroundPredictions.csv")
pp<-presence$Logistic.prediction
testpp<-pp[presence$Test.or.train=="test"]
trainpp<-pp[presence$Test.or.train=="train"]
bb<-background$logistic
combined<-c(testpp,bb)
label<-c(rep(1,length(testpp)),rep(0,length(bb)))
pred<-prediction(combined,label)
perf<-performance(pred,"tpr","fpr")
plot(perf,colorize=TRUE)
performance(pred,"auc")@y.values[[1]] #RETURNS AUC
AUC<-function(p,ind){
pres<-p[ind]
combined<-c(pres,bb)
label<-c(rep(1,length(pres)),rep(0,length(bb)))
predic<-prediction(combined,label)
return(performance(predic,'auc')@y.values[[1]])
}
b1<-boot(testpp,AUC,100) #RETURNS AUC WITH STANDARD ERROR
b1
非常感谢任何建议或意见!谢谢。
答案 0 :(得分:2)
在不知道数据集和应用程序的细节的情况下,
...洙
答案 1 :(得分:1)
包ROCR可以使用fpr.stop=
参数计算部分AUC值。正如约翰所说,这个比率只是这个值除以参考模型的相同计算。