wilcoxon rank sum test R代码

x<-c( x 1 , x 2 , x 3 , . . . , x m x_1, x_2, x_3, ..., x_m )
y<-c( y 1 y_1 , y 2 , y 3 , . . . , y n y_2, y_3, ... ,y_n )

Δ \Delta 代表y的模型与x的模型的位移差,Y = X + Δ \Delta

H 0 H_0 : Δ = 0 \Delta =0 vs H 1 H_1 : Δ 0 \Delta \neq0
wilcox.test(y, x)

H 0 H_0 : Δ = 0 \Delta =0 vs H 1 H_1 : Δ > 0 \Delta >0
wilcox.test(y, x, alternative = "greater")

H 0 H_0 : Δ = 0 \Delta =0 vs H 1 H_1 : Δ < 0 \Delta <0
wilcox.test(y, x, alternative = "less")

如果要计算精确的wilcoxon rank sum statistics的p-value:
wilcox.test(y, x, exact = TRUE)

注意:
R中wilcox.test计算返回的W其实是Mann-Whitney statistic U,U与W之间的关系如下:
W = U + n ( n + 1 ) / 2 W = U+n(n+1)/2
n是 y i y_i 的个数。

P.S: wilcoxon rank sum test的原理之后再来写,希望我到时候还记得。

发布了22 篇原创文章 · 获赞 3 · 访问量 1万+

猜你喜欢

转载自blog.csdn.net/helloworld987456/article/details/102785504