[R-br] Problemas com fixed breaks

[Benilton Carvalho]

Mensagem abaixo em nome da Julianna... (dados disponiveis via datafilehost)

====

E quando fui conferir os breaks, os intervalos estavam certos,mas o mapa
resultante estava errado.
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???

segue o código:

#######################################################
# CONSTRUINDO MAPAS
#######################################################
##Code (Comments are preceded by ##)
## Load required packages

# Limpando memoria.
rm(list=ls(all=TRUE))

library(maps)
library(sp)
library(maptools)
library(spdep)
library(classInt)
library(RColorBrewer)
library(shape)
library(SDMTools)
library(latticeExtra)
library(ReadImages)
library(png)

rosa_dos_ventos = function(loc=c(-34.5, -8), size = 0.15, bearing=0,
cols=1, cex=0.6,...)
{
  cols <- rep(c("white","black"),8)
  radii <- rep(size/c(1,4,2,4),4)
  x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1]
  y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2]
  for (i in 1:15)
  {
    x1 <- c(x[i],x[i+1],loc[1])
    y1 <- c(y[i],y[i+1],loc[2])
    polygon(x1,y1,col=cols[i])
  }
  polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16])
  b <- c("O","S","L","N")
  for (i in 0:3) text((size+par("cxy")[1]-0.40)*cos(bearing+i*pi/2)+loc[1],

(size+par("cxy")[2]-0.42)*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex)
}

## Lendo dados
## http://www.datafilehost.com/download-e5cb7e86.html
dados = read.table("dadosJulianna.csv", sep = ",", header = TRUE)

# Lendo mapa.
mapa = readShapePoly("c:/users/Julianna/Documents/SecSaude/saudemapa.shp")
proj4string(mapa) <- CRS("+init=epsg:4291")
scaleXY(mapa, 1836656)


# arquivo png
#img <- readPNG("c:/users/Julianna/Documents/SecSaude/logo_pb.png") # seu
arquivo png aqui
#set.seed(123)
#par(xpd = TRUE)


colors <- brewer.pal(7, "Reds")
brks <- classIntervals(dados$ano_2006, n = 7,style="fixed",
intervalClosure="right",
                       fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, 10000),
main="Fixed")
brks<- brks$brks

ID =
match(substr(mapa$CODIGO,1,6),dados$t_indicadores_municipio.f_municipio)
dados = dados[ID,]

png(filename="c:/mif2006.png", height=750, width=1500)


plot(mapa, col = colors[findInterval(dados$ano_2006,
brks,all.inside=TRUE)], axes = TRUE, cex = 1.5, mar=5)
##, xlim=c(-39.7, -34.7), ylim=c(-8.5, -6)

##invisible(text(getSpPPolygonsLabptSlots(mapa),
labels=as.character(mapa$NOMEMUN_1), cex=0.65))

title(main = "PORCENTAGEM DE Ã BITOS INVESTIGADOS\ntipo: mulheres em idade
fértil",
      font.main = 2, cex.main = 1.8)

text(-35.7, -6.2, "ANO: 2006", cex = 1.8)

legend(-39.6,-7.6, c("Sem informação", "0", "]0,20]", "]20,40]",
"]40,60]", "]60,80]","Acima de 80"),
       fill = colors, bty="n", x.intersp = 1.1, y.intersp = 1.1, cex = 1.6)

rosa_dos_ventos(loc=c(-34.5,-8),size=0.17,cex=1.20, bearing=0, cols=1)

SpatialPolygonsRescale(layout.scale.bar(), offset = c(-35,-8.30),
                       scale = 1, fill = c("transparent", "black"),
                       plot.grid= FALSE)

text(-34.5, -8.35, "110.3 KM", cex= 1.3)

#rasterImage(img, xleft = -39.5, ybottom = -6.25, xright = -38.9 , ytop =
-6)

dev.off()
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120330/0347b7f5/attachment.html>