En primer lugar, y retomando lo publicado en la anterior entrada, mapeamos exclusivamente en la localidad Heróica Puebla de Zaragoza y con las divisiones por AGEB (Área Geoestadística Básica), con los siguientes comandos:
shp2dta using ageb_urb, database(agebdb) coordinates(agebcoord)
use agebdb, clear
spmap using agebcoord if cve_loc==211140001, id(_ID) fcolor(blue*0.5)
El siguiente paso será mapear el porcentaje de Población Económicamente Activa (PEA) por AGEB, de la siguiente manera:
spmap PEA using agebcoord if cve_loc==211140001, id(_ID) clnumber(20) fcolor(Blues2) ocolor(none ..) legstyle(3) legend(ring(1) position(3)) ndfcolor(red)
En donde:
clnumber(20) = Número de rangos que deseamos en la variable.
fcolor(Blues2) = Esquema de coloreo que deseamos en los polígonos, este tipo en particular admite hasta 99 categorias, y funciona perfectamente en variables continuas.
ocolor(none ..) = Indica que no queremos colorear el contorno de los polígonos.
legstyle(3) = Índica el estilo de la leyenda, este tipo es el indicado para este ejemplo.
ndfcolor(red) = El color para los polígonos donde no se tiene información.
Me encanto..¡¡ Sólo incorporaría el comando para establecer el directorio de trabajo en donde se encuentren las capas de información.
ResponderBorrarGracias, si sería una opción, lo que usualmente hago es colocar el archivo .do (el script de Stata) en donde tengo mis capas de información (.dbf; .prj; .shp; .shx). Pero efectivamente, sería más elegante invocar el directorio de las capas de información con el comando cd. De la siguiente manera:
ResponderBorrarcd C:\Usuario\AlexPlatas\MisCapas
shp2dta using ageb_urb, database(agebdb) coordinates(agebcoord)
use agebdb, clear
Saludos