Los mapas como éste no los hago habitualmente con SAS, pero puede interesaros saber cómo se pueden hacer. Para la realización de mapas a partir de shapefiles suelo usar R; los que seáis habituales del blog sabéis que hay entradas al respecto. Sin embargo, recientemente he conocido QGIS, que interactúa a la perfección con Excel.
El ejemplo de hoy es para aquellos que no conocéis QGIS o R y estáis más habituados a trabajar con SAS. Se trata de importar un shapefile con SAS, poder hacer un sencillo tratamiento de datos y posteriormente realizar un mapa muy simple. Para ilustrar el ejemplo, vamos a realizar un mapa por municipios de una provincia española con SAS y, para ello, necesitamos un fichero .shp que os vais a descargar de http://www.gadm.org/country; seleccionáis los mapas de España y descargáis el ZIP que contiene los archivos.
El primer paso es crear un conjunto de datos SAS a partir del shapefile:
proc mapimport datafile="/directorio/mapa_sas/ESP_adm4.shp"
out=work.espania;
run;
Se emplea el PROC MAPIMPORT, que genera un conjunto de datos SAS compuesto por las coordenadas, identificadores y nombres de los municipios. Hacemos un tratamiento de datos muy sencillo:
data cadiz;
set espania;
if name_2 = "Cádiz";
run;
proc sort data=cadiz out=municipios_lista nodupkey;
by id_4;
run;
data municipios_datos;
set municipios_lista;
aleatorio = ranpoi(3, 4);
run;
Seleccionamos la provincia de Cádiz y, sobre sus municipios, calculamos una variable aleatoria que será la que representemos en SAS con el PROC GMAP:
proc gmap data=municipios_datos map=cadiz;
id id_4;
choro aleatorio / discrete;
run;
quit;
Ya sabéis hacerlo con SAS, pero os recomiendo que vayáis instalando QGIS. Saludos.