[R-es] Mapa de quantiles con spplot

[Rolando Valdez]

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 en 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 <onunez en unex.es> 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
>> R-help-es en r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>> 
> 
> 
> --
> ____________________________________
> 
> 
> Olivier G. Nuñez
> Email: onunez en unex.es
> http://kolmogorov.unex.es/~onunez
> Tel : +34 663 03 69 09
> Departamento de Matemáticas
> Universidad de Extremadura
> 
> 
> ____________________________________
> 
> 
> 

Rolando Valdez