本節STATA 命令摘要
ranksum觀察變量,by(分組變量)
kwallis 觀察變量,by(分組變量)
· 秩 和檢驗(Mann,WhitneyandWilcoxon非參數檢驗)
對於計量資料不滿足正態分布要求或方差不齊性,但樣本資料之間是獨立抽取的,則可以應用秩和檢驗方法進行比較兩組資料的中位數是否有差異。STATA命令為:
ranksum 觀察變量,by(分組變量)
例:研究不同飼料對雌鼠體重增加的關係(摘自醫學統計方法,金丕煥主編,p218)。表中用x表示雌鼠體重增加(克),用group=1表示高蛋白飼料組以及用group=2表示低蛋白飼料組。
x
134
146
104
119
124
161
107
83
113
129
group
1
1
1
1
1
1
1
1
1
1
x
97
123
70
118
101
85
107
132
94
group
1
1
2
2
2
2
2
2
2
無效假設Ho:兩組增加體重的中位數相同。
ranksumx,by(group)
Test:Equalityofmedians(Two-SampleWilcoxonRank-Sum)
①
SumofRanks:49.5(group==2)
ExpectedSum:70②
③
z-statistic-1.73
Prob>|z|0.0832④
①為第二組(低飼料組)的秩的和;②若無效假設成立,則第二組的秩的和期望值為70;③秩和統計檢驗量z;④對於無效假設Ho對應的p值。
在本例中,雖然第二組的秩和為49.5而期望值估計為70,但p值為0.0832,所以根據該資料和統計結果一般不能認為 用高蛋白飼料餵養能明顯增加雌鼠的體重。
· 多組資料中位數比較(完全隨機化設計資料的檢驗)
對於完全隨機化設計資料的比較,若各組資料不全服從正態分布(即:至少有一組的資料均不服從正態分布)或各組的資料方差不齊性,則可以用KruskalandWallis方法進行檢驗(Ho:各組的中位數相同)。STATA命令為:
kwallis 觀察變量,by(分組變量)
例:50隻小鼠隨機分配到5個不同飼料組,每組10隻小鼠。在餵養一定時間後,測定鼠肝中的鐵的含量(mg/g)如表所示:試比較各組鼠肝中鐵的含量是否有顯著性差別(摘自醫學統計方法,金丕煥主編,p220)。用x表示鼠肝中鐵的含量以及用group=1,2,3,4,5分別表示對應的5個組。
x:
2.23
1.14
2.63
1
1.35
2.01
1.64
1.13
1.01
1.70
group:
1
1
1
1
1
1
1
1
1
1
x:
5.59
0.96
6.96
1.23
1.61
2.94
1.96
3.68
1.54
2.59
group:
2
2
2
2
2
2
2
2
2
2
x:
4.5
3.92
10.33
8.23
2.07
4.9
6.84
6.42
3.72
6
group:
3
3
3
3
3
3
3
3
3
3
x:
1.35
1.06
0.74
0.96
1.16
2.08
0.69
0.68
0.84
1.34
group:
4
4
4
4
4
4
4
4
4
4
x:
1.4
1.51
2.49
1.74
1.59
1.36
3
4.81
5.21
5.12
group:
5
5
5
5
5
5
5
5
5
5
kwallis x,by(group)
Test:Equalityofpopulations(Kruskal-WallisTest)
①
group_Obs_RankSum
110188.50
210280.50
310420.00
41095.00
510291.00
chi-squared=27.856with4d.f.②
probability=0.0001③
①為各組的秩和值;②為該統計量的c2檢驗值;③為無效假設檢驗所對應的p值。
本例結果表明:5組的中位數有顯著的差異。即:5 個不同飼料組的小鼠肝臟中鐵的含量有顯著差異,說明小鼠肝臟中鐵的含量與餵養的飼料有關。