Keyword: ダービン・ワトソン検定, 統計量, 有意確率
概要
本サンプルはダービン・ワトソン検定の統計量の計算を行うC言語によるサンプルプログラムです。 本サンプルは以下に示される10個の残差についてダービン・ワトソン検定を行い、統計量と統計量の有意確率の下限、上限を出力します。
※本サンプルはnAG Cライブラリに含まれる関数 nag_durbin_watson_stat() のExampleコードです。本サンプル及び関数の詳細情報は nag_durbin_watson_stat のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
(本関数の詳細はnag_durbin_watson_stat のマニュアルページを参照)1 2 3 4
このデータをダウンロード |
nag_durbin_watson_stat (g02fcc) Example Program Data 2 3.735719 0.912755 0.683626 0.416693 1.9902 -0.444816 -1.283088 -3.666035 -0.426357 -1.918697
- 1行目はタイトル行で読み飛ばされます。
- 2行目に回帰モデルの独立変数の数(p)を指定しています。
- 3から4行目には残差(res)を指定しています。
出力結果
(本関数の詳細はnag_durbin_watson_stat のマニュアルページを参照)1 2 3 4 5
この出力例をダウンロード |
nag_durbin_watson_stat (g02fcc) Example Program Results Durbin-Watson statistic 0.9238 Lower and upper bound 0.0610 0.0060
- 3行目にダービン・ワトソン統計量が出力されています。
- 5行目に統計量の有意確率の下限と上限が出力されています。
ソースコード
(本関数の詳細はnag_durbin_watson_stat のマニュアルページを参照)
※本サンプルソースコードはnAG数値計算ライブラリ(Windows, Linux, MAC等に対応)の関数を呼び出します。
サンプルのコンパイル及び実行方法
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
このソースコードをダウンロード |
/* nag_durbin_watson_stat (g02fcc) Example Program. * * CLL6I261D/CLL6I261DL Version. * * Copyright 2017 Numerical Algorithms Group. * * Mark 26.1, 2017. */ #include <stdio.h> #include <nag.h> #include <nag_stdlib.h> #include <nagg02.h> int main(void) { /* Scalars */ double d, pdl, pdu; Integer exit_status, i, p, n; NagError fail; /* Arrays */ double *res = 0; INIT_FAIL(fail); exit_status = 0; printf("nag_durbin_watson_stat (g02fcc) Example Program Results\n"); /* Skip heading in data file */ scanf("%*[^\n] "); scanf("%ld%*[^\n] ", &p); n = 10; /* Allocate memory */ if (!(res = nAG_ALLOC(n, double))) { printf("Allocation failure\n"); exit_status = -1; goto END; } for (i = 1; i <= n; ++i) scanf("%lf", &res[i - 1]); scanf("%*[^\n] "); /* nag_durbin_watson_stat (g02fcc). * Computes Durbin-Watson test statistic */ nag_durbin_watson_stat(n, p, res, &d, &pdl, &pdu, &fail); if (fail.code != NE_NOERROR) { printf("Error from nag_durbin_watson_stat (g02fcc).\n%s\n", fail.message); exit_status = 1; goto END; } printf("\n"); printf(" Durbin-Watson statistic %10.4f\n\n", d); printf(" Lower and upper bound %10.4f%10.4f\n", pdl, pdu); END: nAG_FREE(res); return exit_status; }