Mest3 U1 A1 Elvc
Mest3 U1 A1 Elvc
Mest3 U1 A1 Elvc
Grupo: MT-MEST3-2201-B1-000
.
Actividad 1
gas = scan('http://verso.mat.uam.es/~joser.berrendero/datos/gas6677.dat')
plot(gas)
Vemos que el gráfico resultante no es el más apropiado para describir una serie temporal. Si
queremos que R trate a un objeto como serie temporal, tenemos que determinar
apropiadamente sus características con el comando ts. Para definir la serie correctamente
escribimos:
gas.ts = ts(gas, start = c(1966,1), frequency = 12)
print(gas.ts)
plot(gas.ts)
Vemos ahora cómo el resultado depende de las características que hemos definido para la
serie. Si queremos comparar la distribución del consumo de gasolina para cada mes, un
gráfico útil es
boxplot(gas.ts ~ cycle(gas.ts))
El comando cycle determina la unidad de tiempo a la que pertenece cada observación de la
serie:
cycle(gas.ts)
Por lo tanto, en el gráfico anterior se ha producido un diagrama de cajas para cada mes del
año.
Según el diagrama en 1978. Porque en cuanto más nos acercamos a la época actual hay más
personas que ya tienen vehículo y por lo tanto tienen que consumir la gasolina
boxplot(gas.ts ~ cycle(gas.ts))
Con R es muy sencillo obtener una descomposición estructural de este tipo. Se usa el
comando decompose:
gas.ts.desc = decompose(gas.ts)
plot(gas.ts.desc, xlab='Año')
Referencias: