聚類(lèi)分析的目的是把分類(lèi)對(duì)象按一定的規(guī)則分成若干類(lèi),,這些類(lèi)不是事先給定的,,而是根據(jù)數(shù)據(jù)的特征確定的,對(duì)類(lèi)的數(shù)目和類(lèi)的結(jié)構(gòu)不必作任何的假定,。在同一類(lèi)里的這些對(duì)象在某種意義上傾向于彼此相似,,而在不同類(lèi)里的對(duì)象傾向于彼此不相似。
聚類(lèi)分析根據(jù)分類(lèi)對(duì)象不同分為Q型聚類(lèi)分析和R型聚類(lèi)分析,。Q型聚類(lèi)是指對(duì)樣品進(jìn)行聚類(lèi),,R型聚類(lèi)分析是指對(duì)變量進(jìn)行聚類(lèi)。 SAS中用于判別分析的過(guò)程主要有以下四個(gè): PROC CLUSTER 譜系聚類(lèi)(Q型聚類(lèi)分析) PROC FASTCLUS K均值快速聚類(lèi),,適用于大樣本(Q型聚類(lèi)分析) PROC MODECLUS 非參數(shù)聚類(lèi)(Q型聚類(lèi)分析) PROC VARCLUS 變量聚類(lèi)(R型聚類(lèi)分析) 1. 用PROC CLUSTER進(jìn)行樣品聚類(lèi)分析(Q型聚類(lèi)分析) [例1] 試根據(jù)1997年信息基礎(chǔ)設(shè)施的發(fā)展情況,,對(duì)世界20個(gè)國(guó)家和地區(qū)進(jìn)行聚類(lèi)分析。描述信息基礎(chǔ)設(shè)施的變量主要有六個(gè): Cal-每千人擁有電話線數(shù),; Cellphone-每千戶(hù)居民移動(dòng)電話數(shù),; Fee-高峰時(shí)期每三分鐘國(guó)際電話成本; Computer-每千人擁有的計(jì)算機(jī)數(shù),; Mips-每千人中計(jì)算機(jī)功率(每秒百萬(wàn)指令),; Net-每千人互聯(lián)網(wǎng)絡(luò)戶(hù)主數(shù)。 data cluster1; infile datalines; input Country $ Cal Cellphone Fee Computer Mips Net; datalines; 美國(guó) 631.6 161.9 0.36 403 26073 35.34 日本 498.4 143.2 3.57 176 10223 6.26 德國(guó) 557.6 70.6 2.18 199 11571 9.48 瑞典 684.1 281.8 1.4 286 16660 29.39 瑞士 644 93.5 1.98 234 13621 22.68 丹麥 620.3 248.6 2.56 296 17210 21.84 新加坡 498.4 147.5 2.5 284 13578 13.49 中國(guó)臺(tái)灣 469.4 56.1 3.68 119 6911 1.72 韓國(guó) 434.5 73 3.36 99 5795 1.66 巴西 81.9 16.3 3.02 19 876 0.52 智利 138.6 8.2 1.4 31 1411 1.28 墨西哥 92.2 9.8 2.61 31 1751 0.35 俄羅斯 174.9 5 5.12 24 1101 0.48 波蘭 169 6.5 3.68 40 1796 1.45 匈牙利 262.2 49.4 2.66 68 3067 3.09 馬來(lái)西亞 195.5 88.4 4.19 53 2734 1.25 泰國(guó) 78.6 27.8 4.95 22 1662 0.11 印度 13.6 0.3 6.28 2 101 0.01 法國(guó) 559.1 42.9 1.27 201 11702 4.76 英國(guó) 521.1 122.5 0.98 248 14461 11.91 ; run; PROC CLUSTER DATA=cluster1 STANDARD METHOD=CENTROID CCC PSEUDO UT=TREE; PROC TREE DATA=TREE HORIZONTAL SPACES=1; RUN; [說(shuō)明] METHOD=的選項(xiàng)可以為: AVERAGE(平均法) CENTROID(重心法) COMPLETE(最長(zhǎng)距離法) DENSITY(非參數(shù)概率密度估計(jì)法) EML(最大似然法) FLEXIBLE(flexible-beta法) MCQUITTY(Mcquitty的相似分析法) MEDIAN(中位數(shù)法) SINGLE(最短距離法) TWOSTAGE(兩階段密度法) WARD(Ward最小方差法) STANDARD 對(duì)變量實(shí)施標(biāo)準(zhǔn)化,。 CCC,、PSEUDO 為了計(jì)算一些統(tǒng)計(jì)量用以判別全部樣品究竟聚成幾類(lèi)較為合適。CCC要求打印聚類(lèi)判別據(jù)的立方及在一致無(wú)效假設(shè)下近似期望值R2,,PSEUDO要求打印偽F(標(biāo)志PSF)和t2(標(biāo)志PST2)統(tǒng)計(jì)量,。當(dāng)分類(lèi)數(shù)目不同時(shí),它們就有不同的取值,,CCC和PSF出現(xiàn)峰值所對(duì)應(yīng)的分類(lèi)數(shù)較合適,、PST2出現(xiàn)峰值的前一行所對(duì)應(yīng)的分類(lèi)數(shù)較合適。 OUT=TREE 產(chǎn)生名為T(mén)REE的輸出數(shù)據(jù)集,它可被TREE過(guò)程用來(lái)輸出聚類(lèi)結(jié)果的樹(shù)狀圖,。HORIZONTAL要求將樹(shù)狀圖水平放置,,SPACES=1要求各樣品之間的間隔為1。 [結(jié)果及其解釋] 如果聚為3類(lèi),,其聚類(lèi)結(jié)果為:一類(lèi)是信息基礎(chǔ)設(shè)施最為發(fā)達(dá)的美國(guó),,一類(lèi)是一些發(fā)達(dá)國(guó)家,其他的國(guó)家和地區(qū)聚為另外一類(lèi),。 2. 用PROC FASTCLUS進(jìn)行大樣本的樣品聚類(lèi)分析(Q型聚類(lèi)分析) 處理大樣本時(shí)一般采用非分層聚類(lèi)法(快速聚類(lèi)法),。聚類(lèi)的個(gè)數(shù)k可以根據(jù)需要事先指定。與分層聚類(lèi)方法相比,,非分層聚類(lèi)方法不必確定距離矩陣,,不必存儲(chǔ)基本數(shù)據(jù),因此適用于處理很大的數(shù)據(jù)集,。 下面仍使用[例1]的數(shù)據(jù),,對(duì)PROC FASTCLUS加以說(shuō)明: PROC FASTCLUS DATA=cluster1 UT=result MAXC=3 CLUSTER=c; RUN; [說(shuō)明] MAXC=3 要求總共聚為3類(lèi) OUT=result 結(jié)果輸出到數(shù)據(jù)集result CLUSTER=c 分類(lèi)標(biāo)志的變量名c,其取值為1,,2,,3 [結(jié)果及其解釋] 聚類(lèi)結(jié)果為:一類(lèi)是信息基礎(chǔ)設(shè)施最為發(fā)達(dá)的美國(guó),一類(lèi)是一些發(fā)達(dá)國(guó)家,,其他的國(guó)家和地區(qū)聚為另外一類(lèi),。 3. 用 PROC VARCLUS 進(jìn)行變量聚類(lèi)分析(R型聚類(lèi)分析) [例 2] 對(duì)1996年全國(guó)30個(gè)省市自治區(qū)經(jīng)濟(jì)發(fā)展基本情況的八項(xiàng)指標(biāo)作變量聚類(lèi)分析: X1:GDP X2:居民消費(fèi)水平 X3:固定資產(chǎn)投資 X4:職工平均工資 X5:貨物周轉(zhuǎn)量 X6:居民消費(fèi)價(jià)格指數(shù) X7:商品零售價(jià)格指數(shù) X8:工業(yè)總產(chǎn)值 data cluster2; infile datalines; input Province $ X1-X8; datalines; 北京 1394.89 2505 519.01 8144 373.9 117.3 112.6 843.43 天津 920.11 2720 345.46 6501 342.8 115.2 110.6 582.51 河北 2849.52 1258 704.87 4839 2033.3 115.2 115.8 1234.85 山西 1092.48 1250 290.9 4721 717.3 116.9 115.6 697.25 內(nèi)蒙 832.88 1387 250.23 4134 781.7 117.5 116.8 419.39 遼寧 2793.37 2397 387.99 4911 1371.1 116.1 114 1840.55 吉林 1129.2 1872 320.45 4430 497.4 115.2 114.2 762.47 黑龍江 2014.53 2334 435.73 4145 824.8 116.1 114.3 1240.37 上海 2462.57 5343 996.48 9279 207.1 118.7 113 1642.95 江蘇 5155.25 1926 1434.95 5943 1025.5 115.8 114.3 2026.64 浙江 3524.79 2249 1006.39 6619 754.4 116.6 113.5 916.59 安徽 2003.58 1254 474 4609 908.3 114.8 112.7 824.14 福建 2160.52 2320 553.97 5857 609.3 115.2 114.4 433.67 江西 1205.11 1182 282.84 4211 411.7 116.9 115.9 571.84 山東 5002.34 1527 1229.55 5145 1196.6 117.6 114.2 2207.69 河南 3002.74 1034 670.35 4344 1574.4 116.5 114.9 1367.92 湖北 2391.42 1527 571.86 4685 849 120 116.6 1220.72 湖南 2195.7 1408 422.61 4797 1011.8 119 115.5 843.83 廣東 5381.72 2699 1639.83 8250 656.5 114 111.6 1396.35 廣西 1606.15 1314 382.59 5105 556 118.4 116.4 554.97 海南 364.17 1814 198.35 5340 232.1 113.5 111.3 64.33 四川 3534 1261 822.54 4645 902.3 118.5 117 1431.81 貴州 630.07 942 150.84 4475 301.1 121.4 117.2 324.72 云南 1206.68 1261 334 5149 310.4 121.3 118.1 716.65 西藏 55.98 1110 17.87 7382 4.2 117.3 114.9 5.57 陜西 1000.03 1208 300.27 4396 500.9 119 117 600.98 甘肅 553.35 1007 114.81 5493 507 119.8 116.5 468.79 青海 165.31 1445 47.76 5753 61.6 118 116.3 105.8 寧夏 169.75 1355 61.98 5079 121.8 117.1 115.3 114.4 新疆 834.57 1469 376.95 5348 339 119.7 116.7 428.76 ; run; PROC VARCLUS DATA=cluster2 CENTROID MAXC=3; VAR x1-x8; RUN; [說(shuō)明] CENTROID 聚類(lèi)方法為重心法,默認(rèn)聚類(lèi)方法為主成分法 MAXC=3 要求總共聚為3類(lèi) [結(jié)果及其解釋] 聚類(lèi)結(jié)果為:第一類(lèi)變量主要反映了生產(chǎn)狀況,;第二類(lèi)變量主要反映了消費(fèi)狀況,,第三類(lèi)變量主要反映了價(jià)格狀況。 window.google_render_ad(); |
|