[R-br] Problemas com fixed breaks

[Benilton Carvalho]

Julianna,

e' possivel que sua variavel 'ano_2006' esteja num formato "diferente"
(isso e' so' uma hipotese maluca, nao tenho evidencias para isso)...
Entretanto, tente usar os dados a partir do arquivo que voce usou
originalmente para essa mensagem... dai' use o seguinte:

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

Veja se o problema se repete... Se sim, diga-nos o uqe vc obtem para:

class(dados$2006)
summary(dados$2006)

No meu caso, usando o seu csv do email anterior, eu obtive:

> class(dados$ano_2006)
[1] "numeric"
> summary(dados$ano_2006)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
 -1.000   0.000   0.000   9.956   0.000 100.000


benilton

On 30 March 2012 15:03, Benilton Carvalho <beniltoncarvalho em gmail.com>wrote:

> 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/fcaaa220/attachment.html>