-
正交多項式迴歸
鎖定
正交多項式迴歸是用正交多項式表安排試驗和迴歸分析處理數據。它與用最小二乘法配製的一般多項式迴歸不同,其迴歸係數的估計是互相獨立的,若統計檢驗某一迴歸係數與零無顯著性差異,只需從迴歸方程中刪去這一項,而無需對其他的迴歸係數重新進行計算。
多項式迴歸雖然是一種有效的統計方法,但這種方法存在着兩個缺點:一是計算量較大,特別是當自變量個數較多,或者自變量冪較高時,計算量迅速增加;二是迴歸係數間存在着相關性,從而剔除一個變量後還必須重新計算求出迴歸係數。
- 中文名
- 正交多項式迴歸
- 外文名
- Orthogonal polynomial regression
- 分 類
- 統計方法
- 學 科
- 數理科學
- 用 途
- 處理數據
目錄
- 1 正交多項式
- 2 數學模型
- 3 迴歸係數的顯著性檢驗
- 4 程序框圖
正交多項式迴歸正交多項式
對於定義在區間上[a,b]的一個函數系
,如果其中任何兩個函數在此區間上的積分為零,而他們之中每個函數自乘的積分不等於零,即
正交多項式迴歸數學模型
設變量y和x的n組觀測數據服從以下k次多項式
若把
看作新的自變量,則(3)式就成為一個k元線性模型,其結構矩陣為
在多項式迴歸中遇到的困難是解正規方程係數矩陣的工作量太大,可以使其對角線上的元素不為零,而其餘元素均為零,從而簡化計算,而且同時消去了係數間的相關性。
對於
我們可以通過選擇係數
使得
在正交多項式迴歸中自變量的選擇是等間隔的,設間隔為h,x0=a, 則
為簡化問題,用
代替x作為自變量。在條件許可時,取自變量x1=1,x2=2,…,xn=n。當x1=1,x2=2,…,xn=n時有
,這時驗證以下多項式是正交的,即
對於正交多項式
有
正交多項式迴歸迴歸係數的顯著性檢驗
正交多項式迴歸方程與迴歸係數的顯著性檢驗可利用正交多項式的性質按表1進行。經檢驗不顯著的高次項可以剔除,將其效應併入殘差平方和,自由度也同時併入,如果對迴歸方程精度不滿意,可以增加高次項,而已經計算出的結果不必重算。
[1]
表1 正交多項式迴歸方差分析表
方差來源 | 平方和 | 自由度 | 均方 | F | 顯著性 |
迴歸 | k | ||||
一次項 | 1 | ||||
二次項 | 1 | ||||
... | ... | ... | ... | ... | ... |
k次項 | 1 | ||||
殘差 | n-k-1 | ||||
總計 | n-1 |
正交多項式迴歸程序框圖
1.數學模型:
2.變量及數組説明:
J-正確讀入數據的控制變量;
N-試驗組數;
M-所取正交多項式項數;
X(I)-存自變量數值;
Y(I)-存因變量數值;
Z(I)-存Y(I)的平方項;
E(I,1)-存在正交多項式一次項
;
E(I,2)-存在正交多項式二次項
;
E(I,3)-存在正交多項式三次項
(其中I=1,…N);
S(J)-結構矩陣逆矩陣元素 J=1,2,3;
B(J)-常數項矩陣B J=1,2,3;
D(J)-迴歸係數 J=0,1,2,3;
Q(J)-偏迴歸平方和 J=0,1,2,3;
S0-剩餘平方和;
S-標準離差;
S1-總平方和;
F(J)-F檢驗值。
- 參考資料
-
- 1. 薛為民.正交多項式迴歸及其應用:化學工業出版社,1989
- 2. 正交多項式 .中國知網[引用日期2017-10-21]