STATA多組計(jì)量比較的非參數(shù)檢驗(yàn)命令與輸出結(jié)果說明

1、本節(jié)STATA?命令摘要ranksum 觀察變量,by分組變量kwallis?觀察變量,by分組變量I秩和檢驗(yàn)Ma nn ， Whit ney and Wilcox on非參數(shù)檢驗(yàn)對于計(jì)量資料不滿足正態(tài)分布要求或方差不齊性，但樣本資料之間是獨(dú)立抽取的，那么可以應(yīng)用秩和檢驗(yàn)方法進(jìn)行比擬兩組資料的中位數(shù)是否有差異。STATA命令為：ranksum?觀察變量,by 分組變量例：研究不同飼料對雌鼠體重增加的關(guān)系。表中用x表示雌鼠體重增加克，用group=1表示高蛋白飼料組以及用group=2表示低蛋白飼料組。x13414610411912416110783113129group1111111111x9

2、7123701181018510713294group112222222無效假設(shè)Ho:兩組增加體重的中位數(shù)相同ran ksum x,? by(group)Test: Equality of media ns Two-Sample Wilcox on Ran k-SumSum of Ran ks: 49.5 group = 2Expected Sum: 70? z-statistic -1.73Prob > |z|? 0.0832?為第二組低飼料組的秩的和； 假設(shè)效假設(shè)成立，那么第二組的秩的和期望值為70;秩和統(tǒng)計(jì)檢驗(yàn)量z :對于無效假設(shè)Ho對應(yīng)的p值。在本例中，雖然第二組的秩和為 49.

3、5而期望值估計(jì)為70,但p值為0.0832，所以根 據(jù)該資料和統(tǒng)計(jì)結(jié)果一般不能認(rèn)為用高蛋白飼料喂養(yǎng)能明顯增加雌鼠的體重。多組資料中位數(shù)比擬完全隨機(jī)化設(shè)計(jì)資料的檢驗(yàn)對于完全隨機(jī)化設(shè)計(jì)資料的比擬，假設(shè)各組資料不全服從正態(tài)分布即：至少有一 組的資料均不服從正態(tài)分布或各組的資料方差不齊性，那么可以用 Kruskal a nd Wallis 方法進(jìn) 行檢驗(yàn)Ho:各組的中位數(shù)相同。STATA命令為：kwallis? 觀察變量,by分組變量例：50只小鼠隨機(jī)分配到5個(gè)不同飼料組，每組10只小鼠。在喂養(yǎng)一定時(shí)間后，測定鼠肝中的鐵的含量mg/g如表所示：試比擬各組鼠肝中鐵的含量是否有顯著性差異。用x?表示鼠肝中

4、鐵的含量以及用group=1，2，3，4，5分別表示對應(yīng)的5個(gè)組。x2.231.142.631.001.352.011.641.131.011.70group1111111111x5.590.966.961.231.612.941.963.681.542.59group2222222222x4.53.9210.338.232.074.96.846.423.726group3333333333x1.351.060.740.961.162.080.690.680.841.34group4444444444x1.41.512.491.741.591.3634.815.215.12group5555555555kwallis? x, by(group)Test: Equality of populatio ns (Kruskal-Wallis Test)group?Obs? _Ran kSum1? I。？ 188.502? io? 280.503? 10? 420.004? 10? 95.005? 10? 291.00chi-squared =? 27.856 with 4 d.f.?probability =? 0.0001?為各組的秩和值；為該統(tǒng)計(jì)量的c2檢驗(yàn)值；為無效假