# Instalamos los paquetes necesarios
install.packages("tidyverse")
install.packages("lubridate")
install.packages("forecast")
install.packages("tsibble")
install.packages("ggplot2")
Estadística Avanzada
FORMULARIO - Segunda Parcial
1 Preparación del entorno
1.1 Instalación paquetes necesarios
1.2 Activar los paquetes
# Activamos los paquetes
library(tidyverse)
library(lubridate)
library(forecast)
library(tsibble)
library(ggplot2)
2 Análisis de Series de tiempo
2.1 Visualización de los datos
%>%
base_datos ggplot(aes(x = fechas, y = Variable)) +
geom_line(color = "darkblue") +
labs(title = "Título", x = "Fecha", y = "Variable")
2.2 Conversión a serie de tiempo
<- ts(base_datos$variable, start = c(año, 1), frequency = 12) temp_ts
2.3 Descomposición
<- decompose(temp_ts)
descomp_temp plot(descomp_temp)
2.4 Ajuste del modelo ETS
<- ets(temp_ts)
modelo_temp_ets summary(modelo_temp_ets)
2.5 Pronóstico con ETS
<- forecast(modelo_temp_ets, h = 12)
pronostico_ets plot(pronostico_ets)
2.6 Ajuste del modelo ARIMA
<- auto.arima(temp_ts)
modelo_temp_arima summary(modelo_temp_arima)
2.7 Pronóstico con ARIMA
<- forecast(modelo_temp_arima, h = 12)
pronostico_arima plot(pronostico_arima)
2.8 Comparación de modelos ETS vs ARIMA
A continuación se presenta una comparación de los principales criterios de información entre los modelos:
Modelo | AIC | AICc | BIC |
---|---|---|---|
ETS(A,N,A) | |||
ARIMA (auto) |
Conclusión:
3 Números índices
3.1 Índice de precio
<- # <- Colocar el precio
precio_anterior <- # <- Colocar el precio
precio_actual <- (precio_actual / precio_anterior) * 100
indice_precio indice_precio
3.2 Precio relativo
<- precio_actual / precio_anterior
precio_relativo precio_relativo
3.3 Cadenas relativas
<- c( , , , ) # <- Colocar los valores
cadena <- cumprod(cadena / 100) * 100
indice_encadenado indice_encadenado
3.4 Cambio de periodo de base
<- c( , , ) # <- Colocar los valores
indice_viejo <- # <- Colocar el valor base
nuevo_base <- (indice_viejo / nuevo_base) * 100
indice_nuevo_base indice_nuevo_base
3.5 Fusión de dos series de número índice
<- c( , , ) # <- Colocar los valores
serie1 <- c( , , ) # <- Colocar los valores
serie2 <- serie1[ ] / serie2[ ] # <- Indicar las referencias en los []
ajuste <- serie2 * ajuste
serie2_ajustada serie2_ajustada
3.6 Índices de precios al consumidor (IPC)
<- c( , , ) # <- Colocar los valores
precios_anterior <- c( , , ) # <- Colocar los valores
precios_actual <- c( , , ) # <- Colocar los pesos
pesos
<- sum((precios_actual / precios_anterior) * pesos) * 100
ipc ipc
3.7 Poder de compra y deflación de los valores
<- # <- Colocar el valor
valor_nominal <- # <- Colocar el valor
indice_precio <- valor_nominal / (indice_precio / 100)
valor_real valor_real
3.8 Análisis de tendencia
<- : # <- Colocar los años
anios <- c( , , , , ) # <- Colocar los precios
precios <- data.frame(anios, precios)
df
<- lm(precios ~ anios, data = df)
modelo summary(modelo)
<- ggplot(df, aes(x = anios, y = precios)) +
plot geom_point() +
geom_smooth(method = "lm", se = FALSE) +
labs(title = "Tendencia de Precios", x = "Año", y = "Índice")
plot