[R-es] Mapa de quantiles con spplot

[Rolando Valdez]

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