No me sale ningún error, simplemente el resultado no es el esperado. 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. El 05/05/2014, a las 05:44, Oscar Perpiñan <[email protected]> escribió: > Hola Rolando, > > Por el código que muestras debiera funcionar, pero falta información > para poder responderte. > Cuando dices que no funcionó, ¿que significa? ¿Obtienes un error? ¿El > resultado no es el esperado? > Por favor, envía los datos que estás utilizando, al menos una muestra, > o algo que sirva para poder reproducir tu código. > > Saludos. > > Oscar. > ----------------------------------------------------------------- > Oscar Perpiñán Lamigueiro > Dpto. Ingeniería Eléctrica (ETSIDI-UPM) > Grupo de Sistemas Fotovoltaicos (IES-UPM) > URL: http://oscarperpinan.github.io > Twitter: @oscarperpinan > > > El día 3 de mayo de 2014, 21:39, Rolando Valdez <[email protected]> escribió: >> 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 Rolando Valdez _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
