（数量化Ｉ類）演習問題Ｂ ４．９表１２の実行結果例 (p.176-177 表１２)
テレビ視聴率を曜日、時間帯
で、予測したい
<< 数量化Ⅰ類(quant1) k134084 関谷 >> Input ファイル名:qut12.dat
タイトル
サンプル数(n):16 アイテム数(ns):2
アイテム毎のカテゴリ数(kn[i]) kn[1]:2 kn[2]:3
ダミー変数の数(kc += kn[i]):5
方程式の実際の変数の数(nn=kc-ns+1):4
変数添え字への変換テーブルjcz[アイテム][カテゴリ]
jcz[1][1]:1 jcz[1][2]:2 jcz[2][1]:0 jcz[2][2]:3 jcz[2][3]:4
方程式の実際の変数の数の確認(k):4
サンプルデータの読み込み
サンプル番号(i) 目的変数(y[i]) アイテムjのカテゴリ番号jz[i][1] jz[i][2]
i: 1 y[i]: 9.000000 jz[ 1][1]:2 jz[ 1][2]:3
i: 2 y[i]: 10.000000 jz[ 2][1]:2 jz[ 2][2]:3
i: 3 y[i]: 6.000000 jz[ 3][1]:1 jz[ 3][2]:2
i: 4 y[i]: 12.000000 jz[ 4][1]:1 jz[ 4][2]:1
i: 5 y[i]: 11.000000 jz[ 5][1]:2 jz[ 5][2]:3
i: 6 y[i]: 18.000000 jz[ 6][1]:2 jz[ 6][2]:1
i: 7 y[i]: 17.500000 jz[ 7][1]:2 jz[ 7][2]:1
i: 8 y[i]: 12.500000 jz[ 8][1]:1 jz[ 8][2]:2
i: 9 y[i]: 14.000000 jz[ 9][1]:2 jz[ 9][2]:3
i:10 y[i]: 13.000000 jz[10][1]:2 jz[10][2]:3
i:11 y[i]: 20.000000 jz[11][1]:1 jz[11][2]:2
i:12 y[i]: 11.000000 jz[12][1]:2 jz[12][2]:2
i:13 y[i]: 6.000000 jz[13][1]:1 jz[13][2]:1
i:14 y[i]: 15.000000 jz[14][1]:1 jz[14][2]:2
i:15 y[i]: 16.500000 jz[15][1]:2 jz[15][2]:2
i:16 y[i]: 8.000000 jz[16][1]:1 jz[16][2]:1
ダミー変数への代入 x[i][(jcz[j][(jz[i][j])])]=1.0;
アイテムNo 1 1 2 2 目的変数y ダミー変数
カテゴリNo 1 2 2 3
<<< (quant1.c)データ行列 X seted [ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
行/列 1 2 3 4 5 6
1 0.000 1.000 0.000 1.000 9.000 1.000
2 0.000 1.000 0.000 1.000 10.000 1.000
3 1.000 0.000 1.000 0.000 6.000 1.000
4 1.000 0.000 0.000 0.000 12.000 1.000
5 0.000 1.000 0.000 1.000 11.000 1.000
6 0.000 1.000 0.000 0.000 18.000 1.000
7 0.000 1.000 0.000 0.000 17.500 1.000
8 1.000 0.000 1.000 0.000 12.500 1.000
9 0.000 1.000 0.000 1.000 14.000 1.000
10 0.000 1.000 0.000 1.000 13.000 1.000
11 1.000 0.000 1.000 0.000 20.000 1.000
12 0.000 1.000 1.000 0.000 11.000 1.000
13 1.000 0.000 0.000 0.000 6.000 1.000
14 1.000 0.000 1.000 0.000 15.000 1.000
15 0.000 1.000 1.000 0.000 16.500 1.000
16 1.000 0.000 0.000 0.000 8.000 1.000
<<< 積和行列 v(x'x) [ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
行/列 1 2 3 4 5 6
1 7.000 0.000 4.000 0.000 79.500 7.000
2 0.000 9.000 2.000 5.000 120.000 9.000
3 4.000 2.000 6.000 0.000 81.000 6.000
4 0.000 5.000 0.000 5.000 57.000 5.000
5 79.500 120.000 81.000 57.000 2751.750 199.500
6 7.000 9.000 6.000 5.000 199.500 16.000
<<< 掃き出し後の v-1 [ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
行/列 1 2 3 4 5 6
1 0.263 0.105 -0.211 -0.105 10.500 1.000
2 0.105 0.342 -0.184 -0.342 15.000 1.000
3 -0.211 -0.184 0.368 0.184 1.500 0.000
4 -0.105 -0.342 0.184 0.542 -3.600 0.000
5 -10.500 -15.000 -1.500 3.600 200.700 0.000
6 -1.000 -1.000 0.000 0.000 0.000 0.000
<< 標準化前の予測誤差など[ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
No. 観測値 予測値 誤差
1 9.00 11.40 -2.40
2 10.00 11.40 -1.40
3 6.00 12.00 -6.00
4 12.00 10.50 1.50
5 11.00 11.40 -0.40
6 18.00 15.00 3.00
7 17.50 15.00 2.50
8 12.50 12.00 0.50
9 14.00 11.40 2.60
10 13.00 11.40 1.60
11 20.00 12.00 8.00
12 11.00 16.50 -5.50
13 6.00 10.50 -4.50
14 15.00 12.00 3.00
15 16.50 16.50 0.00
16 8.00 10.50 -2.50
[標準カテゴリースコア]
アイテムNo. 1 範囲: 4.50
カテゴリNo. 度数 スコア
1 7 -2.53
2 9 1.97
アイテムNo. 2 範囲: 5.10
カテゴリNo. 度数 スコア
1 5 0.56
2 6 2.06
3 5 -3.04
<< 予測誤差など[ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
No. 観測値 予測値 誤差
1 9.00 11.40 -2.40
2 10.00 11.40 -1.40
3 6.00 12.00 -6.00
4 12.00 10.50 1.50
5 11.00 11.40 -0.40
6 18.00 15.00 3.00
7 17.50 15.00 2.50
8 12.50 12.00 0.50
9 14.00 11.40 2.60
10 13.00 11.40 1.60
11 20.00 12.00 8.00
12 11.00 16.50 -5.50
13 6.00 10.50 -4.50
14 15.00 12.00 3.00
15 16.50 16.50 0.00
16 8.00 10.50 -2.50
重相関係数 Ry12...m: 0.490
<<< 標本行列 X(アイテム得点1,2,..,m,y,1.0) [ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
行/列 1 2 3 4
1 15.000 -3.600 9.000 1.000
2 15.000 -3.600 10.000 1.000
3 10.500 1.500 6.000 1.000
4 10.500 0.000 12.000 1.000
5 15.000 -3.600 11.000 1.000
6 15.000 0.000 18.000 1.000
7 15.000 0.000 17.500 1.000
8 10.500 1.500 12.500 1.000
9 15.000 -3.600 14.000 1.000
10 15.000 -3.600 13.000 1.000
11 10.500 1.500 20.000 1.000
12 15.000 1.500 11.000 1.000
13 10.500 0.000 6.000 1.000
14 10.500 1.500 15.000 1.000
15 15.000 1.500 16.500 1.000
16 10.500 0.000 8.000 1.000
<<< 生の積和行列 rs(x'x) [ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
行/列 1 2 3 4
1 2796.750 -162.000 2634.750 208.500
2 -162.000 78.300 -83.700 -9.000
3 2634.750 -83.700 2751.750 199.500
4 208.500 -9.000 199.500 16.000
<<< 積和行列 s(x'x) [ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
行/列 1 2 3 4
1 79.734 -44.719 35.016 -13.031
2 -44.719 73.238 28.519 0.563
3 35.016 28.519 264.234 -12.469
4 13.031 -0.563 12.469 0.063
<<< 相関行列 R(x'x) [ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
行/列 1 2 3
1 1.000 -0.585 0.241
2 -0.585 1.000 0.205
3 0.241 0.205 1.000
<<< 掃き出し後 R-1 [ p.177-表12＿テレビ視聴率(曜日、時間帯） ] >>>
行/列 1 2 3
1 1.918 1.271 -0.723
2 1.271 1.886 -0.693
3 -0.723 -0.693 1.317
<偏相関係数>
pr(y,1): 0.455
pr(y,2): 0.440