Displaying 5 results from an estimated 5 matches for "plotclr".
2014 May 05
3
Mapa de quantiles con spplot
...m1wr4d/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 obt...
2014 May 03
2
Mapa de quantiles con spplot
...a 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(colc...
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 <- findCols(class)
## Es sólo u...
2014 May 07
2
Mapa de quantiles con spplot
...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
>
>> 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 dat...
2010 May 14
2
help color coding map in R
...tate", plot = FALSE, fill = TRUE)
# Calculate the range of the map (with extra margins).
xl <- extendrange(us.map$range[1:2])
yl <- extendrange(us.map$range[3:4])
library(maps)
library(lattice)
library(latticeExtra)
library(RColorBrewer) # creates nice color schemes
library(classInt)
plotclr <- brewer.pal(nclr,"PuRd")
class <- classIntervals(NatSTSummaryHigh.abi$STMean, nclr, style="fisher" )
colcode <- findColours(class, plotclr)
# Plot a multi-panel map of all the states, and colour
xyplot(y~x | NatSTSummaryHigh.abi$PrimaryDX, data = state.center,groups...