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 <rvaldezr en gmail.com> 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