第一節(jié) 各組樣本數(shù)目相同時(shí)的計(jì)算方法
作者:徐榮祥 出版社:中國(guó)科學(xué)技術(shù)出版社 發(fā)行日期:2009年7月示例381健康人與輕度、中度、重度、特重度燒傷病人的血清肌酸(Cre)含量記錄于表381中,問(wèn)燒傷病人各期間的血清肌酸平均含量與健康人有無(wú)顯著差異?
【解題步驟】
第一:先對(duì)這些均數(shù)之間進(jìn)行變異數(shù)分析,即F值檢驗(yàn)。
1計(jì)算各組的平均數(shù)xi,并記錄于表381下框內(nèi)(以下皆同),序號(hào)為③;先將各組資料進(jìn)行合并,求出各組資料的合并值∑xi,序號(hào)為①;將各組個(gè)數(shù)填于表內(nèi)ki欄中,序號(hào)為②;然后根據(jù)求平均值計(jì)算公式求出各組的平均值xi,即為③。
2計(jì)算總平均數(shù)xi,序號(hào)為⑧:將各組的合并值(∑xi )相加,得∑x=33797,序號(hào)為⑥,將各組的個(gè)體數(shù)相加,得n=40,即序號(hào)為⑦;用⑦去除⑥得五個(gè)組資料的總平均數(shù)x=8345,序號(hào)為⑧。
3計(jì)算各組變量平方和∑xi2④及其總值(∑x2)⑨:∑xi2④為各組資料平方后的總和,將它們相加得∑x2=300396,序號(hào)為⑨。
4計(jì)算各組資料和的平方均值(∑xi)2/ki及其總值[∑(x)2/n]:先將各組資料相加,平方后除各組的個(gè)體數(shù),即為各組資料和的平方的均值(∑xi)2/ki,再將它們相加,得:(∑xi)2/ki=288599,序號(hào)為⑩。
5計(jì)算校正數(shù):
6計(jì)算總的離均差平方和:此值是總變異的一個(gè)指標(biāo)。在本題中,引起變異的可能原因有兩個(gè),一個(gè)是由于各期病情的反應(yīng)不同,使組與組之間產(chǎn)生差異(組間差異),另一個(gè)是因組內(nèi)個(gè)體存在的差異(組內(nèi)差異)。其總變異公式為(383):
7自由度的計(jì)算:總變異的自由度為n-1=40-1=39;
組間變異的自由度為n-1,本題為5組,即5-1=4;
組內(nèi)變異的自由度為n(ki-1),本題為5(8-1)=35。
8求比值F:通過(guò)以上計(jì)算,得出了離均差平方和及相應(yīng)的自由度,從而以計(jì)算均方(即變異數(shù)),然后比較兩個(gè)均方,求其比值F,查F值表(表382)判斷差異是否顯著。
組間和組內(nèi)離均差平方和除以相應(yīng)的自由度,即得“組間均方”和“組內(nèi)均方”。以組內(nèi)均方除組間均方所得之商即為F值。方差分析模式于下(表382):
將表381中的數(shù)據(jù)及校正值代入表382中,得(表383):
F值表上端的數(shù)值是較大均方的自由度,其符號(hào)為n1,左側(cè)的數(shù)值是另一個(gè)均方的自由度,其符號(hào)為n2。由本題變異數(shù)分析上表中較大的均方是251199,其自由度是4,故n1=4,又n2=35。查F值概率(P)005表(382),n1=4的縱行與n2=35的橫行相交處的數(shù)值,即265;查F概率(P)001表的數(shù)字為392。因注表中無(wú)35項(xiàng),用內(nèi)插法計(jì)算而得。表382(續(xù))數(shù)據(jù)相當(dāng)于概率(P)005的F值及數(shù)據(jù)相當(dāng)于概率(P)001F值。本例F值為745,大于392,即F>001,故P<001。
9結(jié)果分析:由于P<001,故可認(rèn)為燒傷病人血清肌酸含量與健康人血清肌酸含量有非常顯著性差異。F值顯著性界限(表384)為:
以上為變異數(shù)的計(jì)算,并且已經(jīng)證實(shí)了各組平均數(shù)之間有顯著性差異。但是各組平均數(shù)之間,究竟哪一對(duì)或哪幾對(duì)之間有顯著性差異,單從以上計(jì)算結(jié)果并看不出來(lái),尚需進(jìn)行下一步計(jì)算。
第二:Q值檢驗(yàn):
1列出每?jī)删鶖?shù)相互比較表:該表的繪制原則是必須將各組均數(shù)按大小順序加以排列(表385)第①欄。然后,求出各平均數(shù)相互間的差值(表384)第②欄中。如,4430由10239-5809而得。第③欄是各平均數(shù)與次小平均數(shù)的差數(shù),如2674由9824-7150而得。其余類(lèi)推。
2計(jì)算標(biāo)準(zhǔn)誤:公式為(387),計(jì)算結(jié)果:
3確定各組均數(shù)相差顯著時(shí)所需Q值:Q值是兩均數(shù)之差的標(biāo)準(zhǔn)誤在p=005點(diǎn)的倍數(shù)。從Q值表上看出(表386),表上端為處理數(shù),表左側(cè)為自由度。使用時(shí)先將所要比較的各個(gè)均數(shù)按數(shù)值大小順序排成表,然后兩組相減得出每?jī)蓚€(gè)均數(shù)之差,查Q值表,找出各均數(shù)相差時(shí)所需要的Q值。相鄰兩個(gè)均數(shù)比較時(shí),用處理數(shù)2的Q值,相鄰一個(gè)均數(shù)比較時(shí),用處理數(shù)3的Q值,以此類(lèi)推。本例組內(nèi)變異的自由度為35,在Q值表中的自由度欄內(nèi)沒(méi)有35,故用內(nèi)插法推算,得出各組內(nèi)變異的自由度為35的數(shù)值如下:
4兩均數(shù)相差的界限:將Q值與標(biāo)準(zhǔn)誤相乘所得之積,即為該兩均數(shù)相差的顯著之界限。如:
處理數(shù)為2時(shí)的Qα:35×sx=288×649=1896
處理數(shù)為3時(shí)的Qα:35×sx=347×649=2252
處理數(shù)為4時(shí)的Qα:35×sx=382×649=2479
處理數(shù)為5時(shí)的Qα:35×sx=407×649=2641
不同處理數(shù)計(jì)算結(jié)果填如下(表386):
將所計(jì)算的兩均數(shù)相差界限顯著界值(Q×sx)填入表388中的②、③、④、⑤欄相應(yīng)數(shù)字的下方。
5結(jié)果判斷與分析:根據(jù)表388中的數(shù)據(jù),判斷分析健康人與燒傷病人血清肌酸含量平均數(shù)有無(wú)顯著差異。由表388②和③欄得知,輕度與中度燒傷病人肌酸含量均數(shù)的差值均小于Q×sx值(13,41<1879;753<1879),即健康者同輕度與中度燒傷組之間無(wú)統(tǒng)計(jì)學(xué)差異。由③欄看出,重度與特重度病人肌酸含量均數(shù)的差值均大于Q×sx值(4015>2479;4430>2641),故認(rèn)為他們之間有顯著性差異。由第③欄得知,中度燒傷病人均數(shù)的差值小于Q×sx值,與健康者無(wú)統(tǒng)計(jì)學(xué)差異;而重度與特重度病人肌酸含量之間均數(shù)差值均大于Q×sx值(2674>2252;3089>2479),故認(rèn)為輕度與中度燒傷病人血清肌酸含量與健康人之間無(wú)統(tǒng)計(jì)學(xué)差異,重度與特重度燒傷病人血清肌酸含量與輕度和中度燒傷病人有顯著差異,而特重度病人肌酸含量也顯著高于重度病人。