Keyword: ダービン・ワトソン検定, 統計量, 有意確率
概要
本サンプルはダービン・ワトソン検定の統計量の計算を行うFortranによるサンプルプログラムです。 本サンプルは以下に示される10個の残差についてダービン・ワトソン検定を行い、統計量と統計量の有意確率の下限、上限を出力します。
※本サンプルはnAG Fortranライブラリに含まれるルーチン g02fcf() のExampleコードです。本サンプル及びルーチンの詳細情報は g02fcf のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
（本ルーチンの詳細はg02fcf のマニュアルページを参照）1 2 3 4
| このデータをダウンロード |
G02FCF Example Program Data 10 2 :: N, IP 3.735719 0.912755 0.683626 0.416693 1.990200 -0.444816 -1.283088 -3.666035 -0.426357 -1.918697 :: End of RES
- １行目はタイトル行で読み飛ばされます。
- ２行目に残差の数(n)と回帰モデルの独立変数の数(ip)を指定しています。
- ３～４行目には残差(res)を指定しています。
出力結果
（本ルーチンの詳細はg02fcf のマニュアルページを参照）1 2 3 4 5
| この出力例をダウンロード |
G02FCF Example Program Results Durbin-Watson statistic 0.9238 Lower and upper bound 0.0610 0.0060
- ３行目にダービン・ワトソン統計量が出力されています。
- ５行目に統計量の有意確率の下限と上限が出力されています。
ソースコード
（本ルーチンの詳細はg02fcf のマニュアルページを参照）
※本サンプルソースコードは科学技術・統計計算ライブラリである「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
| このソースコードをダウンロード |
PROGRAM g02fcfe
! G02FCF Example Program Text
! Mark 23 Release. nAG Copyright 2011.
! .. Use Statements ..
USE nag_library, ONLY : g02fcf, nag_wp
! .. Implicit None Statement ..
IMPLICIT NONE
! .. Parameters ..
INTEGER, PARAMETER :: nin = 5, nout = 6
! .. Local Scalars ..
REAL (KIND=nag_wp) :: d, pdl, pdu
INTEGER :: i, ifail, ip, n
! .. Local Arrays ..
REAL (KIND=nag_wp), ALLOCATABLE :: res(:), work(:)
! .. Executable Statements ..
WRITE (nout,*) 'G02FCF Example Program Results'
WRITE (nout,*)
! Skip heading in data file
READ (nin,*)
! Read in the problem size
READ (nin,*) n, ip
ALLOCATE (res(n),work(n))
! Read in the data
READ (nin,*) (res(i),i=1,n)
! Calculate the statistic
ifail = 0
CALL g02fcf(n,ip,res,d,pdl,pdu,work,ifail)
! Display the results
WRITE (nout,99999) ' Durbin-Watson statistic ', d
WRITE (nout,*)
WRITE (nout,99998) ' Lower and upper bound ', pdl, pdu
99999 FORMAT (1X,A,F10.4)
99998 FORMAT (1X,A,2F10.4)
END PROGRAM g02fcfe