Keyword: 順位, 回帰, regression, rank, right-censored
概要
本サンプルは右打ち切りデータが観測値に含まれる場合の順位を使った回帰（Regression using ranks) を行うサンプルプログラムです。 本サンプルは以下に示される標本、説明変数と打ち切り変数を分析対象とします。このサンプルでは回帰分析を行い、スコア統計量、スコア統計量の共分散行列、パラメータ推定、パラメータ推定の共分散行列、パラメータ推定の標準誤差やZ統計量を算出します。
※本サンプルはnAG Toolbox for MATLAB®が提供する関数 g08rb() のExampleコードです。実行にはMATLAB®本体（他社製品）とnAG Toolbox for MATLAB®が必要です。
本サンプル及び関数の詳細情報は g08rb のマニュアルページをご参照ください。
入力データ
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 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87
nv = [int32(40)];
y = [143;
164;
188;
188;
190;
192;
206;
209;
213;
216;
220;
227;
230;
234;
246;
265;
304;
216;
244;
142;
156;
163;
198;
205;
232;
232;
233;
233;
233;
233;
239;
240;
261;
280;
280;
296;
296;
323;
204;
344];
x = [0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
0;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1;
1];
icen = [int32(0);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;1];
gamma = 1e-05;
nmax = int32(40);
tol = 1e-05;
[parvar, irank, zin, eta, vapvec, parest, ifail] = ...
g08rb(nv, y, x, icen, gamma, nmax, tol)
- nv には標本の観測値の数を指定しています。
- y には標本の観測値を指定しています。
- x には説明変数を指定しています。
- icen には打ち切り変数を指定しています。
- gamma には一般化ロジスティック分布を定義するパラメータ値を指定しています。
- nmax には標本サイズの最大値を指定しています。
- tol には観測値の同順位（タイ）についての許容基準を指定しています。
- 最後に本関数を呼び出す構文を指定しています。
出力結果
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 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136
parvar =
7.6526
0.1307
irank =
2
5
6
7
8
9
12
13
14
15
16
17
18
25
28
30
35
15
27
1
3
4
10
11
19
20
21
22
23
24
26
27
29
31
32
33
34
36
10
36
zin =
-0.9494
-0.8682
-0.8250
-0.8250
-0.7800
-0.7487
-0.6462
-0.6092
-0.5707
-0.5307
-0.4873
-0.4418
-0.3942
0.0234
0.2837
0.5198
1.6126
0.4693
1.1837
-0.9750
-0.9230
-0.8960
-0.7164
-0.6820
-0.3179
-0.3179
-0.1440
-0.1440
-0.1440
-0.1440
0.1003
0.1837
0.3948
0.7460
0.7460
1.1543
1.1543
2.1126
0.2836
3.1126
eta =
0.0506
0.1318
0.1750
0.1750
0.2200
0.2513
0.3538
0.3908
0.4293
0.4693
0.5127
0.5582
0.6058
1.0234
1.2837
1.5198
2.6126
0.4693
1.1837
0.0250
0.0770
0.1040
0.2836
0.3180
0.6821
0.6821
0.8560
0.8560
0.8560
0.8560
1.1003
1.1837
1.3948
1.7460
1.7460
2.1543
2.1543
3.1126
0.2836
3.1126
vapvec =
array elided
parest =
4.5840
0.5990
2.7459
0.3615
1.6571
ifail =
0
- parvar はスコア統計量とパラメータ推定の共分散行列を示しています。スコア統計量は上三角行列、パラメータ推定は下三角行列で示されています。
- irank は観測値の順位を示しています。
- zin は順序統計量の関数 g=−f'⁄f （f：確率密度関数、f'：f の導関数）の期待値を示しています。
- eta は順序統計量の関数 g'（関数 g の導関数）の期待値を示しています。
- vapvec は順序統計量の関数 g の分散共分散行列の上三角行列を示しますが、ここでは出力が省略されています。
- parest は以下を示しています。
- １３０行目はスコア統計量を示しています。
- １３１行目はパラメータ推定を示しています。
- １３２行目はカイ二乗統計量を示しています。
- １３３行目はパラメータ推定の標準誤差を示しています。
- １３４行目は近似Z統計量を示しています。
- ifail は関数がエラーを検知しなければ"0"が出力されます。