Archivo por meses: marzo 2015

Hacer un mapa de Nicaragua en R

El gran número de programadores que han apoyado a la iniciativa libre de R, ha influído en muchos de sus usos avanzados, algunas de estos usos apenas con mínima relación con las estadísticas. Ahora en R podemos crear páginas web de informes, hacer presentaciones, elaborar documentos, gráficos dinámicos, analizar canciones, mapas, etc. Lo bueno que nos disponen de sus paquetes para trabajar estos elementos que hace años se pensaban en R y sin conocimientos de programación no llegaríamos tan lejos. Ahora pues veremos como trabajar con mapa.

Los trazos y datos de un mapa están contenidos en un archivo llamado shapefile, existen otras extensiones a shp  pero esta es más conocida en el mundo del trazo de mapas.

Primero necesitamos ese shapefile y segundo, tener claro qué datos vamos a representar en el mapa.

El conjunto de paquetes a usar son:

install.packages(“sp”)

install.packages(“maptools”)

install.packages(“ggplot2″)

install.packages(“maptools”)

Continuar…

Rplot1

Deybi Morales León

Economista

morales.economia@gmail.com

 

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

 

Extraer lista de Sismos de Ineter con R

En este post quiero mostrarle un truco para extraer la lista de los sismos del día desde la página del Ineter utilizando el programa R.

1.Instalar el paquete rvest y activarlo
2.Especificar la página web con el html de los datos que publica www.ineter.com, a este truco le llamaría barrer todo html en busca del origen de la información, es la página web que necesitamos.
3.Visualizar el html de origen de la información para ubicar el prefijo que separa la información que queremos.
4.Por último si quiero imprimir en la consola solo ejecuto datos, a pesar de que mi última línea ha creado un archivo de txt llamado Sismos en mis documentos para que ubiques y lo habras.

##########################################
install.packages(“rvest”)
library(rvest)
Sismos <- html(“http://webserver2.ineter.gob.ni/geofisica/sis/events/sismos.php”)
datos <- Sismos %>%
html_nodes(“pre”) %>%
html_text()
write.table(datos, ‘Sismos.txt’, sep=’\t’)
############################################

Para visualizar los html ejecute:
html(“http://www.ineter.gob.ni“)
html(“http://webserver2.ineter.gob.ni/geofisica/sis/events/sismos.php”)