Hola, no me funciona: Te contesto lo mismo que le contesté a otro compañero de la lista, espero ser más claro.
Les comparto la carpeta con la capa y los datos que estoy usando. https://www.dropbox.com/s/umy0evt3qm1wr4d/dissolve.zip Esto es lo que estoy haciendo: library(rgdal) library(maptools) library(sp) library(classInt) library(RColorBrewer) zm <- readOGR(“.”, “zmdis”) data <- read.csv(“part_pot.csv”, header=T, sep=“,”, dec=“.”) zm@data=data.frame(data) plotvar88 <- zm$part88 nclr <- 8 plotclr <- brewer.pal(nclr, “Blues”) class <- classIntervals(plotvar88, nclr, style=“quantile”) colcode <- findColours(class, plotclr) plot(zm) plot(zm, col=colcode, add=T) legend(3777189,2249615, legend=names(attar(colcode, “table”)), fill=attr(colcode, “palette”), cex=0.6, bty=“n”) Entonces obtengo el mapa como lo quiero (Me falta la división política, pero para este ejemplo la omití): https://www.dropbox.com/s/w3ps61o4jr1ciy5/part88.tiff Cuando lo hago con spplot: spplot(zm, "part88", col.regions=plotclr, cuts=nclr, key.space="bottom”) Me regresa lo siguiente: https://www.dropbox.com/s/2cgiua5ra2tx4a0/part88_spplot.tiff Entonces lo que quiero es usar el spplot para pintar dos variables c(“part88”, “part93”), que el mismo ‘legend’ aplique para ambas. Saludos cordiales El 05/05/2014, a las 10:39, Olivier Nuñez <[email protected]> escribió: > No resulta fácil contestarte con la información que proporcionas. > Sobre la base de las variables que definiste, intuyó que lo siguiente debería > funcionar: > > zm$clases <- class > spplot(zm, "clases" , col.regions=plotclr) > > Un saludo. Olivier > >> Hola, >> >> Intento representar en un mapa participaciones porcentuales de los sectores >> económicos y no logro hacerlo con spplot. >> >> He intentado con spplot(zm, c(“part88”, “part93”), cuts=4, >> col.regions=brewer.pal(4, >> “Set3”)) >> >> Lo pude hacer utilizando el base graphics de R, definiendo: >> >>> plotvar88 <- zm$part88 >>> nclr <- 8 >>> plotclr <- brewer.pal(nclr, "PuOr") >>> plotclr <- plotclr[nclr:1] >>> class <- classIntervals(plotvar88, nclr, style="quantile") >>> colcode <- findColours(class, plotclr) >>> plot(zm) >>> plot(zm, col=colcode, add=T) >>> legend(3377189,2249615, legend=names(attr(colcode, "table")), >>> fill=attr(colcode, >>> "palette"), cex=0.6, bty="n”) >> >> Después intenté meter algunos de esos argumentos utilizando spplot >> >> spplot(zm, c(“part88”, “part93”), col.regions=plotclr, cuts=nclr) >> >> y tampoco funcionó. >> >> Tengo algunas nociones de cómo hacerlo en spplot, pero no sé cómo, pues >> desconozco >> la forma en utilizar los argumentos por ejemplo usando “cuts”, “at”. He >> revisado la >> ayuda que provee el paquete y sigo sin poder concretarlo. >> >> Agradecería un poco de orientación. >> >> Saludos cordiales. >> >> Rolando Valdez >> >> _______________________________________________ >> R-help-es mailing list >> [email protected] >> https://stat.ethz.ch/mailman/listinfo/r-help-es >> > > > -- > ____________________________________ > > > Olivier G. Nuñez > Email: [email protected] > http://kolmogorov.unex.es/~onunez > Tel : +34 663 03 69 09 > Departamento de Matemáticas > Universidad de Extremadura > > > ____________________________________ > > > Rolando Valdez _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
