[R-br] Problemas com fixed breaks

[Mauro Sznelwar]

Tentei rodar de novo, mas deu dois problemas:

library(maptools)         
Carregando pacotes exigidos: foreign
Carregando pacotes exigidos: lattice
Checking rgeos availability: FALSE
        Note: when rgeos is not available, polygon geometry     computations in maptools depend on gpclib,
        which has a restricted licence. It is disabled by default;
        to enable gpclib, type gpclibPermit()

e depois


mapa = readShapePoly("saudemapa.shp")
Erro em getinfo.shape(filen) : Error opening SHP file

E ainda faltou disponibilizar o arquivo
logo_pb.png
   

  Em 1 de abril de 2012 21:05, Mauro Sznelwar <sznelwar em uol.com.br> escreveuo:

    Estava tentando rodar este código, e vi que falta este arquivoc:/users/Julianna/Documents/SecSaude/saudemapa.shp
    Daria para enviar?.

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