Práctica 2

Juan Carlos Senent Núñez

Ejercicio 1

Realiza un breve análisis descriptivo numérico y gráfico. Convierte los datos en formato de serie temporal y dibuja la serie.

##Lectura de los datos##
##df<- read.csv("unemployment.csv")
Tabla<-scan("https://www.r-exercises.com/wp-content/uploads/2017/04/unemployment.csv", skip=1)
##Crear la serie temporal##
unem<-ts(Tabla,start = 2012, freq=12)
##Dibujar dicha serie##
plot(unem)

Ejercicio 2

Usa la función ses del paquete forecast para obtener la predicción basada en un suavizado esponencial simple para los próximos 12 meses, y dibuja esa predicción.

library(forecast)
##suavizado exponencial simple para los proximos 12 meses##
forecast_ses<-ses(unem, h =12)
##dibujar la predicción anterior##
plot(forecast_ses)

Ejercicio 3

Estima un modelo de suavizado exponencial usando la función ets con los parámetros por defecto. Luego pasa el modelo como input para una predicción (usa la funcion textsl{forecast}) de lso próximos 12 meses, y dibuja la predicción.

##modelo suavizado con ets##
fit<-ets(unem)
##modelo con input para una predicción##
forecast_ets<-forecast(fit,h=12)
##dibujando predicción##
plot(forecast_ets)

Ejercicio 4

Realiza un “summary” del modelo estimado en el ejercicio anterior y encuentra la estructura de la estimación. ¿Incluye esta los componentes de tendcia y estacionalidad? Si estas componentes están presentes, ¿el modelo es aditivo o multiplicativo?

##sumario del modelo suaviado con ets##
summary(fit)

## ETS(M,A,M) 
## 
## Call:
##  ets(y = unem) 
## 
##   Smoothing parameters:
##     alpha = 0.5717 
##     beta  = 1e-04 
##     gamma = 1e-04 
## 
##   Initial states:
##     l = 8.4353 
##     b = -0.0564 
##     s=0.9567 0.9453 0.957 0.9669 1.0245 1.0591
##            1.0365 0.9642 0.9403 1.0224 1.0542 1.0729
## 
##   sigma:  0.021
## 
##      AIC     AICc      BIC 
## 37.38299 50.98299 73.81628 
## 
## Training set error measures:
##                        ME      RMSE      MAE        MPE     MAPE      MASE
## Training set -0.009791689 0.1271141 0.102154 -0.1209349 1.687472 0.1322298
##                    ACF1
## Training set 0.08649403

Respuesta El modelo tiene tendencia y estacionalidad, por otro lado el comportamiento del gráfico corresponde a un modelo Aditivo.

Ejercicio 5

Usa la función ets para estimar dos modelos, uno con tendencia y otro sin capturar la tendencia de esta serie. Realiza la predicción para los próximos 12 meses y dibuja las predicciones de las estimaciones.

##Con tendencia##
fit_1 <- ets(unem, damped = TRUE)
fcast_ets_1 <- forecast(fit_1, h = 12)
plot(fcast_ets_1)

##Sin tendencia##
fit_2 <- ets(unem, model = "ZNZ")
fcast_ets_2 <- forecast(fit_2, h = 12)
plot(fcast_ets_2)