ウィルコクソンの符号付順位和検定の統計量の分布     Last modified: May 15, 2002

 統計数値表編集委員会: 簡約統計数値表. 日本規格協会,東京,1984. にも,n = 50 までの値は出ていますが,代表的な有意水準に対応する統計量という形です。

 以下のスクリプトでは,あっという間に結果が出るでしょう。

 別の所に示した u_table.awk は,この延長線上にあった。いわば,mp_table.awk が1次元とすると,u_table.awk は2次元的な広がりを持つということ。


# mp_table.awk

BEGIN {
   LIMIT = 50
   denom = table[0] = 1
   sum = 0
   for (n = 1; n <= LIMIT; sum += n++) {
       for (i = sum; i >= 0; i--) {
           table[i+n] += table[i]
       }
       denom /= 2
       lim = (n+1)*n/2+1
       cum2 = 1
       cum1 = 0
       printf "\n\nn = %i\n\n", n
       printf "%4s:%15s  %13s%13s%13s\n",
              "S", "度数", "確率", "累積1", "累積2"
       for (i = 0; i < lim; i++) {
           x = table[i]
           y = x*denom
           cum1 += y
           printf "%4i:%15.15g  %13.10f%13.10f%13.10f\n",
                  i, x, y, cum1, cum2
           cum2 -= y
       }
   }
}

・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI