関連情報

メディアン検定(サイズの異なる2標本)

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

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

概要

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

メディアン検定のデータ 

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

入力データ

(本ルーチンの詳細は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

  • 1行目はタイトル行で読み飛ばされます。
  • 2行目に2つの標本のサイズの合計(n)と1つ目の標本のサイズ(n1)を指定しています。
  • 3〜8行目に一つめの標本と2つのめの標本(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

  • 7行目〜8行目は読み込まれた1つ目の標本のデータがそのまま出力されています。
  • 10行目〜12行目は読み込まれた2つ目の標本のデータがそのまま出力されています。
  • 14行目には1つめの標本の中で中央値より小さい値のデータ数が出力されています。
  • 15行目には2つめの標本の中で中央値より小さい値のデータ数が出力されています。
  • 17行目には有意確率が出力されています。

ソースコード

(本ルーチンの詳細は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


Results matter. Trust NAG.

Privacy Policy | Trademarks