クラスカル・ウォリス検定

ホーム > 製品 > nAG数値計算ライブラリ > nAG Fortranライブラリ > サンプルソースコード集 > クラスカル・ウォリス検定

Keyword: クラスカル・ウォリス, Kruskal-Wallis Test, 検定, 一元配置分散分析

概要

本サンプルはクラスカル・ウォリスの一元配置分散分析（クラスカル・ウォリス検定／Kruskal-Wallis Test) を行うFortranによるサンプルプログラムです。本サンプルは以下に示される５つの標本を分析対象とし、検定統計量と自由度、有意確率を出力します。

クラスカル・ウォリス検定のデータ　

※本サンプルはnAG Fortranライブラリに含まれるルーチン g08aff() のExampleコードです。本サンプル及びルーチンの詳細情報は g08aff のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで

入力データ

（本ルーチンの詳細はg08aff のマニュアルページを参照）

このデータをダウンロード

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 のマニュアルページを参照）

この出力例をダウンロード

 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ライブラリ」のルーチンを呼び出します。
サンプルのコンパイル及び実行方法

このソースコードをダウンロード

    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

Fortranによるサンプルソースコード : 使用ルーチン名：g08aff

概要

入力データ

出力結果

ソースコード