Mensagem abaixo em nome da Julianna... (dados disponiveis via datafilehost) <div><br></div><div>====</div><div><br></div><div><div>E quando fui conferir os breaks, os intervalos estavam certos,mas o mapa resultante estava errado.</div>
<div>Porque quando eu não importo o csv e uso o dataframe resultante da função, os breaks ficam ok(da maneira que eu especifico em FixedBreaks) e, do contrario, não???</div><div><br></div><div>segue o código:</div><div><br>
</div><div>#######################################################</div><div># CONSTRUINDO MAPAS</div><div>#######################################################</div><div>##Code (Comments are preceded by ##)</div><div>
## Load required packages</div>
<div><br></div><div># Limpando memoria.</div><div>rm(list=ls(all=TRUE))</div><div><br></div><div>library(maps)</div><div>library(sp)</div><div>library(maptools) </div><div>library(spdep) </div><div>library(classInt) </div>
<div>library(RColorBrewer) </div><div>library(shape)</div><div>library(SDMTools)</div><div>library(latticeExtra)</div><div>library(ReadImages)</div><div>library(png)</div><div><br></div><div>rosa_dos_ventos = function(loc=c(-34.5, -8), size = 0.15, bearing=0, cols=1, cex=0.6,...)</div>
<div>{</div><div> cols <- rep(c("white","black"),8)</div><div> radii <- rep(size/c(1,4,2,4),4)</div><div> x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1]</div><div> y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2]</div>
<div> for (i in 1:15)</div><div> {</div><div> x1 <- c(x[i],x[i+1],loc[1])</div><div> y1 <- c(y[i],y[i+1],loc[2])</div><div> polygon(x1,y1,col=cols[i])</div><div> }</div><div> polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16])</div>
<div> b <- c("O","S","L","N")</div><div> for (i in 0:3) text((size+par("cxy")[1]-0.40)*cos(bearing+i*pi/2)+loc[1],</div><div> (size+par("cxy")[2]-0.42)*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex)</div>
<div>}</div><div><br></div><div>## Lendo dados</div><div>## <a href="http://www.datafilehost.com/download-e5cb7e86.html">http://www.datafilehost.com/download-e5cb7e86.html</a></div><div>dados = read.table("dadosJulianna.csv", sep = ",", header = TRUE)</div>
<div><br></div><div># Lendo mapa.</div><div>mapa = readShapePoly("c:/users/Julianna/Documents/SecSaude/saudemapa.shp")</div><div>proj4string(mapa) <- CRS("+init=epsg:4291")</div><div>scaleXY(mapa, 1836656)</div>
<div><br></div><div><br></div><div># arquivo png</div><div>#img <- readPNG("c:/users/Julianna/Documents/SecSaude/logo_pb.png") # seu arquivo png aqui</div><div>#set.seed(123)</div><div>#par(xpd = TRUE)</div>
<div>
<br></div><div><br></div><div>colors <- brewer.pal(7, "Reds")</div><div>brks <- classIntervals(dados$ano_2006, n = 7,style="fixed", intervalClosure="right",</div><div> fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, 10000), main="Fixed")</div>
<div>brks<- brks$brks</div><div><br></div><div>ID = match(substr(mapa$CODIGO,1,6),dados$t_indicadores_municipio.f_municipio)</div><div>dados = dados[ID,]</div><div><br></div><div>png(filename="c:/mif2006.png", height=750, width=1500)</div>
<div><br></div><div><br></div><div>plot(mapa, col = colors[findInterval(dados$ano_2006, brks,all.inside=TRUE)], axes = TRUE, cex = 1.5, mar=5)</div><div>##, xlim=c(-39.7, -34.7), ylim=c(-8.5, -6)</div><div><br></div><div>
##invisible(text(getSpPPolygonsLabptSlots(mapa), labels=as.character(mapa$NOMEMUN_1), cex=0.65))</div><div><br></div><div>title(main = "PORCENTAGEM DE à BITOS INVESTIGADOS\ntipo: mulheres em idade fértil",</div>
<div> font.main = 2, cex.main = 1.8)</div><div><br></div><div>text(-35.7, -6.2, "ANO: 2006", cex = 1.8)</div><div><br></div><div>legend(-39.6,-7.6, c("Sem informação", "0", "]0,20]", "]20,40]", "]40,60]", "]60,80]","Acima de 80"),</div>
<div> fill = colors, bty="n", x.intersp = 1.1, y.intersp = 1.1, cex = 1.6)</div><div><br></div><div>rosa_dos_ventos(loc=c(-34.5,-8),size=0.17,cex=1.20, bearing=0, cols=1)</div><div><br></div><div>SpatialPolygonsRescale(layout.scale.bar(), offset = c(-35,-8.30),</div>
<div> scale = 1, fill = c("transparent", "black"),</div><div> plot.grid= FALSE)</div><div><br></div><div>text(-34.5, -8.35, "110.3 KM", cex= 1.3)</div>
<div><br></div><div>#rasterImage(img, xleft = -39.5, ybottom = -6.25, xright = -38.9 , ytop = -6)</div><div><br></div><div>dev.off()</div><div><br></div></div>