メディアン検定 - サンプルソースコード Fortran

ホーム > 製品 > nAG数値計算ライブラリ > nAG Fortranライブラリ > サンプルソースコード集 > メディアン検定

Keyword: メディアン, Median Test, 検定

概要

本サンプルはメディアン検定 (Median Test) を行うFortranによるサンプルプログラムです。本サンプルは以下に示されるサイズの異なる２つの標本を分析対象とし、中央値より小さい値のデータの個数と有意確率を出力します。

メディアン検定のデータ　

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

入力データ

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

G08ACF Example Program Data
39  16                              :: N,N1
13.0  6.0 12.0  7.0 12.0  7.0 10.0
 7.0 10.0  7.0 10.0  7.0 10.0  8.0 
 9.0  8.0 17.0  6.0 16.0  8.0 15.0  
 8.0 15.0 10.0 15.0 10.0 14.0 10.0 
14.0 11.0 14.0 11.0 13.0 12.0 13.0 
12.0 13.0 12.0 12.0                 :: End of X

１行目はタイトル行で読み飛ばされます。
２行目に２つの標本のサイズの合計(n)と１つ目の標本のサイズ(n1)を指定しています。
３～８行目に一つめの標本と２つのめの標本(x)を指定しています。

出力結果

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

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

 G08ACF Example Program Results

 Median test

 Data values

     Group 1   13.  6. 12.  7. 12.  7. 10.  7.
               10.  7. 10.  7. 10.  8.  9.  8.

     Group 2   17.  6. 16.  8. 15.  8. 15. 10.
               15. 10. 14. 10. 14. 11. 14. 11.
               13. 12. 13. 12. 13. 12. 12.

     13 scores below median in group 1
      6 scores below median in group 2

      Significance   0.00088

７行目～８行目は読み込まれた１つ目の標本のデータがそのまま出力されています。
１０行目～１２行目は読み込まれた２つ目の標本のデータがそのまま出力されています。
１４行目には１つめの標本の中で中央値より小さい値のデータ数が出力されています。
１５行目には２つめの標本の中で中央値より小さい値のデータ数が出力されています。
１７行目には有意確率が出力されています。

ソースコード

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

※本サンプルソースコードは科学技術・統計計算ライブラリである「nAG Fortranライブラリ」のルーチンを呼び出します。
サンプルのコンパイル及び実行方法

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

    PROGRAM g08acfe

!      G08ACF Example Program Text

!      Mark 23 Release. nAG Copyright 2011.

!      .. Use Statements ..
       USE nag_library, ONLY : g08acf, nag_wp
!      .. Implicit None Statement ..
       IMPLICIT NONE
!      .. Parameters ..
       INTEGER, PARAMETER              :: nin = 5, nout = 6
!      .. Local Scalars ..
       REAL (KIND=nag_wp)              :: p
       INTEGER                         :: i1, i2, ifail, n, n1
!      .. Local Arrays ..
       REAL (KIND=nag_wp), ALLOCATABLE :: w(:), x(:)
!      .. Executable Statements ..
       WRITE (nout,*) 'G08ACF Example Program Results'
       WRITE (nout,*)

!      Skip heading in data file
       READ (nin,*)

!      Read in problem size
       READ (nin,*) n, n1

       ALLOCATE (x(n),w(n))

!      Read in data
       READ (nin,*) x(1:n)

!      Display title
       WRITE (nout,*) 'Median test'
       WRITE (nout,*)

!      Output data
       WRITE (nout,*) 'Data values'
       WRITE (nout,*)
       WRITE (nout,99999) '    Group 1  ', x(1:n1)
       WRITE (nout,*)
       WRITE (nout,99999) '    Group 2  ', x((n1+1):n)

!      Perform median test
       ifail = 0
       CALL g08acf(x,n,n1,w,i1,i2,p,ifail)

!      Display results
       WRITE (nout,*)
       WRITE (nout,99998) i1, ' scores below median in group 1'
       WRITE (nout,99998) i2, ' scores below median in group 2'
       WRITE (nout,*)
       WRITE (nout,99997) '     Significance  ', p

99999  FORMAT (1X,A,8F4.0/(14X,8F4.0))
99998  FORMAT (1X,I6,A)
99997  FORMAT (1X,A,F8.5)
    END PROGRAM g08acfe

メディアン検定（サイズの異なる２標本）

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

概要

入力データ

出力結果

ソースコード