11/28/2006

11.12.2 通化線預測值glmval

利用通化式預測值時,可用指令glmval進行計算,並與特定函數連結,其語法如下:


yfit = glmval(B,X,'link')
[yfit, dlo, dhi] = glmval(B,X,'link', stats, clev)
[yfit, dlo, dhi] = glmval(B, X, 'link', stats, clev, N, offset, 'const')

這是配合通化線適配指令glmfit之輸出參數B及連結函數'link',可以利用glmval指令,輸入預測值X,得到對應之新的觀測值yfit。B及連結函數'link'應與原glmfit指令使用的相同。yfit值為利用連結函數中之反函數依X*B的方式求得。

在輸入參數中,大體上其定義與glmfit指令之參數相同,B為利用glmfit求得,clev則為信任水準,其預設值為0.95。其上下限 [yfit-dlo, yfit+dhi]為信任水準之範圍,對應於特定之X值。此範圍僅適用於適配曲線,不能應用於新的觀察值。 N 為二項式N參數,配合glmfit之二項式分配曲線使用。offset及'const'的定義則與前面同。

就上述之汽車里程測定為例:

w = [2100:200:4300]';
poor = [1 2 0 3 8 8 14 17 19 15 17 21]';
total = [48 42 31 34 31 21 23 23 21 16 17 21]';
[b2,d2,s2] = glmfit([w w.^2],[poor total],'binomial');
b2'
wnew = (3000:100:4000)';
[yfit,dlo,dhi] = glmval(b2,[wnew wnew.^2],'logit',s2,0.95,30);
errorbar(wnew,yfit,dlo,dhi);

ans =
-7.3109 0.0002 0.0000