1) Identificar una especie y buscar literatura sobre las condiciones del clima aptas para su presencia

2)Identificar un pai representativo en el cual se conozca que vive la especie y Filtrar el shape global con este país

3) Cortar las imagenes de clima para el país

4) Crear un mapa de aptitud de clima por medio de los humbrales identificados

5) Interpretación de los mapas de aptitud

6) Calcular el % de áreas que son aptas en el país (Todos los meses aptos)

Exploración con la variable “Temperatura”

Dynastes satanas

Especie de escarabajo que vive en alturas desde los 900 hasta los 2000 metros. Las características climatológicas de su hábitat son temperaturas entre 7 y 24 ºC, con una precipitación media anual de aproximadamente 1500 a 1600 mm.

Bolivia

require(raster)
## Loading required package: raster
## Warning: package 'raster' was built under R version 4.0.5
## Loading required package: sp
## Warning: package 'sp' was built under R version 4.0.5
Map_global<- shapefile("shape_global/g2008_0.shp") # Archivo que almacena las entidades geométricas del mapa con el que se trabajará.
Map_global
## class       : SpatialPolygonsDataFrame 
## features    : 534 
## extent      : -180, 180, -55.72333, 83.62742  (xmin, xmax, ymin, ymax)
## crs         : +proj=longlat +datum=WGS84 +no_defs 
## variables   : 11
## names       :             AREA,       PERIMETER, G2008_0_, G2008_0_ID, ADM0_CODE,   ADM0_NAME, LAST_UPDAT, CONTINENT,                    REGION, STR_YEAR0, EXP_YEAR0 
## min values  : 0.10069705556089,   1.40729013556,    10019,          1,         1, Afghanistan,   20050415,    Africa, Australia and New Zealand,         0,         0 
## max values  :  2825.6935742546, 969.26620675757,     9999,       9998,        98,    Zimbabwe,   20081118,   Oceania,            Western Europe,      2006,         0
plot(Map_global) #Gráfica del mapa 

Tabla<- Map_global@data #Creación de una tabla con toda la información que contiene el mapa.
View(Tabla)
Selec<- which(Map_global$ADM0_NAME=="Bolivia") #Selección del país de interes 
Bolv<- Map_global[Selec,] 
plot(Bolv)

Temperatura en Bolivia

Temp_maps<- list.files("C:/Users/sebas/Desktop/Angie Herrera/2021-1/Bioestadistica/Mapas/Actividad temperatura/", full.names = TRUE) #Cargar los mapas de los meses del año con variable temperatura

Temp_global<- stack(Temp_maps) #Transforma los datos disponibles como columnas separadas en un marco de datos o una lista en una sola columna para que puedan utilizarse en los análisis

names(Temp_global)= month.name #Pone los meses del año en los respectivos mapas
Temp_Bolv<- crop(Temp_global, Bolv) #Para cortar el país y unirlo con la temperatura
Temp_Bolv<-mask(Temp_Bolv,Bolv)
plot(Temp_Bolv)

Aptitud de Dynastes satanas en Bolivia

Nicho_DSatanas<- Temp_Bolv>7&Temp_Bolv<24
plot(Nicho_DSatanas)

Por medio de las gráficas podemos observar que el rango geográfico de distribución del escarabajo en Bolivia, según la temperatura,es a lo largo de la cordillera de los Andes. Sin embargo, aunque el rango cambie un poco en algunos meses el año, siempre se tomará una parte de los Andes Bolivianos lo cual, puede indicarnos que en esa región hay mayor densidad poblacional.

Porcentaje de aptitud

cell_size<-area(Nicho_DSatanas, na.rm=TRUE, weights=FALSE)
SinNA<-is.na(cell_size)
Vector<-as.vector(SinNA)
Tabla<-table(Vector)
prop.table(Tabla)
## Vector
##     FALSE      TRUE 
## 0.5729149 0.4270851

El 57% del país tiene las condiciones aptas para que Dynastes satanas habite en él.

Exploración con la variable “Precipitación”

Precp_maps<- list.files("C:/Users/sebas/Desktop/Angie Herrera/2021-1/Bioestadistica/Mapas/Actividad precipitación/", full.names = TRUE)
Precp_global<- stack(Precp_maps) 
names(Precp_global)= month.name 
Precp_Bolv<- crop(Precp_global, Bolv)
Precp_Bolv<-mask(Precp_Bolv,Bolv)
plot(Precp_Bolv)

Nicho2_DSatanas<- Precp_Bolv>1500&Precp_Bolv<1600
plot(Nicho2_DSatanas)