Hola, Lo hice de la forma en que propones y no me sale, no sé si algo esté mal con R, porque no me respeta ni el ‘key.space’, me lo sigue poniendo en la parte izquierda del gráfico
Saludos El 06/05/2014, a las 02:21, Olivier Nuñez <[email protected]> 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 > >> 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@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 <[email protected]> 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 >>>> [email protected] >>>> https://stat.ethz.ch/mailman/listinfo/r-help-es >>>> >>> >>> >>> -- >>> ____________________________________ >>> >>> >>> Olivier G. Nuñez >>> Email: [email protected] >>> http://kolmogorov.unex.es/~onunez >>> Tel : +34 663 03 69 09 >>> Departamento de Matemáticas >>> Universidad de Extremadura >>> >>> >>> ____________________________________ >>> >>> >>> >> >> Rolando Valdez >> >> >> >> > > > -- > ____________________________________ > > > Olivier G. Nuñez > Email: [email protected] > http://kolmogorov.unex.es/~onunez > Tel : +34 663 03 69 09 > Departamento de Matemáticas > Universidad de Extremadura > > > ____________________________________ > > > Rolando Valdez _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
