Displaying 4 results from an estimated 4 matches for "part88".
2014 May 05
3
Mapa de quantiles con spplot
...pbox.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, b...
2014 May 06
3
Mapa de quantiles con spplot
Hola,
El problema con la propuesta de Olivier es que los intervalos son
diferentes para cada variable.
La forma sencilla y rápida es:
spplot(zm["part88"], col.regions=plotclr, at=class$brks)
Pero para que quede más elegante hay que dar algunos pasos más:
## Intervalos en forma character
op <- options(digits=4)
tab <- print(class)
options(op)
intChar <- names(tab)
## Indice de correspondencia entre elemento e intervalo
idxP88 <-...
2014 May 03
2
Mapa de quantiles con spplot
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,...
2014 May 07
2
Mapa de quantiles con spplot
...fico
Saludos
El 06/05/2014, a las 02:21, Olivier Nuñez <onunez en unex.es> escribió:
> El segundo argumento en spplot ha de definir el número/nombre del intervalo al que
> pertenece el area a pintar.
> Este código debería funcionar:
>
> zm$class88=findCols(classIntervals(zm$part88, nclr, style="quantile"))
> zm$class93=findCols(classIntervals(zm$part93, nclr, style="quantile"))
>
> spplot(zm, c("class88","class93"), col.regions=plotclr, cuts=(nclr-1),
> key.space="bottom")
>
> Un saludo. Olivier
>
>...