廈門大學(xué)《應(yīng)用多元統(tǒng)計(jì)分析》第04章判別分析課件

第四章第四章 判別分析判別分析第一節(jié)第一節(jié) 引言引言 第二節(jié)第二節(jié) 距離判別法距離判別法 第三節(jié)第三節(jié) 貝葉斯（貝葉斯（Bayes）判別法）判別法 第四節(jié)第四節(jié) 費(fèi)歇（費(fèi)歇（Fisher）判別法）判別法 第五節(jié)第五節(jié) 實(shí)例分析與計(jì)算機(jī)實(shí)現(xiàn)實(shí)例分析與計(jì)算機(jī)實(shí)現(xiàn) 第一節(jié)第一節(jié) 引言引言n在我們的日常生活和工作實(shí)踐中，常常會(huì)遇到判別分析問題，在我們的日常生活和工作實(shí)踐中，常常會(huì)遇到判別分析問題，即根據(jù)歷史上劃分類別的有關(guān)資料和某種最優(yōu)準(zhǔn)則，確定一即根據(jù)歷史上劃分類別的有關(guān)資料和某種最優(yōu)準(zhǔn)則，確定一種判別方法，判定一個(gè)新的樣本歸屬哪一類。例如，某醫(yī)院種判別方法，判定一個(gè)新的樣本歸屬哪一類。例如，某醫(yī)院有部分患有肺炎、肝炎、冠心病、糖尿病等病人的資料，記有部分患有肺炎、肝炎、冠心病、糖尿病等病人的資料，記錄了每個(gè)患者若干項(xiàng)癥狀指標(biāo)數(shù)據(jù)?，F(xiàn)在想利用現(xiàn)有的這些錄了每個(gè)患者若干項(xiàng)癥狀指標(biāo)數(shù)據(jù)。現(xiàn)在想利用現(xiàn)有的這些資料找出一種方法，使得對(duì)于一個(gè)新的病人，當(dāng)測(cè)得這些癥資料找出一種方法，使得對(duì)于一個(gè)新的病人，當(dāng)測(cè)得這些癥狀指標(biāo)數(shù)據(jù)時(shí)，能夠判定其患有哪種病。又如，在天氣預(yù)報(bào)狀指標(biāo)數(shù)據(jù)時(shí)，能夠判定其患有哪種病。又如，在天氣預(yù)報(bào)中，我們有一段較長時(shí)間關(guān)于某地區(qū)每天氣象的記錄資料中，我們有一段較長時(shí)間關(guān)于某地區(qū)每天氣象的記錄資料（晴陰雨、氣溫、氣壓、濕度等），現(xiàn)在想建立一種用連續(xù)（晴陰雨、氣溫、氣壓、濕度等），現(xiàn)在想建立一種用連續(xù)五天的氣象資料來預(yù)報(bào)第六天是什么天氣的方法。這些問題五天的氣象資料來預(yù)報(bào)第六天是什么天氣的方法。這些問題都可以應(yīng)用判別分析方法予以解決。都可以應(yīng)用判別分析方法予以解決。n把這類問題用數(shù)學(xué)語言來表達(dá)，可以敘述如下：設(shè)有把這類問題用數(shù)學(xué)語言來表達(dá)，可以敘述如下：設(shè)有n個(gè)樣個(gè)樣本，對(duì)每個(gè)樣本測(cè)得本，對(duì)每個(gè)樣本測(cè)得p項(xiàng)指標(biāo)（變量）的數(shù)據(jù)，已知每個(gè)樣項(xiàng)指標(biāo)（變量）的數(shù)據(jù)，已知每個(gè)樣本屬于本屬于k個(gè)類別（或總體）個(gè)類別（或總體）G1，G2， ，Gk中的某一類，且中的某一類，且它們的分布函數(shù)分別為它們的分布函數(shù)分別為F1(x)，F(xiàn)2(x)， ，F(xiàn)k(x)。我們希望。我們希望利用這些數(shù)據(jù)，找出一種判別函數(shù)，使得這一函數(shù)具有某種利用這些數(shù)據(jù)，找出一種判別函數(shù)，使得這一函數(shù)具有某種最優(yōu)性質(zhì)，能把屬于不同類別的樣本點(diǎn)盡可能地區(qū)別開來，最優(yōu)性質(zhì)，能把屬于不同類別的樣本點(diǎn)盡可能地區(qū)別開來，并對(duì)測(cè)得同樣并對(duì)測(cè)得同樣p項(xiàng)指標(biāo)（變量）數(shù)據(jù)的一個(gè)新樣本，能判定項(xiàng)指標(biāo)（變量）數(shù)據(jù)的一個(gè)新樣本，能判定這個(gè)樣本歸屬于哪一類。這個(gè)樣本歸屬于哪一類。n 判別分析內(nèi)容很豐富，方法很多。判斷分析按判別的總體判別分析內(nèi)容很豐富，方法很多。判斷分析按判別的總體數(shù)來區(qū)分，有兩個(gè)總體判別分析和多總體判別分析；按區(qū)分?jǐn)?shù)來區(qū)分，有兩個(gè)總體判別分析和多總體判別分析；按區(qū)分不同總體所用的數(shù)學(xué)模型來分，有線性判別和非線性判別；不同總體所用的數(shù)學(xué)模型來分，有線性判別和非線性判別；按判別時(shí)所處理的變量方法不同，有逐步判別和序貫判別等。按判別時(shí)所處理的變量方法不同，有逐步判別和序貫判別等。判別分析可以從不同角度提出問題，因此有不同的判別準(zhǔn)則，判別分析可以從不同角度提出問題，因此有不同的判別準(zhǔn)則，如馬氏距離最小準(zhǔn)則、如馬氏距離最小準(zhǔn)則、Fisher準(zhǔn)則、平均損失最小準(zhǔn)則、最準(zhǔn)則、平均損失最小準(zhǔn)則、最小平方準(zhǔn)則、最大似然準(zhǔn)則、最大概率準(zhǔn)則等等，按判別準(zhǔn)小平方準(zhǔn)則、最大似然準(zhǔn)則、最大概率準(zhǔn)則等等，按判別準(zhǔn)則的不同又提出多種判別方法。本章僅介紹常用的幾種判別則的不同又提出多種判別方法。本章僅介紹常用的幾種判別分析方法：距離判別法、分析方法：距離判別法、Fisher判別法、判別法、Bayes判別法和逐判別法和逐步判別法。步判別法。第二節(jié)第二節(jié) 距離判別法距離判別法一一 馬氏距離的概念馬氏距離的概念 二二 距離判別的思想及方法距離判別的思想及方法 三三 判別分析的實(shí)質(zhì)判別分析的實(shí)質(zhì) 一、馬氏距離的概念一、馬氏距離的概念n 圖圖4.1n為此，我們引入一種由印度著名統(tǒng)計(jì)學(xué)家馬哈拉諾比斯為此，我們引入一種由印度著名統(tǒng)計(jì)學(xué)家馬哈拉諾比斯（Mahalanobis, 1936）提出的）提出的“馬氏距離馬氏距離”的概念。的概念。n 二、距離判別的思想及方法二、距離判別的思想及方法 1、兩個(gè)總體的距離判別問題、兩個(gè)總體的距離判別問題n 問題：設(shè)有協(xié)方差矩陣問題：設(shè)有協(xié)方差矩陣相等的兩個(gè)總體相等的兩個(gè)總體G1和和G2，其均值，其均值分別是分別是 1和和 2，對(duì)于一個(gè)新的樣品，對(duì)于一個(gè)新的樣品X，要判斷它來自哪個(gè)總，要判斷它來自哪個(gè)總體。體。n 一般的想法是計(jì)算新樣品一般的想法是計(jì)算新樣品X到兩個(gè)總體的馬氏距離到兩個(gè)總體的馬氏距離D2（X，G1）和和D2（X，G2），并按照如下的判別規(guī)則進(jìn)行判斷，并按照如下的判別規(guī)則進(jìn)行判斷n這個(gè)判別規(guī)則的等價(jià)描述為：求新樣品這個(gè)判別規(guī)則的等價(jià)描述為：求新樣品X到到G1的距離與到的距離與到G2的距離之差，如果其值為正，的距離之差，如果其值為正，X屬于屬于G2；否則；否則X屬于屬于G1。 n我們考慮我們考慮 n n n n這里我們應(yīng)該注意到：這里我們應(yīng)該注意到：2、多個(gè)總體的距離判別問題、多個(gè)總體的距離判別問題n n n n 三、判別分析的實(shí)質(zhì)三、判別分析的實(shí)質(zhì)n我們知道，判別分析就是希望利用已經(jīng)測(cè)得的變量數(shù)據(jù)，找我們知道，判別分析就是希望利用已經(jīng)測(cè)得的變量數(shù)據(jù)，找出一種判別函數(shù)，使得這一函數(shù)具有某種最優(yōu)性質(zhì)，能把屬出一種判別函數(shù)，使得這一函數(shù)具有某種最優(yōu)性質(zhì)，能把屬于不同類別的樣本點(diǎn)盡可能地區(qū)別開來。為了更清楚的認(rèn)識(shí)于不同類別的樣本點(diǎn)盡可能地區(qū)別開來。為了更清楚的認(rèn)識(shí)判別分析的實(shí)質(zhì)，以便能靈活的應(yīng)用判別分析方法解決實(shí)際判別分析的實(shí)質(zhì)，以便能靈活的應(yīng)用判別分析方法解決實(shí)際問題，我們有必要了解問題，我們有必要了解“劃分劃分”這樣概念。這樣概念。n設(shè)設(shè)R1，R2，Rk是是p維空間維空間R p的的k個(gè)子集，如果它們互不個(gè)子集，如果它們互不 相交，且它們的和集為相交，且它們的和集為R p，則稱，則稱R1，R2， ，Rk為為R p的一的一個(gè)劃分。個(gè)劃分。n n 這樣我們將會(huì)發(fā)現(xiàn)，判別分析問題實(shí)質(zhì)上就是在某種意義這樣我們將會(huì)發(fā)現(xiàn)，判別分析問題實(shí)質(zhì)上就是在某種意義上，以最優(yōu)的性質(zhì)對(duì)上，以最優(yōu)的性質(zhì)對(duì)p維空間維空間R p構(gòu)造一個(gè)構(gòu)造一個(gè)“劃分劃分”，這個(gè)，這個(gè)“劃分劃分”就構(gòu)成了一個(gè)判別規(guī)則。這一思想將在后面的各節(jié)就構(gòu)成了一個(gè)判別規(guī)則。這一思想將在后面的各節(jié)中體現(xiàn)的更加清楚。中體現(xiàn)的更加清楚。 第三節(jié)第三節(jié) 貝葉斯（貝葉斯（Bayes）判別法）判別法一一 Bayes判別的基本思想判別的基本思想 二二 Bayes判別的基本方法判別的基本方法 n從上節(jié)看距離判別法雖然簡單，便于使用。但是該方法也有從上節(jié)看距離判別法雖然簡單，便于使用。但是該方法也有它明顯的不足之處。它明顯的不足之處。第一，判別方法與總體各自出現(xiàn)的概率的大小無關(guān)；第一，判別方法與總體各自出現(xiàn)的概率的大小無關(guān)；第二，判別方法與錯(cuò)判之后所造成的損失無關(guān)。第二，判別方法與錯(cuò)判之后所造成的損失無關(guān)。Bayes判別判別法就是為了解決這些問題而提出的一種判別方法。法就是為了解決這些問題而提出的一種判別方法。一、一、Bayes判別的基本思想判別的基本思想n n n n 二、二、Bayes判別的基本方法判別的基本方法n n如果已知樣品如果已知樣品X來自總體來自總體Gi 的先驗(yàn)概率為的先驗(yàn)概率為qi ， ，則在規(guī)則則在規(guī)則R下，由（下，由（4.12）式知，誤判的總平均損失為）式知，誤判的總平均損失為 ki, 2, 1n n n n 第四節(jié)第四節(jié) 費(fèi)歇（費(fèi)歇（Fisher）判別法）判別法一一 Fisher判別的基本思想判別的基本思想 二二 Fisher判別函數(shù)的構(gòu)造判別函數(shù)的構(gòu)造 三三 線性判別函數(shù)的求法線性判別函數(shù)的求法 nFisher判別法是判別法是1936年提出來的，該方法的主要思想是通過年提出來的，該方法的主要思想是通過將多維數(shù)據(jù)投影到某個(gè)方向上，投影的原則是將總體與總體將多維數(shù)據(jù)投影到某個(gè)方向上，投影的原則是將總體與總體之間盡可能的放開，然后再選擇合適的判別規(guī)則，將新的樣之間盡可能的放開，然后再選擇合適的判別規(guī)則，將新的樣品進(jìn)行分類判別。品進(jìn)行分類判別。一、一、Fisher判別的基本思想判別的基本思想n 二、二、Fisher判別函數(shù)的構(gòu)造判別函數(shù)的構(gòu)造1、針對(duì)兩個(gè)總體的情形、針對(duì)兩個(gè)總體的情形n 2、針對(duì)多個(gè)總體的情形、針對(duì)多個(gè)總體的情形n n 三、線性判別函數(shù)的求法三、線性判別函數(shù)的求法n n n n這里值得注意的是，本書有幾處利用極值原理求極值時(shí)，只這里值得注意的是，本書有幾處利用極值原理求極值時(shí)，只給出了不要條件的數(shù)學(xué)推導(dǎo)，而有關(guān)充分條件的論證省略了，給出了不要條件的數(shù)學(xué)推導(dǎo)，而有關(guān)充分條件的論證省略了，因?yàn)樵趯?shí)際問題中，往往根據(jù)問題本身的性質(zhì)就能肯定有最因?yàn)樵趯?shí)際問題中，往往根據(jù)問題本身的性質(zhì)就能肯定有最大值（或最小值），如果所求的駐點(diǎn)只有一個(gè)，這時(shí)就不需大值（或最小值），如果所求的駐點(diǎn)只有一個(gè)，這時(shí)就不需要根據(jù)極值存在的充分條件判定它是極大還是極小而就能肯要根據(jù)極值存在的充分條件判定它是極大還是極小而就能肯定這唯一的駐點(diǎn)就是所求的最大值（或最小值）。為了避免定這唯一的駐點(diǎn)就是所求的最大值（或最小值）。為了避免用較多的數(shù)學(xué)知識(shí)或數(shù)學(xué)上的推導(dǎo)，這里不追求數(shù)學(xué)上的完用較多的數(shù)學(xué)知識(shí)或數(shù)學(xué)上的推導(dǎo)，這里不追求數(shù)學(xué)上的完整性。整性。n 第五節(jié)第五節(jié) 實(shí)例分析與計(jì)算機(jī)實(shí)現(xiàn)實(shí)例分析與計(jì)算機(jī)實(shí)現(xiàn)n這一節(jié)我們利用這一節(jié)我們利用SPSS對(duì)對(duì)Fisher判別法和判別法和Bayes判別法進(jìn)行計(jì)判別法進(jìn)行計(jì)算機(jī)實(shí)現(xiàn)。算機(jī)實(shí)現(xiàn)。n為研究某地區(qū)人口死亡狀況，已按某種方法將為研究某地區(qū)人口死亡狀況，已按某種方法將15個(gè)已知地區(qū)個(gè)已知地區(qū)樣品分為樣品分為3類，指標(biāo)含義及原始數(shù)據(jù)如下。試建立判別函數(shù)，類，指標(biāo)含義及原始數(shù)據(jù)如下。試建立判別函數(shù)，并判定另外并判定另外4個(gè)待判地區(qū)屬于哪類？個(gè)待判地區(qū)屬于哪類？ X1 ： 0歲組死亡概率歲組死亡概率 X 4 ： 55歲組死亡概率歲組死亡概率 X 2 ：1歲組死亡概率歲組死亡概率 X5 ： 80歲組死亡概率歲組死亡概率 X 3 ： 10歲組死亡概率歲組死亡概率 X6 ： 平均預(yù)期壽命平均預(yù)期壽命 表表4.1 各地區(qū)死亡概率表各地區(qū)死亡概率表 (一一) 操作步驟操作步驟1. 在在SPSS窗口中選擇窗口中選擇AnalyzeClassifyDiscriminate，調(diào)，調(diào)出判別分析主界面，將左邊的變量列表中的出判別分析主界面，將左邊的變量列表中的“group”變量選變量選入分組變量中，將入分組變量中，將變量選入自變量中，并選擇變量選入自變量中，并選擇Enter independents together單選按鈕，即使用所有自變量進(jìn)行判單選按鈕，即使用所有自變量進(jìn)行判別分析。別分析。圖圖4.2 判別分析主界面判別分析主界面2. 點(diǎn)擊點(diǎn)擊Define Range按鈕，定義分組變量的取值范圍。本例按鈕，定義分組變量的取值范圍。本例中分類變量的范圍為中分類變量的范圍為1到到3，所以在最小值和最大值中分別輸，所以在最小值和最大值中分別輸入入1和和3。單擊。單擊Continue按鈕，返回主界面。按鈕，返回主界面。3. 單擊單擊Statistics按鈕，指定輸出的描述統(tǒng)計(jì)量和判別函數(shù)按鈕，指定輸出的描述統(tǒng)計(jì)量和判別函數(shù)系數(shù)。選中系數(shù)。選中Function Coefficients欄中的欄中的Fishers和和Unstandardized。這兩個(gè)選項(xiàng)的含義如下：。這兩個(gè)選項(xiàng)的含義如下：Fishers：給出：給出Bayes判別函數(shù)的系數(shù)。（注意：這個(gè)選項(xiàng)不是判別函數(shù)的系數(shù)。（注意：這個(gè)選項(xiàng)不是要給出要給出Fisher判別函數(shù)的系數(shù)。這個(gè)復(fù)選框的名字之所以為判別函數(shù)的系數(shù)。這個(gè)復(fù)選框的名字之所以為Fishers，是因?yàn)榘磁袆e函數(shù)值最大的一組進(jìn)行歸類這種思想，是因?yàn)榘磁袆e函數(shù)值最大的一組進(jìn)行歸類這種思想是由是由Fisher提出來的。這里極易混淆，請(qǐng)讀者注意辨別。）提出來的。這里極易混淆，請(qǐng)讀者注意辨別。）Unstandardized：給出未標(biāo)準(zhǔn)化的：給出未標(biāo)準(zhǔn)化的Fisher判別函數(shù)（即典型判判別函數(shù)（即典型判別函數(shù)）的系數(shù)（別函數(shù)）的系數(shù)（SPSS默認(rèn)給出標(biāo)準(zhǔn)化的默認(rèn)給出標(biāo)準(zhǔn)化的Fisher判別函數(shù)系判別函數(shù)系數(shù)）。數(shù)）。n單擊單擊Continue按鈕，返回主界面。按鈕，返回主界面。圖圖4.3 Statistics子對(duì)話框子對(duì)話框4. 單擊單擊Classify按鈕，定義判別分組參數(shù)和選擇輸出結(jié)果。按鈕，定義判別分組參數(shù)和選擇輸出結(jié)果。選擇選擇Display欄中的欄中的Casewise results，輸出一個(gè)判別結(jié)果表，輸出一個(gè)判別結(jié)果表，包括每個(gè)樣品的判別分?jǐn)?shù)、后驗(yàn)概率、實(shí)際組和預(yù)測(cè)組編號(hào)包括每個(gè)樣品的判別分?jǐn)?shù)、后驗(yàn)概率、實(shí)際組和預(yù)測(cè)組編號(hào)等。其余的均保留系統(tǒng)默認(rèn)選項(xiàng)。單擊等。其余的均保留系統(tǒng)默認(rèn)選項(xiàng)。單擊Continue按鈕。按鈕。圖圖4.4 Classify子對(duì)話框子對(duì)話框5. 單擊單擊Save按鈕，指定在數(shù)據(jù)文件中生成代表判別分組結(jié)果按鈕，指定在數(shù)據(jù)文件中生成代表判別分組結(jié)果和判別得分的新變量，生成的新變量的含義分別為：和判別得分的新變量，生成的新變量的含義分別為：Predicted group membership：存放判別樣品所屬組別的值；：存放判別樣品所屬組別的值； Discriminant scores：存放：存放Fisher判別得分的值，有幾個(gè)典型判別得分的值，有幾個(gè)典型判別函數(shù)就有幾個(gè)判別得分變量；判別函數(shù)就有幾個(gè)判別得分變量；Probabilities of group membership：存放樣品屬于各組的：存放樣品屬于各組的Bayes后驗(yàn)概率值。后驗(yàn)概率值。n將對(duì)話框中的三個(gè)復(fù)選框均選中，單擊將對(duì)話框中的三個(gè)復(fù)選框均選中，單擊Continue按鈕返回。按鈕返回。6. 返回判別分析主界面，單擊返回判別分析主界面，單擊OK按鈕，運(yùn)行判別分析過程。按鈕，運(yùn)行判別分析過程。圖圖4.5 Save子對(duì)話框子對(duì)話框（二）（二） 主要運(yùn)行結(jié)果解釋主要運(yùn)行結(jié)果解釋1. Standardized Canonical Discriminant Function Coefficients（給出標(biāo)準(zhǔn)化的典型判別函數(shù)系數(shù)）（給出標(biāo)準(zhǔn)化的典型判別函數(shù)系數(shù)）標(biāo)準(zhǔn)化的典型判別函數(shù)是由標(biāo)準(zhǔn)化的自變量通過標(biāo)準(zhǔn)化的典型判別函數(shù)是由標(biāo)準(zhǔn)化的自變量通過Fisher判別法判別法得到的，所以要得到標(biāo)準(zhǔn)化的典型判別得分，代入該函數(shù)的自得到的，所以要得到標(biāo)準(zhǔn)化的典型判別得分，代入該函數(shù)的自變量必須是經(jīng)過標(biāo)準(zhǔn)化的。變量必須是經(jīng)過標(biāo)準(zhǔn)化的。2. Canonical Discriminant Function Coefficients（給出未標(biāo)（給出未標(biāo)準(zhǔn)化的典型判別函數(shù)系數(shù)）準(zhǔn)化的典型判別函數(shù)系數(shù)）未標(biāo)準(zhǔn)化的典型判別函數(shù)系數(shù)由于可以將實(shí)測(cè)的樣品觀測(cè)值直未標(biāo)準(zhǔn)化的典型判別函數(shù)系數(shù)由于可以將實(shí)測(cè)的樣品觀測(cè)值直接代入求出判別得分，所以該系數(shù)使用起來比標(biāo)準(zhǔn)化的系數(shù)要接代入求出判別得分，所以該系數(shù)使用起來比標(biāo)準(zhǔn)化的系數(shù)要方便一些。見表方便一些。見表4.2（a）。）。由此表可知，兩個(gè)由此表可知，兩個(gè)Fisher判別函數(shù)分別為：判別函數(shù)分別為：實(shí)際上兩個(gè)函數(shù)式計(jì)算的是各觀測(cè)值在各個(gè)維度上的坐標(biāo)，這實(shí)際上兩個(gè)函數(shù)式計(jì)算的是各觀測(cè)值在各個(gè)維度上的坐標(biāo)，這樣就可以通過這兩個(gè)函數(shù)式計(jì)算出各樣品觀測(cè)值的具體空間位樣就可以通過這兩個(gè)函數(shù)式計(jì)算出各樣品觀測(cè)值的具體空間位置。置。 1123456212345674.991.8611.6560.8770.7980.0981.57929.4820.8671.1550.3560.0890.0540.69yXXXXXXyXXXXXX 表表4.2（a） 未標(biāo)準(zhǔn)化的典型判別函數(shù)系數(shù)未標(biāo)準(zhǔn)化的典型判別函數(shù)系數(shù)3. Functions at Group Centroids（給出組重心處的（給出組重心處的Fisher判判別函數(shù)值）別函數(shù)值）如表如表4.2 (b) 所示，實(shí)際上為各類別重心在空間中的坐標(biāo)位置。所示，實(shí)際上為各類別重心在空間中的坐標(biāo)位置。這樣，只要在前面計(jì)算出各觀測(cè)值的具體坐標(biāo)位置后，再計(jì)算這樣，只要在前面計(jì)算出各觀測(cè)值的具體坐標(biāo)位置后，再計(jì)算出它們分別離各重心的距離，就可以得知它們的分類了。出它們分別離各重心的距離，就可以得知它們的分類了。Functions at Group Centroids-2.5941.0139.194-.257-6.600-.756GROUP1.002.003.0012FunctionUnstandardized canonical discriminantfunctions evaluated at group means表表4.2（b） 組重心處的組重心處的Fisher判別函數(shù)值判別函數(shù)值 4. Classification Function Coefficients（給出（給出Bayes判別函數(shù)判別函數(shù)系數(shù)）系數(shù)）如表如表4.3所示，所示，GROUP欄中的每一列表示樣品判入相應(yīng)列的欄中的每一列表示樣品判入相應(yīng)列的Bayes判別函數(shù)系數(shù)。在本例中，各類的判別函數(shù)系數(shù)。在本例中，各類的Bayes判別函數(shù)如下：判別函數(shù)如下：第一組：第一組：第二組：第二組：第三組：第三組： 11234565317.2143.9153.190.153.011.0189.3FXXXXXX 21234566202.2164.7171.2100.062.512.1207.0FXXXXXX 31234564982.9134.9144.585.950.010.5181.7FXXXXXX 將各樣品的自變量值代入上述三個(gè)將各樣品的自變量值代入上述三個(gè)Bayes判別函數(shù)，得到三個(gè)判別函數(shù)，得到三個(gè)函數(shù)值。比較這三個(gè)函數(shù)值，哪個(gè)函數(shù)值比較大就可以判斷該函數(shù)值。比較這三個(gè)函數(shù)值，哪個(gè)函數(shù)值比較大就可以判斷該樣品判入哪一類。例如，將第一個(gè)待判樣品的自變量值分別代樣品判入哪一類。例如，將第一個(gè)待判樣品的自變量值分別代入函數(shù)，得到：入函數(shù)，得到： F1=3793.77， F2=3528.32， F3=3882.48比較三個(gè)值，可以看出最大，據(jù)此得出第一個(gè)待判樣品應(yīng)該屬比較三個(gè)值，可以看出最大，據(jù)此得出第一個(gè)待判樣品應(yīng)該屬于第三組。于第三組。Classification Function Coefficients-143.851-164.691-134.862153.137171.185144.462-90.088-99.976-85.94553.00962.52549.97211.00812.09410.520189.261207.003181.714-5317.234-6202.158-4982.880X 1X 2X 3X 4X 5X 6(Constant)1.002.003.00G RO U PFishers linear discriminant functions表表4.3 Bayes判別法的輸出結(jié)果判別法的輸出結(jié)果5. Casewise Statistics（給出個(gè)案觀察結(jié)果）（給出個(gè)案觀察結(jié)果）在在Casewise Statistics輸出表針對(duì)每個(gè)樣品給出了了大部分的判輸出表針對(duì)每個(gè)樣品給出了了大部分的判別結(jié)果，其中包括：實(shí)際類（別結(jié)果，其中包括：實(shí)際類（Actual Group）、預(yù)測(cè)類）、預(yù)測(cè)類（Predicted Group）、）、Bayes判別法的后驗(yàn)概率、與組重心的判別法的后驗(yàn)概率、與組重心的馬氏距離（馬氏距離（Squared Mahalanobis Distance to Centroid）以及）以及Fisher判別法的每個(gè)典型判別函數(shù)的判別得分（判別法的每個(gè)典型判別函數(shù)的判別得分（Discriminant Scores）。出于排版要求，這里給出結(jié)果表的是經(jīng)過加工的，）。出于排版要求，這里給出結(jié)果表的是經(jīng)過加工的，隱藏了其中的一些項(xiàng)目，如表隱藏了其中的一些項(xiàng)目，如表4.4所示。從表中可以看出四個(gè)所示。從表中可以看出四個(gè)待判樣本依次被判別為第三組、第一組、第二組和第三組。待判樣本依次被判別為第三組、第一組、第二組和第三組。Casewise Statistics111.000.297-2.1771.364111.000.236-2.2701.375111.000.117-2.7411.32311.998.507-3.199.638111.000.418-2.582.366221.000.4699.674.231221.000.8688.332-.613221.0005.98510.128-2.518221.0004.7938.3421.760221.000.1019.491-.145331.000.139-6.687-.394331.000.322-7.163-.685331.0005.365-8.655-1.82333.8793.384-4.766-.60833.995.998-5.727-.270ungrouped31.000361.567-20.714-13.498ungrouped1.998.558-3.319.831ungrouped21.00028.66814.0082.086ungrouped31.0001.982-7.595-1.752CaseNumber12345678910111213141516171819ActualGroupPredictedGroupP(G=g| D=d)SquaredMahalanobisDistance toCentroidHighest GroupFunction1Function2Discriminant Scores表表4.4 個(gè)案觀察結(jié)果表個(gè)案觀察結(jié)果表6. 由于我們?cè)谟捎谖覀冊(cè)赟ave子對(duì)話框中選擇了生成表示判別結(jié)果的新子對(duì)話框中選擇了生成表示判別結(jié)果的新變量，所以在數(shù)據(jù)編輯窗口中，可以觀察到產(chǎn)生的新變量。變量，所以在數(shù)據(jù)編輯窗口中，可以觀察到產(chǎn)生的新變量。其中，變量其中，變量dis-1存放判別樣品所屬組別的值，變量存放判別樣品所屬組別的值，變量dis1-1和和dis2-1分別代表將樣品各變量值代入第一個(gè)和第二個(gè)判別函分別代表將樣品各變量值代入第一個(gè)和第二個(gè)判別函數(shù)所得的判別分?jǐn)?shù)，變量數(shù)所得的判別分?jǐn)?shù)，變量dis1-2、dis2-2和和dis3-2分別代表樣分別代表樣品分別屬于第品分別屬于第1組、第組、第2組和第組和第3組的組的Bayes后驗(yàn)概率值。后驗(yàn)概率值。本章結(jié)束本章結(jié)束