函數 > 實驗設計法 > 遞歸分析 > 範例:擬合的品質
範例:擬合的品質
使用 polyfitpolyfitstat 函數實行線性遞歸及變異數分析,以測試擬合的品質。
1. 定義聚合物加工的實驗數據表。反應溫度 t 及催化劑進給率 fr 會影響聚合物的黏度 vy
按一下以複製此運算式
擬合的品質
2. 呼叫 polyfit 函數建立數據的線性遞歸模型。
按一下以複製此運算式
按一下以複製此運算式
3. 計算每項溫度與進給率預測的黏度。
按一下以複製此運算式
按一下以複製此運算式
按一下以複製此運算式
4. 計算殘差 (已計算之模型值與量測值之間的差分)。
按一下以複製此運算式
5. 對照觀測到的黏度、溫度及進給率,繪製殘差。
按一下以複製此運算式
按一下以複製此運算式
按一下以複製此運算式
殘差圖指出,當黏度及溫度的量值分別增加時,觀測黏度及溫度的變異數也會增加。
6. 呼叫 polyfitstat 計算用於線性模型的各種統計。顯示第 8 列 polyfitstat 所傳回的變數分析(ANOVA) 矩陣。
按一下以複製此運算式
按一下以複製此運算式
在變數分析(ANOVA) 矩陣中,變異數來源會分在遞歸與殘差元件之間。遞歸元件會在每個遞歸係數之間進一步分割。但您無法分辨殘差的缺適性及純誤差,因為實驗結果 vy 沒有複製品。
計算與使用遞歸所用的變數分析(ANOVA) 表。
1. 計算因誤差的平方和 (SSE)。
按一下以複製此運算式
SSE 等於 χ2,適合擬合優點的一般公制。此為計算最小平方解時最小化的數量。此誤差是模型與數據擬合程度的量測。其顯示遞歸未說明的偏差有多少。
2. 定義誤差 df_error 的自由度 (與總自由度 df_total 有關),以及參數 df_param 的自由度。自由度是扣除擬合參數數目的數據長度。
按一下以複製此運算式
按一下以複製此運算式
按一下以複製此運算式
3. 定義與總平方和 (SST) 有關的遞歸平方和 (SSR)。
按一下以複製此運算式
按一下以複製此運算式
4. 定義均方誤差 (MSE) 及遞歸均方 (MSR)。將誤差除以適當的自由度。
按一下以複製此運算式
按一下以複製此運算式
5. 建立變異數分析表,以描繪擬合的特性。
平方和
DF
均方
F 係數
Regression
按一下以複製此運算式
按一下以複製此運算式
按一下以複製此運算式
按一下以複製此運算式
Error
按一下以複製此運算式
按一下以複製此運算式
按一下以複製此運算式
Total
按一下以複製此運算式
按一下以複製此運算式
您可比較上表與 polyfitstat 變數分析(ANOVA) 矩陣。
按一下以複製此運算式
6. 預測模型與數據擬合的程度:
按一下以複製此運算式
按一下以複製此運算式
這表示線性遞歸模型可說明黏度 92.7% 的變化性。
7. 定義假設檢驗的顯著等級,以檢驗模型是否符合數據。
按一下以複製此運算式
8. 計算臨界 F 值。
按一下以複製此運算式
9. 檢驗模型擬合數據的假設。
按一下以複製此運算式
接受假設。使用此線性遞歸模型,可預測聚合物的黏度。
參照
Montgomery, D.C., Design and Analysis of Experiments, 5th ed., John Wiley & Sons, New York, 2001, pp.398
這是否有幫助?