R 在物种分布建模中,如何获得每个模型和集合模型的AUC分数

R 在物种分布建模中,如何获得每个模型和集合模型的AUC分数,r,auc,R,Auc,我想得到评估指标,即每个模型的AUC分数,以及SDM的集成模型的AUC分数。使用R中的BIOMOD2软件包进行了分布建模。通过get_评估,我得到了测试值、数据、截止值、敏感性和特异性,但我希望每个模型的表中有AUC分数,以便我可以有一个箱线图。我知道command models_scores_graph提供了一个图表,但我希望它是表格形式。代码直至模型评分图如下所示: DataSpecies <- read.csv(system.file("external/species/m

我想得到评估指标,即每个模型的AUC分数,以及SDM的集成模型的AUC分数。使用R中的BIOMOD2软件包进行了分布建模。通过get_评估,我得到了测试值、数据、截止值、敏感性和特异性,但我希望每个模型的表中有AUC分数,以便我可以有一个箱线图。我知道command models_scores_graph提供了一个图表,但我希望它是表格形式。代码直至模型评分图如下所示:

DataSpecies <- read.csv(system.file("external/species/mammals_table.csv",
                            package="biomod2"), row.names = 1)



myRespName <- 'GuloGulo'

myResp <- as.numeric(DataSpecies[,myRespName])

myRespXY <- DataSpecies[,c("X_WGS84","Y_WGS84")]

myExpl = raster::stack( system.file( "external/bioclim/current/bio3.grd", 
                             package="biomod2"),
                system.file( "external/bioclim/current/bio4.grd", 
                             package="biomod2"), 
                system.file( "external/bioclim/current/bio7.grd", 
                             package="biomod2"),  
                system.file( "external/bioclim/current/bio11.grd", 
                             package="biomod2"), 
                system.file( "external/bioclim/current/bio12.grd", 
                             package="biomod2"))



 myBiomodData <- BIOMOD_FormatingData(resp.var = myResp,
                             expl.var = myExpl,
                             resp.xy = myRespXY,
                             resp.name = myRespName)

myBiomodOption <- BIOMOD_ModelingOptions()

myBiomodModelOut <- BIOMOD_Modeling( myBiomodData, 
                             models = c('SRE','CTA','RF'), 
                             models.options = myBiomodOption, 
                             NbRunEval=1, 
                             DataSplit=80, 
                             Yweights=NULL, 
                             VarImport=3, 
                             models.eval.meth = c('TSS'),
                             SaveObj = TRUE,
                             rescal.all.models = FALSE,
                             do.full.models = FALSE,
                             modeling.id='test')

myBiomodModelOut_scores <- get_evaluations(myBiomodModelOut)

models_scores_graph(myBiomodModelOut, by = "models" ,
                   metrics = c("ROC","TSS"))

DataSpecies要获得AUC输出,需要在
BIOMOD\u Modeling()
中将其指定为评估方法。正如您所拥有的,您只使用TSS,因此它仅显示为
Testing.data

您可以完全删除
models.eval.meth
参数(这将显示Kappa、TSS和AUC),也可以指定所需的方法(如下添加ROC)

myBiomodModelOut
myBiomodModelOut <- BIOMOD_Modeling( myBiomodData, 
                                     models = c('SRE','CTA','RF'), 
                                     models.options = myBiomodOption, 
                                     NbRunEval=1, 
                                     DataSplit=80, 
                                     Yweights=NULL, 
                                     VarImport=3, 
                                     models.eval.meth = c('ROC', 'TSS'),
                                     SaveObj = TRUE,
                                     rescal.all.models = FALSE,
                                     do.full.models = FALSE,
                                     modeling.id='test')