Скачать с ютуб SPEI. Cálculo del Índice Estandarizado de Precipitación-Evapotranspiración в хорошем качестве

SPEI. Cálculo del Índice Estandarizado de Precipitación-Evapotranspiración

######## 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 Pp-ETP Datos$BAL = Pp-ETP_har One and tvelwe-months 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