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