出展:パーソナルソフトウェア技法 Watts S. Humphrey 著 ISBN4-320-02929-1

PSP課題 10A_VB : 重回帰による90%と70%の予測

前提条件

11章までと A9,A10,A11 節と課題 06A_VB

要求

n 個の5組の実数(番号、新規LOC、再利用LOC、変更LOC、時間) をファイルから読込み1列目を除いてリストに格納せよ。
これを時間に対する 3 変数の重回帰パラメータを求めよ。 70%、90% の信頼区間を出力せよ。
ベースは 06A_VB とせよ。

テスト

入力ファイル期待される出力
テスト110A_VB_test_1.xml結果1
テスト210A_VB_test_2.xml結果2

最初の行の3つの数字は、推定用の新規LOCと再利用LOCと変更LOCを表す。

参考

z_k = β0 + w_k β1 + x_k β2 + y_kβ3
としたとき、
β0 * n + β1 * sum_i(w_i) + β2 * sum_i(x_i) + β3 * sum_i(y_i) = sum_i(z_i)
β0 * sum_i(w_i) + β1 * sum_i(w_i^2) + β2 * sum_i(w_i*x_i) + β3 * sum_i(w_i*y_i) = sum_i(w_i*z_i)
β0 * sum_i(x_i) + β1 * sum_i(x_i*w_i) + β2 * sum_i(x_i^2) + β3 * sum_i(x_i*y_i) = sum_i(x_i*z_i)
β0 * sum_i(y_i) + β1 * sum_i(y_i*w_i) + β2 * sum_i(y_i*x_i) + β3 * sum_i(y_i^2) = sum_i(y_i*z_i) を解いて β を求める
δ^2 := 1/(n-4) * sum_i((z_i - β0 - β1*w_i - β2*x_i - β3*y_i)^2)<br/>
range := t(α/2,n-4) * δ * sqrt(1 + 1/n + (w_k-Ave(w))^2/sum_i((w_i-Ave(w))^2)
 + (x_k-Ave(x))^2/sum_i((x_i-Ave(x))^2) + (y_k-Ave(y))^2/sum_i((y_i-Ave(y))^2))

解答例