Keyword: クラスカル・ウォリス, Kruskal-Wallis Test, 検定, 一元配置分散分析
概要
本サンプルはクラスカル・ウォリスの一元配置分散分析（クラスカル・ウォリス検定／Kruskal-Wallis Test) を行うFortranによるサンプルプログラムです。 本サンプルは以下に示される５つの標本を分析対象とし、検定統計量と自由度、有意確率を出力します。
※本サンプルはnAG Fortranライブラリに含まれるルーチン g08aff() のExampleコードです。本サンプル及びルーチンの詳細情報は g08aff のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
（本ルーチンの詳細はg08aff のマニュアルページを参照）1 2 3 4 5 6 7 8 9
| このデータをダウンロード |
G08AFF Example Program Data 5 :: K 5 8 6 8 8 :: L 23.0 27.0 26.0 19.0 30.0 29.0 25.0 33.0 36.0 32.0 28.0 30.0 31.0 38.0 31.0 28.0 35.0 33.0 36.0 30.0 27.0 28.0 22.0 33.0 34.0 34.0 32.0 31.0 33.0 31.0 28.0 30.0 24.0 29.0 30.0 :: End of X
- １行目はタイトル行で読み飛ばされます。
- ２行目は標本の数(k)を指定しています。。
- ３行目に各標本に含まれる観測値の個数(l)を指定しています。
- ４～９行目に標本の観測値(x)を指定しています。
出力結果
（本ルーチンの詳細はg08aff のマニュアルページを参照）1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
| この出力例をダウンロード |
G08AFF Example Program Results
Kruskal-Wallis test
Data values
Group Observations
1 23. 27. 26. 19. 30.
2 29. 25. 33. 36. 32. 28. 30. 31.
3 38. 31. 28. 35. 33. 36.
4 30. 27. 28. 22. 33. 34. 34. 32.
5 31. 33. 31. 28. 30. 24. 29. 30.
Test statistic 10.537
Degrees of freedom 4
Significance 0.032
- ８行目に読み込まれた１つ目の標本の観測値が出力されています。
- ９行目に読み込まれた２つ目の標本の観測値が出力されています。
- １０行目に読み込まれた３つ目の標本の観測値が出力されています。
- １１行目に読み込まれた４つ目の標本の観測値が出力されています。
- １２行目に読み込まれた５つ目の標本の観測値が出力されています。
- １４行目には検定統計量が出力されています。
- １５行目には自由度が出力されています。
- １６行目には有意確率が出力されています。
ソースコード
（本ルーチンの詳細はg08aff のマニュアルページを参照）
※本サンプルソースコードは科学技術・統計計算ライブラリである「nAG Fortranライブラリ」のルーチンを呼び出します。
サンプルのコンパイル及び実行方法
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72
| このソースコードをダウンロード |
PROGRAM g08affe
! G08AFF Example Program Text
! Mark 23 Release. nAG Copyright 2011.
! .. Use Statements ..
USE nag_library, ONLY : g08aff, nag_wp
! .. Implicit None Statement ..
IMPLICIT NONE
! .. Parameters ..
INTEGER, PARAMETER :: nin = 5, nout = 6
! .. Local Scalars ..
REAL (KIND=nag_wp) :: h, p
INTEGER :: i, ifail, k, lx, nhi, nlo
! .. Local Arrays ..
REAL (KIND=nag_wp), ALLOCATABLE :: w(:), x(:)
INTEGER, ALLOCATABLE :: l(:)
! .. Intrinsic Functions ..
INTRINSIC sum
! .. Executable Statements ..
WRITE (nout,*) 'G08AFF Example Program Results'
WRITE (nout,*)
! Skip heading in data file
READ (nin,*)
! Read in problem size
READ (nin,*) k
ALLOCATE (l(k))
! Read in number of observations in each sample
READ (nin,*) l(1:k)
! Calculate total number of observations
lx = sum(l(1:k))
ALLOCATE (x(lx),w(lx))
! Read in data
READ (nin,*) x(1:lx)
! Display title
WRITE (nout,*) 'Kruskal-Wallis test'
WRITE (nout,*)
! Display input data
WRITE (nout,*) 'Data values'
WRITE (nout,*)
WRITE (nout,*) ' Group Observations'
nlo = 1
DO i = 1, k
nhi = nlo + l(i) - 1
WRITE (nout,99999) i, x(nlo:nhi)
nlo = nlo + l(i)
END DO
! Perform ANOVA
ifail = 0
CALL g08aff(x,lx,l,k,w,h,p,ifail)
! Display results
WRITE (nout,*)
WRITE (nout,99998) 'Test statistic ', h
WRITE (nout,99997) 'Degrees of freedom ', k - 1
WRITE (nout,99998) 'Significance ', p
99999 FORMAT (1X,I5,5X,10F4.0)
99998 FORMAT (1X,A,F9.3)
99997 FORMAT (1X,A,I9)
END PROGRAM g08affe