Wilcoxon検定 のバックアップ差分(No.1) - アールメカブ

アールメカブ


Wilcoxon検定 のバックアップ差分(No.1)


  • 追加された行はこの色です。
  • 削除された行はこの色です。
ウィルコクスンの符合付き順位検定

統計検定量は以下のように求める.
標本番号 1 ( &mathml{x_1};) の対応するふたつの数値の差を求め,その絶対値を &mathml{|d_1|}; とする.

 # ウィルコクスンの符合付き順位和検定 
 tmpC <- c(54,49,73,65,89,85,78,93,60,79,84)
 tmpD <- c(61,47,64,60,77,85,70,95,56,73,79)
 tmpCD <- tmpC -tmpD
 

いま求めた絶対値を元に順番に並び替える.


 rankCD <- rank(abs(tmpCD[tmpCD != 0]))
 length(rankCD)
 tmpFlag <- ifelse( tmpCD[tmpCD != 0] > 0, -1, 1)
 length(tmpCD)
 rankCD <- rankCD * tmpFlag
 sum(rankCD) / sqrt(sum(rankCD^2)) 
 
 pnorm(-1.939179,0,1)
 qnorm(0.01,0,1)