Esta pregunta ha sido convertida a Wiki comunitaria y wiki bloqueada porque es un ejemplo de pregunta que busca una lista de respuestas y parece ser lo suficientemente popular como para protegerla del cierre. Es debe ser tratado como un caso especial y no debe ser visto como el tipo de pregunta que se fomenta en este, o cualquier sitio de Stack Exchange pero si deseas contribuir con más contenido, siéntete libre de hacerlo editando esta respuesta. de hacerlo editando esta respuesta.
Esto no es tanto un truco como spplot()
de la empresa. spplot()
La capacidad de la leyenda de escalar las muestras (para que coincidan con los rangos de clasificación) sirve como una herramienta pedagógica útil cuando se discute la distribución de los datos de atributos y los tipos de clasificación. La combinación de gráficos de distribución acumulativa con los mapas ayuda en esta tarea.
Los estudiantes sólo tienen que modificar unos pocos parámetros de script para explorar los tipos de clasificación y los efectos de la transformación de datos. Esta suele ser su primera incursión en R en lo que es un curso centrado principalmente en ArcGIS.
Aquí tienes un fragmento de código:
library(rgdal) # Loads SP package by default
NE = readOGR(".", "NewEngland") # Creates a SpatialPolygonsDataFrame class (sp)
library(classInt)
library(RColorBrewer)
pal = brewer.pal(7,"Greens")
brks.qt = classIntervals(NE$Frac_Bach, n = 7, style = "quantile")
brks.jk = classIntervals(NE$Frac_Bach, n = 7, style = "jenks")
brks.eq = classIntervals(NE$Frac_Bach, n = 7, style = "equal")
# Example of one of the map plots
spplot(NE, "Frac_Bach",at=brks.eq$brks,col.regions=pal, col="transparent",
main = list(label="Equal breaks"))
# Example of one of the cumulative dist plots
plot(brks.eq,pal=pal,main="Equal Breaks")
Ref: Applied Spatial Data Analysis with R (R. Bivand, E Pebesma & V. Gomez-Rubio)