Descomposición estacional en R (Parte I)

He estado escribiendo estos tutoriales para el blog institucional de la UCA. Agradezco en nombre de la iniciativa de Usuarios R Nicaragua el espacio otorgado. El blog está en: http://blog.uca.edu.ni/rnica/ . Puede ser consultado, comentado y quienes deseen agregar aportes, con gusto. El blog está enfocado a la programación para el tratamiento gráfico y estadísticos de datos.

Veamos la descomposición estacional de series de tiempo. En este ejemplo vamos a usar una base de datos puesta a disposición por Paul S.P. Cowpertwait. Tasa de desempleo de los Estados Unidos desde enero 1996 hasta octubre del 2006.

Especificando la dirección web de los datos
USdesempleo <- 'http://staff.elena.aut.ac.nz/Paul-Cowpertwait/ts/USunemp.dat'

Importando los datos
US.mensual <- read.table(USdesempleo, header = T)

Visualizando los datos
View(US.mensual)

Se nota no es una serie de tiempo. Entonces hacemos la conversión.
attach(US.mensual)
US.mensual.ts <- ts(USun, start=c(1996,1), end=c(2006,10), freq=12)
US.mensual.ts


Note hemos importado una serie de tiempo.


     Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1996 5.6 5.5 5.5 5.6 5.6 5.3 5.5 5.1 5.2 5.2 5.4 5.4
1997 5.3 5.2 5.2 5.1 4.9 5.0 4.9 4.8 4.9 4.7 4.6 4.7
1998 4.6 4.6 4.7 4.3 4.4 4.5 4.5 4.5 4.6 4.5 4.4 4.4
1999 4.3 4.4 4.2 4.3 4.2 4.3 4.3 4.2 4.2 4.1 4.1 4.0
2000 4.0 4.1 4.0 3.8 4.0 4.0 4.0 4.1 3.9 3.9 3.9 3.9
2001 4.2 4.2 4.3 4.4 4.3 4.5 4.6 4.9 5.0 5.3 5.5 5.7
2002 5.7 5.7 5.7 5.9 5.8 5.8 5.8 5.7 5.7 5.7 5.9 6.0
2003 5.8 5.9 5.9 6.0 6.1 6.3 6.2 6.1 6.1 6.0 5.9 5.7
2004 5.7 5.6 5.7 5.5 5.6 5.6 5.5 5.4 5.4 5.4 5.4 5.4
2005 5.2 5.4 5.1 5.1 5.1 5.0 5.0 4.9 5.1 4.9 5.0 4.9
2006 4.7 4.8 4.7 4.7 4.6 4.6 4.8 4.7 4.6 4.4 
 
 


Si ha quedado como el ejemplo, ha sido un éxitos. 

Obtenemos opcional una visualización gráfica de la serie:

plot(US.mensual.ts, col = "Blue", ylab="Tasas de desempleo", xlab="Meses", lwd = 3)
 
Para una descomposición multiplicativa aplicamos las líneas:

US.mensual.ts.decom <- decompose(US.mensual.ts, type = "multiplicative")
plot(US.mensual.ts.decom)

 


Para una descomposición aditiva aplicamos las líneas:



 US.mensual.ts.decom <- decompose(US.mensual.ts, type = "additive")

plot(US.mensual.ts.decom)




En la parte 2 se enseñará a extraer los datos de las gráficas. 

Deybi Morales
Economista
morales.economia@gmail.com




Comentarios

Envíanos un mensaje

Nombre

Correo electrónico *

Mensaje *