当统计分析的数据不服从正态分布或分布未知时,可使用非参数统计方法对数据进行分析,本文将介绍几种常用的非参数统计方法。
5 v0 y1 P' q3 w# n, G& [: \1 D: W
1.符号秩和检验
% ?8 | i/ [8 b2 l/ X( A5 t类似于成对数据的t检验。当试验前后数据的差值非正态时,可采用符号检验或者符号秩和检验比较试验前后的数据是否有统计学差异。
! E! ^2 [/ k& R9 m. ^ L2 ] 统计假设为H0:差值的中位数Me=0 vs H1: Me≠0。
! o3 k9 ~) C. w* t% |! \ 计算步骤如下:
' X& _1 x) ^3 r. P
1)求配对数据(xi,yi)的差值di=xi-yi。
) M. @% {+ _* y Q" I; v% \! o" @2)对差值的绝对值|di|由小到大排序,记为Ri。若差值的绝对值相等,取其平均秩。引入虚拟变量Ui,当di>0时,Ui=1;当di≤0时,Ui=0。
\( @$ |8 a) ^" T n5 f
3)计算统计量W+。如果Me>0,则总体X的分布关于正数Me对称。从而有P(X>0)>P(X>Me)=1/2,P(X<0)<1/2,故观察到的取正值的样本个数比较多。因此,当W+较大即W+≥c,c=inf{c*:P(W+≥c)≤α}时,拒绝原假设,认为差值大于0。
( j4 {; y0 G$ M+ V- l! u(1)6 b+ K2 b) e% N1 z W- n
1 a$ t. t/ @- D- t0 A9 ]4 x) B6 I
4)当n>50且H0成立时,W+近似服从正态分布,正态统计量为U。取α=0.05,当|U|>U0.975=1.96时,拒绝原假设,认为差值的中位数不等于0。
& l" s6 o ]/ G0 k: j, D/ S
6 p w5 Z5 P, L$ A(2)
5 P" B/ r* q( _6 L9 C3 I, ]- M2. Wilcoxon秩和检验
& c# [7 K) j3 z 类似于团体t检验,在非参数统计中,考虑用Wilcoxon秩和检验比较两组间的统计学差异。设两总体分别为X、Y。
' a2 h# Z- B& b2 v: J
统计假设为H0:X与Y同分布 vs H1:X与Y分布不同
2 X' o u F! p) O4 Z, i
记yi在合样本(x1,…,xm,y1,…,yn)中的秩为Ri,检验统计量为WY。
* k% i' |: ]- x8 } X: [! J: d$ x(3)
" \! s: u* `! X" E - }1 E( v6 e2 {4 q# T* y& g
当WY过大或过小时,拒绝原假设,认为两总体之间有显著差异。当N=m+n>50且H0成立时,WY近似服从正态分布,正态统计量为U。取α=0.05,当|U|>U0.975=1.96时,拒绝原假设,认为X与Y的分布不同。
! _1 t* h- j9 |0 D6 w1 ~
(4)
M6 p! Q, f0 A* b: x# x) F% `
3. K-W秩和检验
* S4 k& `/ M( r1 M2 Z D H# Q类似于方差分析。当有多组(>两组)按完全随机化设计进行比较时,若不满足方差分析的条件,可用K-W秩和检验,
3 r9 G1 M) o$ p0 y 统计假设为H0:各组的均值相等 vs H1: 各组的均值不全相等。
: d' i) d+ i5 ^ E- Q, f/ K0 ? 假设数据为连续型分布,记Rij是xij在合样本(x11,…,x1n1,…,xk1,…,xknk)中的秩。
1 \ F3 Q( r5 ~) e! |) T" V
(5)
% B* p5 l" Z$ q" t q
(6)
4 d+ S' q* h |7 Q/ \% I+ m(3)、(4)中N=n1+n2+…+nk。而SST=SSB+SSW,由于SST为定值,此时的方法分析不同于一般的方差分析(一般的方差分析法需通过计算SSB和SSW的值计算F统计量,再进行假设检验),这里只需计算SSB的值即可。
! v. l# C5 L/ b4 @" {
(7)
9 h% T% t6 K! F9 W4 y5 ~
3 g# c M, ^+ @Kruskal和Wallis提出统计量H,并证明它近似服从卡方分布。
# T6 L; {7 [0 i Z* @) z' ^3 `" E, d
(8)
) g$ z; g# J; M- L4 b0 [2 C. e) e- x" H
取α=0.05,当H>χ20.95时,拒绝原假设,认为各组均值间有统计学差异。