######## https://spei.csic.es/
setwd("G:/UAAAN/MATERIAS/2021/OHF") # Ajuste de mi directorio
library(readxl) # Libreria usada para caragr datos de xls
Datos = read_excel("Pp_Tmax_Tmin.xlsx", # leyendo el archivo
sheet = "media") # Nombre de la hoja de excel
View(Datos) # Viendo los datos en una tabla

head(Datos) # Ver primeros datos

attach(Datos) # para elegir las variables directamente

#install.packages(SPEI, dependencies = TRUE) # Instalar la libreria
library(SPEI) # cargar la libreria para usarla

PET according to hargreaves # Calculo de la ETP
Datos$ETP_har=hargreaves(Tmin, Tmax,lat=25.443726)

View(Datos) # Viendo los datos en la tabla

attach(Datos) # para elegir las variables directamente

Balance # haciendo el BALANCE PpETP
Datos$BAL = PpETP_har

One and tvelwemonths SPEI
spei1 = spei(Datos[,'BAL'], 1) # Obteniendo el SPE, un mes
plot(spei1) # Haciendo la figura de SPEI

spei12 = spei(Datos[,'BAL'], 12) # Obteniendo el SPE, doce meses
plot(spei12) # Haciendo la figura de SPEI

########### VER LOS ESTADISTICOS DE AJUSTE
class(spei1)

Extract information from spei object: summary, call function, fitted values, and coefficients
summary(spei1)
names(spei1)
spei1$call
spei1$fitted
spei1$coefficients

jpeg(filename="SPEI_Arteaga.jpeg", # Nombre del archivo y extension
width = 33, # Anchura
height = 19, # Altura
res= 300, # Resolucion 72ppi es un estandar
units = "cm") # Unidades.
plot(spei12)
#despues de correr estas lienas, vaya a las lineas donde se genera la plot, y despues
corra la sig linea para cerrar la plot

dev.off() # Cierre del archivo