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. Tasas 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

 

Deja un comentario

Tu dirección de correo electrónico no será publicada. Los campos necesarios están marcados *


*

Puedes usar las siguientes etiquetas y atributos HTML: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong>

Protected by WP Anti Spam