1、主成分與聚類分析首先通過SPSS軟件對環境污染的相應指標進行主成分分析，得到：提取Y1、Y2、Y3和Y4四個主成分，其累積貢獻率已經達到，超過80%，代表所有環境污染指標的絕大部分信息。Y1偏向于解釋工業氫氧化物排放量，Y2偏向于解釋生活煙塵排放量，Y3偏向于解釋生活廢水排放量，Y4偏向于解釋工業二氧化硫排放量。然后，根據主成分分析結果，用Z=*Y1+*Y2+*Y3+0.06519*Y4計算綜合得分，見下表1。表1環境污染地區的主成分綜合得分表序號地區Z排名序號地區Z排名1北京517武漢132天津418長沙283石家莊619廣州194太原820南寧245呼和浩特1221海口316沈陽1722重

2、慶17長春1623成都208哈爾濱224貴陽189上海325昆明2610南京1526拉薩3011杭州927西安712合肥2128蘭州2313福州2529西寧1114南昌2930銀川2715濟南1031烏魯木齊2216鄭州14最后將環境污染的綜合得分作為個案進行層次聚類分析，將31個地區分為5類，如表2。表2 各地區污染分類分類污染情況地區1輕度污染海口、拉薩2比較輕度污染合肥、烏魯木齊、福州、南寧、蘭州、，昆明、成都、銀川、南昌、長沙、沈陽、長春、南京、廣州、貴陽、鄭州、武漢、濟南、西寧、呼和浩特3污染情況一般太原、杭州、石家莊、西安4污染比較嚴重北京、天津5污染十分嚴重上海、哈爾濱、重慶主成

3、分分析和聚類分析在SPSS中的操作過程打開SPSS，“文件-打開-數據”，選中excel，如下圖結果。首先將變量標準化，“分析-描述統計-描述”，將變量全部選入對話框，點上“將標準化得分另存為變量（Z）”，結果如下。在做主成分分析，“分析-降維-因子分析”，將為標準化的變量選入對話框。選擇右側“描述”，在彈出來對話框中點上“系數、顯著性水平、KMO和Bartlett的球形度檢驗”，點“繼續”。選擇右側“抽取”，在彈出來對話框中將方法改成“主成分”，“基于特征值”改成“0.6”，點“繼續”。選擇右側“旋轉”，在彈出來對話框中，將方法改成“最大方差法”，點“繼續”。首先輸出的各自變量的相關矩陣，由

4、下圖可見，各自變量間存在顯著性相關。說明有必要進行因子分析。接著輸出的是KMO和Bartlett的檢驗，結果如下圖，由表可見，KMO是Kaiser-Meyer-Olkin的取樣適當性度量，當KM0越大時，表示變量間的共同因素越多，越適合進行因子分析，此處KMO大于，表示可進行因子分析。此外，Bartlett球形檢驗的近似卡方分布為(自由度為45)達到顯著，代表有共同因素存在，適合進行因素分析。同時，Bartlett球度檢驗給出的相伴概率為，小于顯著性水平，故適合進行因子分析。隨后輸出的解釋的總方差，由下圖可見，主成分個數為4時，累積貢獻率已經達到82.037%，超過80%。故本次因子分析中提取

5、4個因子，我們分別稱之為F1，F2，F3，F4。在下面的成分矩陣圖中可以看出，F1主要是和工業氫氧化物排放量有關，F2主要和生活煙塵排放量有關，F3主要和生活廢水排放量有關，F4主要和工業二氧化硫排放量有關。隨后將成分矩陣表中的系數復制到數據窗口，并將變量名命名為a1，a2，a3，a4。再在數據窗口點擊“轉換-計算變量”，將目標變量t1=a1/SQRT(4.323)。（括號內為對應主成分的初始特征值），依次算出t2，t3，t4。則可知道Y1=0.382*Z工業廢水排放量*Z工業化學需氧量排放量*Z工業二氧化硫排放量*Z工業煙塵排放量*工業氫氧化物排放量*Z生活廢水排放量+0.357*Z生活化學需氧量排放量+0.320*Z生活二氧化硫排放量+0.136*Z生活煙塵排放量+0.261*Z生活氫氧化物排放量，依次依據t2，t3，t4得到Y2，Y3，Y4。再依據Z=0.43226*Y1+0.21911*Y2+0.10380*Y3+0.06519*Y4得到綜合得分，見表1。、然后做層次聚類分析，“分析-分類-系統-系統聚類”，在彈出來對話框把綜合得分選入變量欄，在右側“繪制”中點上“樹狀圖”，點確定。如下圖。依據下面的樹狀圖，可將序號12,31,13,20,28,25,23,30,14,18,6,7,10,