IPC_men <- Serie_IPC_1_2010_a_2_2023
¿que bibliotecas son necesarias en R para analizar series temporales?
inflacion = ts(IPC_men, start = c(2010,1), frequency = 12)
plot(inflacion)
lnIPC <- log(IPC.ts)
Se crea la variable de diferencias
DlnIPC = diff(lnIPC)
grafico
plot(DlnIPC, main=”variacion porcentual del IPC”, col= “blue”)
modelo1 <- dynlm(lnIPC~L(lnIPC))
o Respuesta: modelo4 <- dynlm(lnIPC~L(lnIPC, c(1,2,5,8)))
tend=seq_along(lnIPC)
modelo5 <- dynlm(lnIPC~L(lnIPC, 1:2 ) + trend(lnIPC))
modelo6 <- dynlm(lnIPC~L(lnIPC, 1:2) + trend(lnIPC) + season(lnIPC))
como identifica si un modelo tiene auto correlacion(varianza) y correlacion parcial (covarianza)?
acf(Padj)
pacf(Padj)
como generar un grafico de tendencia lineal?
tsplot(Padj)
para que sirve el comando acf?
detecta autocovarianza y autocorrelacion. Los límites de error son límites aproximados de ruido blanco
pp.test(Padj)
prueba de Phillips Perron
adf.test(Padj)
Dickey-Fuller Test
DPadj=diff(Padj)
la serie temporal pasa a ser en diferencias, se puede estimar mejor con una variables logarítmica para incluir niveles poblacionales
pp.test(DPadj)
adf.test(DPadj)
A travez de metodo formal se verifica si la raiz es unitaria, o sea con media 0 y varianza unica
windows()
acf2(Dpadj)
Genera un gráfico de estimaciones donde muestra los valores de autocorrelación y autocorrelación parcial para diferentes rezagos.
Las barras significativas en la ACF indican correlaciones significativas en esos rezagos. Si están dentro del intervalo de confianza (linea azul) no hay evidencia significativa de autocorrelación.
ndiffs(DPadj)
propone un modelo significativo
model7 <- sarima(DPadj, 2, 0,2)
model7$ttable
La función sarima ajusta modelos ARIMA extendidos con componentes
1. estacionales
2. tendencias,
3. variables exógenas y
4 raíces arbitrarias en el círculo unitario.
Puede utilizarse para
realizar regresiones con errores autocorrelacionados proporciona estimaciones de:
1. parámetros,
2. errores estándar,
3. criterios de información (AIC, AICc, BIC)
4. diagnósticos.
model8 <- arima(DPadj, order = c(1,0,2))
windows()
tsdiag(model8)
el comando ARIMA se utiliza para ajustar modelos de series temporales, y además, puede aplicarse en regresiones donde los errores están autocorrelacionados. Los diagnósticos proporcionan información sobre la calidad del ajuste del modelo.
Box.test(residuals(model8), type = “Ljung-Box”)
Test de ruido blanco Ho: La prueba nos indica que los errores son de ruido blanco