Kruskal wallis檢驗(yàn)
一、案例介紹
在一項(xiàng)動(dòng)物實(shí)驗(yàn)中,研究者欲研究A、B兩種菌對(duì)小鼠巨噬細(xì)胞吞噬功能的激活作用,將59只小鼠隨機(jī)分為三組,其中一組為生理鹽水對(duì)照組,用常規(guī)巨噬細(xì)胞吞噬細(xì)胞吞噬功能的檢測(cè)方法。試比較不同實(shí)驗(yàn)條件下小鼠巨噬細(xì)胞的吞噬率有無(wú)差別?如果有差別具體進(jìn)行分析下。
二、問(wèn)題分析
想要比較不同實(shí)驗(yàn)條件下小鼠巨噬細(xì)胞的吞噬率有無(wú)差別,可以考慮使用差異性檢驗(yàn),由于共有小鼠三組,所以可以使用Kruskal wallis檢驗(yàn)進(jìn)行分析如果有差異使用事后多重比較進(jìn)一步分析。
三、軟件操作及結(jié)果解讀
(一) 數(shù)據(jù)導(dǎo)入
1.數(shù)據(jù)格式
首先將數(shù)據(jù)整理成正確的格式,一般X一列,Y為一例,并且分析的數(shù)據(jù)帶有數(shù)據(jù)標(biāo)簽的,需要另添加一個(gè)表格進(jìn)行說(shuō)明,數(shù)據(jù)格式如下:
2.導(dǎo)入數(shù)據(jù)
將整理好的數(shù)據(jù)上傳至SPSSAU系統(tǒng)內(nèi),如下:
上傳結(jié)果如下:
(二) Kruskal wallis分析
1. 軟件操作
Kruskal wallis檢驗(yàn)屬于非參數(shù)檢驗(yàn),分析路徑為點(diǎn)擊【通用方法】→【非參數(shù)檢驗(yàn)】然后進(jìn)行分析:
2. 結(jié)果解讀
從上表可以得出,三組的吞噬率統(tǒng)計(jì)量H值為32.807,p值小于0.05.所以具有顯著性差異。Kruskal wallis檢驗(yàn)的統(tǒng)計(jì)量計(jì)算過(guò)程如下:
式中,R_i^2為各組的秩和,n_i為各組對(duì)應(yīng)的例數(shù),N=sum
n_{i},本例的數(shù)據(jù)具體計(jì)算過(guò)程如下:
其中樣本共有59個(gè),所以計(jì)算如下:
由于相同秩次較多,所以統(tǒng)計(jì)量需要進(jìn)行校正。
式中,m為結(jié)集的數(shù)量,T_{i}=
t_{i}^3- t_{i}。
既然模型顯著,進(jìn)行進(jìn)一步分析每?jī)山M之間是否存在差異?使用Nemenyi方法進(jìn)行分析,分析結(jié)果如下:
C=1-rac{(2^{3}-2)}{8*(4^{3}-4)}=0.9875
從結(jié)果可以看出,A菌組和B菌組p值為0.955>0.05,所以A菌組和B菌組之間并沒(méi)有顯著性差異,但是A菌組和對(duì)照組的p值遠(yuǎn)小于0.05,所以A菌組和對(duì)照組具有顯著性差異,B菌組和對(duì)照組也同理。
四、結(jié)論
利用Kruskal wallis檢驗(yàn)發(fā)現(xiàn)三組的吞噬率統(tǒng)計(jì)量H值為32.807,p值小于0.05.所以具有顯著性差異。進(jìn)一步分析,使用Nemenyi方法,發(fā)現(xiàn)A菌組和B菌組沒(méi)有顯著性差異,但是A菌組和對(duì)照組、B菌組和對(duì)照組均具有顯著性差異。
五、知識(shí)小貼士
1、關(guān)于多重比較問(wèn)題?
如果進(jìn)行非參數(shù)檢驗(yàn)Kruskal-Wallis時(shí)發(fā)現(xiàn)呈現(xiàn)出顯著性,因而可考慮繼續(xù)對(duì)比兩兩組別之間的差異性,選中“Nemenyi兩兩比較”即可輸出【具體這里的多重比較為Nemenyi方法多重比較】;如果Kruskal-Wallis檢驗(yàn)顯示沒(méi)有差異性,則不需要進(jìn)行兩兩比較
2、Nemenyi法和Dunn t檢驗(yàn)如何進(jìn)行選擇?
在非參數(shù)檢驗(yàn)Kruskal-Wallis后,如果共有K個(gè)組別,需要進(jìn)行事后多重比較,如果有明確的實(shí)驗(yàn)組,即希望使用實(shí)驗(yàn)組與其它組別分別進(jìn)行兩兩比較(比較次數(shù)為K-1),此時(shí)使用Dunn t檢驗(yàn)較為適合,與此同時(shí),也可使用Dunn t檢驗(yàn)校正p?值,其是在Dunn t檢驗(yàn)基礎(chǔ)上使用bonferroni法對(duì)p?值進(jìn)行校正。如果僅僅是兩兩組別之間任意的比較,那么使用Nemenyi法較為適合。