SASLibraries/Maps下有一個檔案Taiwan, 記錄了台灣各縣市的經緯度,透過SASproc gmap,我們可以將表格內容地圖化,以民國100年各縣市之標準化死亡率為例,透過地圖化,可以很快地看出死亡率的高低與各縣市的都市化程度有關。[參見台灣地圖]


[程式一] 建立民國100年各縣市之標準化死亡率資料檔,SDR為標準化死亡率。[: 因為各縣市人口結果不同,故以2000年世界標準人口調整後的標準化率作比較]


[程式二] 為不同等級的死亡率設定在地圖上的色塊,例如 pattern5 v=ms c=brown; 其中V=MS指明區塊為實心,c=brown指明區塊顏色為棕色。


[程式三] PROC FORMAT/VALUE介定區塊對應的標準化率值,並形成一個格式化名稱degfmt,任何一個變項套用這個格式,數值在300-400間的就被定義為'350-400',其他以此類推。


[程式四] proc gmap畫地圖,map=maps.taiwan說明引用的地圖經緯來自目錄館為MAPS下的資料檔TAIWAN,繪圖的來源檔為data=aa以各縣市(ID)標示區塊,並將將SDRdegfmt. 格式套用。


所以圖中出現白色與綠色區塊部份, 表示台北人的標準化死亡率最低,台東縣最高,其次花蓮縣與屏東縣。

[
程式一];

data aa;
input id city $ sdr/*--
標準化死亡率*/;
cards;
1
彰化縣 473.6
2
嘉義縣 506.6
3
基隆市 474.2
4
新竹縣 473.3
5
花蓮縣 578.8
6
宜蘭縣 478.8
7
高雄市 499.9
8
高雄市 499.9
9
苗栗縣 496.1
10
南投縣 532.3
11
澎湖縣 496.2
12
屏東縣 592.9
13
台中市 471.6
14
台中市 471.6
15
台南市 496.6
16
台南市 496.6
17
新北市 423.4
18
台北市 343.5
19
台東縣 639.2
20
桃園縣 450.2
21
雲林縣 525.7
;
proc print;
run;

[
程式二];
goptions reset=all;
pattern1 v=ms c=white;
pattern2 v=ms c=green;
pattern3 v=ms c=yellow;
pattern4 v=ms c=orange;
pattern5 v=ms c=brown;
pattern6 v=ms c=black;

[
程式三];
proc format;
value degfmt 
300-400='<400'
400-449.9='400-449.9'
450-499.9='450-499.9'
500-549.9='500-549.9'
550-599.9='550-599.9'
600-649.9='600-649.9';

  [
程式四];
proc gmap map=maps.taiwan data=aa all;
id id; /*-
以各縣市標示區塊*/
choro sdr/discrete coutline=black;
format sdr degfmt.; /*
SDRdegfmt. 格式套用*/
run;
quit;

 

[結果]

 

Taiwan

arrow
arrow

    estat 發表在 痞客邦 留言(0) 人氣()