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 <comunello.eder@gmail.com> 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 <comunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]