Estoy tratando de modelar una especie vegetal (compuesta por 298 ocurrencias) con BIOMOD2 utilizando la función cobertura terrestre consensuada como variables explicativas. Sin embargo, R me muestra un mensaje de error durante el formateo de los datos de entrada que no entiendo en absoluto:
Error en .local(sp, env, ...) : objeto 'mask.tmp' no encontrado
Tengo la sensación de que el problema está relacionado con el área de entrenamiento o algo similar, porque si ejecuto el análisis en otras áreas de estudio (con las mismas variables explicativas, pero cambiando las ocurrencias y el área de entrenamiento), ¡el código funciona perfectamente!
¿Podría alguien ayudarme a solucionar este error y a entender de dónde procede?
El siguiente código es un ejemplo de mi código utilizado:
#Loading libraries
library(biomod2)
library(raster)
library(rgdal)
#Loading occurrences
sic.occ.EUR<- as.data.frame(readOGR(dsn = "C:\\Users\\Invader\\Desktop\\minba1", layer = "occ_JAP"))
sic.occ.EUR <- sic.occ.EUR[,c(4,5)]
sic.occ.EUR$Sycios <- 1
names(sic.occ.EUR) <- c("X_WGS84","Y_WGS84","Sycios")
#Loading explanatory variables
Tbio <- raster("C:\\Users\\Invader\\Desktop\\minba1\\ consensus_full_class_4_crop.tif")
Pbio <- raster("C:\\Users\\Invader\\Desktop\\minba1\\ consensus_full_class_7_crop.tif")
myExpl.EUR <- stack(Tbio,Pbio)
myRespName.EUR <- "Sycios"
myExpl.proj.EUR <- stack(Tbio,Pbio)
final.EUR.myResp <- as.numeric(sic.occ.EUR[,myRespName.EUR])
myRespXY.EUR <- sic.occ.EUR[,c("X_WGS84","Y_WGS84")]
#Formating input data
myBioData.EUR <- BIOMOD_FormatingData(resp.var=final.EUR.myResp,
expl.var=myExpl.EUR,
resp.xy=myRespXY.EUR,
resp.name=myRespName.EUR,
PA.nb.rep = 5,
PA.nb.absences = 10000,
PA.strategy ='random',
PA.dist.min = NULL,
PA.dist.max = NULL,
PA.sre.quant = 0.025,
PA.table = NULL,
na.rm = TRUE)
