Análisis básico de series de tiempo
para este primer ejemplo usaremos datos del consumo de combustible en españa de los años 1966 hasta 1977
scan("http://verso.mat.uam.es/~joser.berrendero/datos/jj.dat")
## [1] 0.710000 0.630000 0.850000 0.440000 0.610000 0.690000 0.920000
## [8] 0.550000 0.720000 0.770000 0.920000 0.600000 0.830000 0.800000
## [15] 1.000000 0.770000 0.920000 1.000000 1.240000 1.000000 1.160000
## [22] 1.300000 1.450000 1.250000 1.260000 1.380000 1.860000 1.560000
## [29] 1.530000 1.590000 1.830000 1.860000 1.530000 2.070000 2.340000
## [36] 2.250000 2.160000 2.430000 2.700000 2.250000 2.790000 3.420000
## [43] 3.690000 3.600000 3.600000 4.320000 4.320000 4.050000 4.860000
## [50] 5.040000 5.040000 4.410000 5.580000 5.850000 6.570000 5.310000
## [57] 6.030000 6.390000 6.930000 5.850000 6.930000 7.740000 7.830000
## [64] 6.120000 7.740000 8.910000 8.280000 6.840000 9.540000 10.260000
## [71] 9.540000 8.729999 11.880000 12.060000 12.150000 8.910000 14.040000
## [78] 12.960000 14.850000 9.990000 16.200000 14.670000 16.020000 11.610000
gas <- scan("http://verso.mat.uam.es/~joser.berrendero/datos/gas6677.dat")
Graficamos los datos
plot(gas)
Para tratar este vector númerico como una serie de tiempo, utilizaremos el comando ts (time-series objects)
gas.ts = ts(gas, start=c(1966,1),frequency = 12)
print(gas.ts)
## Jan Feb Mar Apr May Jun Jul Aug Sep
## 1966 92.718 91.380 111.643 118.888 119.432 127.796 158.943 178.013 143.385
## 1967 113.661 108.224 142.256 129.835 150.735 149.554 185.792 201.758 166.565
## 1968 135.951 126.615 146.647 165.822 163.365 169.294 215.538 233.427 184.402
## 1969 154.844 143.552 171.573 188.322 192.756 195.296 252.288 268.379 218.810
## 1970 179.759 173.821 211.387 210.551 218.371 232.057 294.173 312.700 251.891
## 1971 193.916 188.375 236.187 249.037 235.957 258.980 321.085 334.562 276.932
## 1972 225.010 225.742 265.159 271.986 290.953 285.108 362.687 386.347 314.205
## 1973 268.578 256.063 312.041 326.741 315.157 353.016 403.662 451.098 356.811
## 1974 289.186 296.881 302.589 334.091 325.790 337.782 423.297 454.172 353.727
## 1975 317.760 298.188 363.429 350.203 372.149 371.877 472.458 485.517 406.223
## 1976 352.200 334.938 372.891 397.388 385.657 416.961 492.480 512.209 411.514
## 1977 363.367 342.979 384.936 421.718 402.877 427.615 538.254 528.007
## Oct Nov Dec
## 1966 127.179 114.403 124.900
## 1967 148.048 131.581 141.315
## 1968 178.432 155.179 163.355
## 1969 203.545 172.148 198.381
## 1970 235.560 202.876 224.383
## 1971 258.269 233.532 251.755
## 1972 292.124 261.740 291.810
## 1973 352.566 305.580 410.614
## 1974 353.413 315.272 341.902
## 1975 377.262 329.794 384.350
## 1976 392.380 369.671 400.243
## 1977
Ahora que tenemos una variable que es un objeto orientado a tiempo, podemos tener una grafica en la cual se entienda la perodicidad de los aumentos del consumo de gas
plot(gas.ts, type="o", bg="blue", pch = 21, col="red", cex=0.7, xlab="anios",ylab="consumo")
Ahora haremos una comparacion interanual del aumento del consumo del combustible en españa
boxplot(gas.ts ~ cycle(gas.ts), bg="blue", pch=21, col="#1500a7c9", xlab="Meses", ylab="Combustible")
Entendamos los ciclos de comportamiento de consumo de gasolina per capita
cycle(gas.ts)
## Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
## 1966 1 2 3 4 5 6 7 8 9 10 11 12
## 1967 1 2 3 4 5 6 7 8 9 10 11 12
## 1968 1 2 3 4 5 6 7 8 9 10 11 12
## 1969 1 2 3 4 5 6 7 8 9 10 11 12
## 1970 1 2 3 4 5 6 7 8 9 10 11 12
## 1971 1 2 3 4 5 6 7 8 9 10 11 12
## 1972 1 2 3 4 5 6 7 8 9 10 11 12
## 1973 1 2 3 4 5 6 7 8 9 10 11 12
## 1974 1 2 3 4 5 6 7 8 9 10 11 12
## 1975 1 2 3 4 5 6 7 8 9 10 11 12
## 1976 1 2 3 4 5 6 7 8 9 10 11 12
## 1977 1 2 3 4 5 6 7 8