library(readr)
datos_temp <- read_csv("met_atemp.csv")
## Parsed with column specification:
## cols(
## JOBJBMET_ATemp = col_double(),
## JOBJBMET_RH = col_double()
## )
head(datos_temp)
## # A tibble: 6 x 2
## JOBJBMET_ATemp JOBJBMET_RH
## <dbl> <dbl>
## 1 23.7 72
## 2 21.8 76
## 3 22 75
## 4 21.9 75
## 5 21.9 76
## 6 22.1 76
JOBJBMET_ATemp
plot(tstemp)
Figura 1. Serie de tiempo para datos meteorológicos: Temperatura (inicio:2005, fin:2020, frecuencia:50)
tstemp2 <- ts(datos_temp$JOBJBMET_ATemp)
plot(tstemp2)
Figura 2. Serie de tiempo para datos meteorológicos: Temperatura (sin especificaciones)
#descomposición de la ts
fit.datos_temp <- stl(tstemp, s.window = "periodic")
plot(fit.datos_temp)
#gráficas para ver detalles de meses y estacionalidad
library(forecast)
## Registered S3 method overwritten by 'quantmod':
## method from
## as.zoo.data.frame zoo
Figura 3. Descomposición de la serie de tiempo para datos meteorológicos: Temperatura (inicio:2005, fin:2020, frecuencia:50)
monthplot(tstemp, xlab = "", ylab = "")
Figura 4. Serie de tiempo (monthplot) para datos meteorológicos: Temperatura (inicio:2005, fin:2020, frecuencia:50)
seasonplot(tstemp, year.labels = TRUE, main = "", type = "l")
Figura 5. Serie de tiempo (seasonplot) para datos meteorológicos: Temperatura (inicio:2005, fin:2020, frecuencia:50)
JOBJBMET_RH
tshr <- ts(datos_temp$JOBJBMET_RH, start = c(2005), end = c(2020), frequency = 50)
plot(tshr)
Figura 6. Serie de tiempo para datos meteorológicos: Humedad relativa(inicio:2005, fin:2020, frecuencia:50)
tshr2 <- ts(datos_temp$JOBJBMET_RH)
plot(tshr2)
Figura 7. Serie de tiempo para datos meteorológicos: Humedad relativa (sin especificaciones)
#descomposición de la ts
fit.datos_temp <- stl(tshr, s.window = "periodic")
plot(fit.datos_temp)
#gráficas para ver detalles de meses y estacionalidad
library(forecast)
Figura 8. Descomposición de la serie de tiempo para datos meteorológicos: Humedad relativa (inicio:2005, fin:2020, frecuencia:50)
monthplot(tshr, xlab = "", ylab = "")
Figura 9. Serie de tiempo (monthplot) para datos meteorológicos: Humedad relativa (inicio:2005, fin:2020, frecuencia:50)
seasonplot(tshr, year.labels = TRUE, main = "", type = "l")
Figura 10. Serie de tiempo (seasonplot) para datos meteorológicos: Humedad relativa (inicio:2005, fin:2020, frecuencia:50)