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指令，輸入預測值Ｘ，得到對應之新的觀測值yfit。Ｂ及連結函數'link'應與原glmfit指令使用的相同。yfit值為利用連結函數中之反函數依X*B的方式求得。

在輸入參數中，大體上其定義與glmfit指令之參數相同，Ｂ為利用glmfit求得，clev則為信任水準，其預設值為0.95。其上下限 [yfit-dlo, yfit+dhi]為信任水準之範圍，對應於特定之Ｘ值。此範圍僅適用於適配曲線，不能應用於新的觀察值。 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