29 dic 2011

Grafico de variables pareadas


Otra manera de hacer graficos elgantes  con una  extención de libreria ggplot2 que es GGally de R (extraído de SAS and R).

library(GGally)
ds = read.csv("http://www.math.smith.edu/r/data/help.csv")
ds$sex = as.factor(ifelse(ds$female==1, "female", "male"))
ds$housing = as.factor(ifelse(ds$homeless==1, "homeless", "housed"))
smallds = subset(ds, select=c("housing", "sex", "i1", "cesd"))
ggpairs(smallds, diag=list(continuous="density", discrete="bar"), axisLabels="show")

19 dic 2011

Mapa de Antioquia en R

# El mapa con los municiopios se puede bajar de aqui http://sigotn.igac.gov.co/sigotn/default.aspx ##
library(maptools)
library(gplots)
library(RColorBrewer)
Municipios<- readShapePoly("Municipios.shp")
rgb.palette <- colorRampPalette(c("red", "orange","blue"), space = "Lab")
Antioquia<-Municipios[Municipios$NOM_DEPART=="ANTIOQUIA",]
plot(Antioquia,border="black",col=colorRampPalette(brewer.pal(9,"Blues"))(90),lwd=0.9,asp=0,axes =T, las = 2)

Mapa de Colombia en R

library(maptools)
library(ggplot2)
library(lattice)
library(sp)
#El archivo .shp se puede bajar de http://www.gadm.org/country #
xx <- readShapePoly("COL_adm1.shp",IDvar="NAME_1", proj4string=CRS("+proj=longlat +ellps=clrk66"))
mapa<-function(xx,depto,sati)
{
  col<-matrix(NA,ncol=2,nrow=length(depto))
  color<-numeric(length(depto))
  for(i in 1:length(depto)){
  ifelse(sati[i]<73.3,color[i]<-"red", ifelse(sati[i]<85.5, color[i]<-"yellow", color[i]<-"darkgreen"))}
  col[,1]<-depto;col[,2]<-color
  col<-col[order(col[,1]),]
  posi<-numeric(length(depto))
  xxx <-xx[xx$NAME_1 %in%depto,]
  plot(xx,border="gray",col="white",lwd=0.9,asp=0,axes =F, las = 1)
  plot(xxx, border="gray", col=col[,2],lwd=1,asp=0,add=T)
  text(coordinates(xxx), labels=sapply(slot(xxx, "polygons"),function(i) slot(i, "ID")), cex=0.65,font=2, pos=1, col="black")
}
mapa(xx,c("Norte de Santander","Quindío","Bolívar","Santander","Atlántico","Valle del Cauca","Cundinamarca","Antioquia"),c(83,78,76.9,75.6,74.3,73.8,70.1,68.6))