廈門大學(xué)《應(yīng)用多元統(tǒng)計分析》

《廈門大學(xué)《應(yīng)用多元統(tǒng)計分析》》由會員分享，可在線閱讀，更多相關(guān)《廈門大學(xué)《應(yīng)用多元統(tǒng)計分析》（63頁珍藏版）》請在裝配圖網(wǎng)上搜索。

1、第四章 判別分析第一節(jié) 引言 第二節(jié) 距離判別法 第三節(jié) 貝葉斯（Bayes）判別法 第四節(jié) 費歇（Fisher）判別法 第五節(jié) 實例分析與計算機實現(xiàn) 第一節(jié) 引言n在我們的日常生活和工作實踐中，常常會遇到判別分析問題，即根據(jù)歷史上劃分類別的有關(guān)資料和某種最優(yōu)準則，確定一種判別方法，判定一個新的樣本歸屬哪一類。例如，某醫(yī)院有部分患有肺炎、肝炎、冠心病、糖尿病等病人的資料，記錄了每個患者若干項癥狀指標數(shù)據(jù)。現(xiàn)在想利用現(xiàn)有的這些資料找出一種方法，使得對于一個新的病人，當測得這些癥狀指標數(shù)據(jù)時，能夠判定其患有哪種病。又如，在天氣預(yù)報中，我們有一段較長時間關(guān)于某地區(qū)每天氣象的記錄資料（晴陰雨、氣溫、氣

2、壓、濕度等），現(xiàn)在想建立一種用連續(xù)五天的氣象資料來預(yù)報第六天是什么天氣的方法。這些問題都可以應(yīng)用判別分析方法予以解決。 n把這類問題用數(shù)學(xué)語言來表達，可以敘述如下：設(shè)有n個樣本，對每個樣本測得p項指標（變量）的數(shù)據(jù)，已知每個樣本屬于k個類別（或總體）G1，G2， ，Gk中的某一類，且它們的分布函數(shù)分別為F1(x)，F(xiàn)2(x)， ，F(xiàn)k(x)。我們希望利用這些數(shù)據(jù)，找出一種判別函數(shù)，使得這一函數(shù)具有某種最優(yōu)性質(zhì)，能把屬于不同類別的樣本點盡可能地區(qū)別開來，并對測得同樣p項指標（變量）數(shù)據(jù)的一個新樣本，能判定這個樣本歸屬于哪一類。 n 判別分析內(nèi)容很豐富，方法很多。判斷分析按判別的總體數(shù)來區(qū)分，有兩

3、個總體判別分析和多總體判別分析；按區(qū)分不同總體所用的數(shù)學(xué)模型來分，有線性判別和非線性判別；按判別時所處理的變量方法不同，有逐步判別和序貫判別等。判別分析可以從不同角度提出問題，因此有不同的判別準則，如馬氏距離最小準則、Fisher準則、平均損失最小準則、最小平方準則、最大似然準則、最大概率準則等等，按判別準則的不同又提出多種判別方法。本章僅介紹常用的幾種判別分析方法：距離判別法、Fisher判別法、Bayes判別法和逐步判別法。 第二節(jié) 距離判別法一 馬氏距離的概念 二 距離判別的思想及方法 三 判別分析的實質(zhì) n為此，我們引入一種由印度著名統(tǒng)計學(xué)家馬哈拉諾比斯（Mahalanobis, 19

4、36）提出的“馬氏距離”的概念。n 二、距離判別的思想及方法 1、兩個總體的距離判別問題n 問題：設(shè)有協(xié)方差矩陣相等的兩個總體G1和G2，其均值分別是1和 2，對于一個新的樣品X，要判斷它來自哪個總體。n 一般的想法是計算新樣品X到兩個總體的馬氏距離D2（X，G1）和D2（X，G2），并按照如下的判別規(guī)則進行判斷n這個判別規(guī)則的等價描述為：求新樣品X到G 1的距離與到G2的距離之差，如果其值為正，X屬于G2；否則X屬于G1。 n我們考慮 n n n n這里我們應(yīng)該注意到： 2、多個總體的距離判別問題n n n n 三、判別分析的實質(zhì)n我們知道，判別分析就是希望利用已經(jīng)測得的變量數(shù)據(jù)，找出一種判

5、別函數(shù)，使得這一函數(shù)具有某種最優(yōu)性質(zhì)，能把屬于不同類別的樣本點盡可能地區(qū)別開來。為了更清楚的認識判別分析的實質(zhì)，以便能靈活的應(yīng)用判別分析方法解決實際問題，我們有必要了解“劃分”這樣概念。n設(shè)R1，R2，Rk是p維空間R p的k個子集，如果它們互不 相交，且它們的和集為R p，則稱R1，R2， ，Rk為R p的一個劃分。 n n 這樣我們將會發(fā)現(xiàn)，判別分析問題實質(zhì)上就是在某種意義上，以最優(yōu)的性質(zhì)對p維空間R p構(gòu)造一個“劃分”，這個“劃分”就構(gòu)成了一個判別規(guī)則。這一思想將在后面的各節(jié)中體現(xiàn)的更加清楚。 第三節(jié) 貝葉斯（Bayes）判別法一 Bayes判別的基本思想 二 Bayes判別的基本方法

6、n從上節(jié)看距離判別法雖然簡單，便于使用。但是該方法也有它明顯的不足之處。第一，判別方法與總體各自出現(xiàn)的概率的大小無關(guān)；第二，判別方法與錯判之后所造成的損失無關(guān)。Bayes判別法就是為了解決這些問題而提出的一種判別方法。 一、Bayes判別的基本思想n n n n 二、Bayes判別的基本方法n n如果已知樣品X來自總體G i 的先驗概率為qi ， ，則在規(guī)則R下，由（4.12）式知，誤判的總平均損失為 ki ,2,1 n n n n 第四節(jié) 費歇（Fisher）判別法一 Fisher判別的基本思想 二 Fisher判別函數(shù)的構(gòu)造 三 線性判別函數(shù)的求法 n Fisher判別法是1936年提出來

7、的，該方法的主要思想是通過將多維數(shù)據(jù)投影到某個方向上，投影的原則是將總體與總體之間盡可能的放開，然后再選擇合適的判別規(guī)則，將新的樣品進行分類判別。 一、Fisher判別的基本思想n 二、Fisher判別函數(shù)的構(gòu)造1、針對兩個總體的情形n 2、針對多個總體的情形n n 三、線性判別函數(shù)的求法n n n n這里值得注意的是，本書有幾處利用極值原理求極值時，只給出了不要條件的數(shù)學(xué)推導(dǎo)，而有關(guān)充分條件的論證省略了，因為在實際問題中，往往根據(jù)問題本身的性質(zhì)就能肯定有最大值（或最小值），如果所求的駐點只有一個，這時就不需要根據(jù)極值存在的充分條件判定它是極大還是極小而就能肯定這唯一的駐點就是所求的最大值（或

8、最小值）。為了避免用較多的數(shù)學(xué)知識或數(shù)學(xué)上的推導(dǎo)，這里不追求數(shù)學(xué)上的完整性。n 第五節(jié) 實例分析與計算機實現(xiàn)n這一節(jié)我們利用SPSS對Fisher判別法和Bayes判別法進行計算機實現(xiàn)。n為研究某地區(qū)人口死亡狀況，已按某種方法將15個已知地區(qū)樣品分為3類，指標含義及原始數(shù)據(jù)如下。試建立判別函數(shù)，并判定另外4個待判地區(qū)屬于哪類？ 表4.1 各地區(qū)死亡概率表 (一) 操作步驟1. 在SPSS窗口中選擇AnalyzeClassifyDiscriminate，調(diào)出判別分析主界面，將左邊的變量列表中的“group”變量選入分組變量中，將變量選入自變量中，并選擇Enter independents tog

9、ether單選按鈕，即使用所有自變量進行判別分析。圖4.2 判別分析主界面 2. 點擊Define Range按鈕，定義分組變量的取值范圍。本例中分類變量的范圍為1到3，所以在最小值和最大值中分別輸入1和3。單擊Continue按鈕，返回主界面。3. 單擊Statistics按鈕，指定輸出的描述統(tǒng)計量和判別函數(shù)系數(shù)。選中Function Coefficients欄中的Fishers和Unstandardized。這兩個選項的含義如下：Fishers：給出Bayes判別函數(shù)的系數(shù)。（注意：這個選項不是要給出Fisher判別函數(shù)的系數(shù)。這個復(fù)選框的名字之所以為Fishers，是因為按判別函數(shù)值最大

10、的一組進行歸類這種思想是由Fisher提出來的。這里極易混淆，請讀者注意辨別。）Unstandardized：給出未標準化的Fisher判別函數(shù)（即典型判別函數(shù)）的系數(shù)（SPSS默認給出標準化的Fisher判別函數(shù)系數(shù)）。 n單擊Continue按鈕，返回主界面。圖4.3 Statistics子對話框輸出Bayes 判別系數(shù)標準化的費希爾判別系數(shù) 4. 單擊Classify按鈕，定義判別分組參數(shù)和選擇輸出結(jié)果。選擇Display欄中的Casewise results，輸出一個判別結(jié)果表，包括每個樣品的判別分數(shù)、后驗概率、實際組和預(yù)測組編號等。其余的均保留系統(tǒng)默認選項。單擊Continue按鈕。

11、圖4.4 Classify子對話框各組概率相等根據(jù)各組樣品數(shù)確定先驗概率使用聯(lián)合協(xié)方差矩陣，此選項表明各總體協(xié)方差矩陣相等各個總體協(xié)方差矩陣不等 5. 單擊Save按鈕，指定在數(shù)據(jù)文件中生成代表判別分組結(jié)果和判別得分的新變量，生成的新變量的含義分別為：Predicted group membership：存放判別樣品所屬組別的值； Discriminant scores：存放Fisher判別得分的值，有幾個典型判別函數(shù)就有幾個判別得分變量；Probabilities of group membership：存放樣品屬于各組的Bayes后驗概率值。n將對話框中的三個復(fù)選框均選中，單擊Contin

12、ue按鈕返回。 6. 返回判別分析主界面，單擊OK按鈕，運行判別分析過程。圖4.5 Save子對話框各組先驗概率未標準化的費希爾判別得分成員屬于各組的后驗概率 Standardized Canonical Discriminant Function （二） 主要運行結(jié)果解釋1. Standardized Canonical Discriminant Function Coefficients（給出標準化的典型判別函數(shù)系數(shù)）標準化的典型判別函數(shù)是由標準化的自變量通過Fisher判別法得到的，所以要得到標準化的典型判別得分，代入該函數(shù)的自變量必須是經(jīng)過標準化的。2. Canonical Discr

13、iminant Function Coefficients（給出未標準化的典型判別函數(shù)系數(shù)）未標準化的典型判別函數(shù)系數(shù)由于可以將實測的樣品觀測值直接代入求出判別得分，所以該系數(shù)使用起來比標準化的系數(shù)要方便一些。見表4.2（a）。 *2*2 *2*1 76543 76543 791.1709.0689.0499.0746.9575.7 100.4284.1173.6232.1969.13268.16 XXXXXXy XXXXXXy Standardized Canonical Discriminant Function 標準化的費希爾判別函數(shù) 2 Canonical Discriminant F

14、unction Coefficients（給出未標準化的典型判別函數(shù)系數(shù)） 1 1 2 3 4 5 6 2 1 2 3 4 5 6 74.99 1.861 1.656 0.877 0.798 0.098 1.579 29.482 0.867 1.155 0.356 0.089 0.054 0.69 y X X X X X X y X X X X X X 未標準化的費希爾判別函數(shù)典型判別函數(shù)并非費希爾判別函數(shù)，二者相差一個常數(shù) )756.0,600.6(,)257.0,194.9(,)013.1,594.2( 321 GGG Functions at Group Centroids（給出組重心處

15、的Fisher判別函數(shù)值）各個組對應(yīng)的典型判別函數(shù)的質(zhì)心坐標， 即 基于費希爾判別函數(shù)的判別規(guī)則本教材只是給出了費希爾判別函數(shù)的求法，并未給出具體的判別準則，以下給出兩個可行的準則1. 計算各個樣本點到各組質(zhì)心處的歐氏距離， 根據(jù)樣品到各個總體歐氏距離的大小判定它屬于哪個總體2. 對于各個樣品的費希爾判別函數(shù)，以及各個組的組別，進一步進行bayes判別。 Bayes 判別函數(shù)系數(shù) Classification Function Coefficients -143.851 -164.691 -134.862 153.137 171.185 144.462 -90.088 -99.976 -85.

16、945 53.009 62.525 49.972 11.008 12.094 10.520 189.261 207.003 181.714 -5317.234 -6202.158 -4982.880 X1 X2 X3 X4 X5 X6 (Constant) 1.00 2.00 3.00 GROUP Fishers linear discriminant functions 111121111 ln)()(F PxSxxxSx 11121 ln- PxSx 其中， S是聯(lián)合協(xié)方差矩陣， P1是G1的先驗概率- 4. Classification Function Coefficients（給出B

17、ayes判別函數(shù)系數(shù)）如表4.3所示，GROUP欄中的每一列表示樣品判入相應(yīng)列的Bayes判別函數(shù)系數(shù)。在本例中，各類的Bayes判別函數(shù)如下：第一組：第二組：第三組： 1 1 2 3 4 5 65317.2 143.9 153.1 90.1 53.0 11.0 189.3F X X X X X X 2 1 2 3 4 5 66202.2 164.7 171.2 100.0 62.5 12.1 207.0F X X X X X X 3 1 2 3 4 5 64982.9 134.9 144.5 85.9 50.0 10.5 181.7F X X X X X X 將各樣品的自變量值代入上述三個B

18、ayes判別函數(shù)，得到三個函數(shù)值。比較這三個函數(shù)值，哪個函數(shù)值比較大就可以判斷該樣品判入哪一類。例如，將第一個待判樣品的自變量值分別代入函數(shù)，得到： F1=3793.77， F2=3528.32， F3=3882.48比較三個值，可以看出 F3 最大，據(jù)此得出第一個待判樣品應(yīng)該屬于第三組。 5. Casewise Statistics（給出個案觀察結(jié)果）在Casewise Statistics輸出表針對每個樣品給出了了大部分的判別結(jié)果，其中包括：實際類（Actual Group）、預(yù)測類（Predicted Group）、Bayes判別法的后驗概率、與組重心的馬氏距離（Squared Maha

19、lanobis Distance to Centroid）以及Fisher判別法的每個典型判別函數(shù)的判別得分（Discriminant Scores）。出于排版要求，這里給出結(jié)果表的是經(jīng)過加工的，隱藏了其中的一些項目，如表4.4所示。從表中可以看出四個待判樣本依次被判別為第三組、第一組、第二組和第三組。 Casewise Statistics 1 1 1.000 .297 -2.177 1.364 1 1 1.000 .236 -2.270 1.375 1 1 1.000 .117 -2.741 1.323 1 1 .998 .507 -3.199 .638 1 1 1.000 .418 -2

20、.582 .366 2 2 1.000 .469 9.674 .231 2 2 1.000 .868 8.332 -.613 2 2 1.000 5.985 10.128 -2.518 2 2 1.000 4.793 8.342 1.760 2 2 1.000 .101 9.491 -.145 3 3 1.000 .139 -6.687 -.394 3 3 1.000 .322 -7.163 -.685 3 3 1.000 5.365 -8.655 -1.823 3 3 .879 3.384 -4.766 -.608 3 3 .995 .998 -5.727 -.270 ungrouped 3

21、 1.000 361.567 -20.714 -13.498 ungrouped 1 .998 .558 -3.319 .831 ungrouped 2 1.000 28.668 14.008 2.086 ungrouped 3 1.000 1.982 -7.595 -1.752 Case Number 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 Actual Group Predicted Group P(G=g | D=d) Squared Mahalanobis Distance to Centroid Highest Group Fu

22、nction 1 Function 2 Discriminant Scores 表4.4 個案觀察結(jié)果表實際組預(yù)測組到距離它最近一組，即到預(yù)測組的廣義馬氏距離屬于該組的后驗概率判別得分，由于我們設(shè)定非標準化的費希爾判別，此處是非標準的費希爾判別得分，若無此選擇，則顯示標準化的費希爾判別得分 6. 由于我們在Save子對話框中選擇了生成表示判別結(jié)果的新變量，所以在數(shù)據(jù)編輯窗口中，可以觀察到產(chǎn)生的新變量。其中，變量dis-1存放判別樣品所屬組別的值，變量dis1-1和dis2-1分別代表將樣品各變量值代入第一個和第二個判別函數(shù)所得的判別分數(shù)，變量dis1-2、dis2-2和dis3-2分別代表樣品分別屬于第1組、第2組和第3組的Bayes后驗概率值。 本 章 結(jié) 束