Valeu Éder. Vou ver oq e você sugeriu. Muito obrigado!!!
Tenho uma outra pergunta, o terraview, usa escala de cores par mostrar essas
quantidades. Como controlar isso no R?.
[ ]'s.
Edson Lira
Estatístico
Manaus-Amazonas
Em Sexta-feira, 23 de Maio de 2014 13:17, Éder Comunello
<[email protected]> escreveu:
Edson, boa tarde!
Acredito que você deveria usar um "mapa de zonas" no lugar do "mapa de
BAIRROS".
No link que segue estão os arquivos (shapefile) para uso no script (ao final):
http://dl.dropbox.com/s/5dpo1ad1oez2752/manausZon.zip
### <code r>
# rm(list=ls())
library(maptools)
setwd(choose.dir());getwd() ### defina o diretório de trabalho
mapZonas <- readShapePoly("manausZon.shp")
mapZonas@data
{plot(mapZonas, col=heat.colors(8), axes=T)
with(mapZonas@data, text(X1, X2, ZONAS, cex=.5))
title("Manaus - Zonas")}
apAll <- read.csv2("doa_apto_00_12.csv"); names(apAll)
ap2011 <- subset(apAll, TRUE, select = c("ZONAS", "ano2011")); ap2011
dados <- merge(mapZonas@data, ap2011)
mapZonas@data <- dados
quant <- quantile(dados$ano2011, include.lowest=TRUE); quant
breaks <- c(floor(quant[-length(quant)]/10),
ceiling(quant[length(quant)]/10))*10; breaks
dcut <- cut(dados$ano2011, breaks, right=F, dig.lab=4); dcut
levels(dcut)
cores <- c("#779999", "#99BBBB", "#BBDDDD", "#DDFFFF") ### cores <- palette()
# cores <- rev(heat.colors(4))
# x11()
plot(mapZonas, axes=T, col=cores[dcut])
legend(legend=levels(dcut), "topright", bty="n", fill=cores, cex=.8)
title("Doações aptas em Manaus por zona \n 2011")
# with(mapZonas@data, text(X1, X2, ZONAS, cex=.5)) text(X1, X2, ZONAS, cex=.5))
tmp <- with(mapZonas@data, paste(ZONAS, '\n', ano2011))
with(mapZonas@data, text(X1, X2, tmp, cex=.5))
### </code>
Éder Comunello <[email protected]>
Dourados, MS - [22 16.5'S, 54 49'W]_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
