範例:多項遞歸預測保留分析
使用 DOE 函數與程式求解統計上顯著的遞歸。使用會測試所有係數組合的強迫方法。這是保留給進階統計應用程式處理的困難工作。在 PTC Mathcad 中,您可利用遞歸函數確保會納入每種組合,但同時程式的長度仍維持得相當短。
1. 定義一組實驗數據。
2. 從 Data 萃取係數數目。
此實驗的可能遞歸數目如下:
3. 萃取標題列。
4. 萃取失敗率 (Y) 和位勢係數 (X)。
例如,若是係數 "A"、"B" 及 "C",則項數為:
◦ List_Terms 的第一個引數為向量,開頭為代表遞歸常數的空白字串,之後包含所有的係數名稱。
◦ 第二個引數必須以 0 開頭。程式迴圈會以遞歸方式使用此引數,直到取得想要的結果為止。
◦ 傳回的項數為選取的所有可能係數組合。例如,會省略 CBA 或 ACB,因為用於遞歸中會產生與 ABC 相同的結果。
6. 使用函數
rows 和
vec2str,建立將項數對應至單一字元的程式。傳回的清單會用作為參照,以供稍後擷取項數。
傳回清單的第一欄包含字元,第二欄包含對應的項數。
8. 建立決定
polyfitstat 函數產生的遞歸是否符合您需求的程式,包含顯著性
(P) 及遞歸描述資料的詳細程度
(R2)。
9. 建立格式化遞歸結果的程式。此程式會傳回原始項數名稱以及 R2 與 P,以在找出所有遞歸之後進一步排序。
10. 建立依整體遞歸顯著性 (P) 排序結果的程式。
11. 使用函數 polyfitstat 建立處理遞歸分析的程式。此簡短但功能強大的程式可執行下列工作:
◦ 求解所有預測組合。
◦ 反覆運算所有預測組合,並為每項組合建立遞歸。
◦ 測試遞歸是否符合所提供的需求。
◦ 格式化有效的遞歸。
◦ 排序已格式化的遞歸。
12. 呼叫程式 Reg 處理遞歸分析。
此程式在傳回結果之前可能會花掉大量的時間,特別是遞歸數目增加時 (例如顯著程度增加時)。在此情況下,此程式會傳回符合 α 與 R2 所定義之需求的兩項遞歸。
13. 檢視某些遞歸分析的結果。
14. 定義 f 為標準遞歸方程式。
15. 重新定義 f,以採用 Reg 所傳回的遞歸係數。
16. 繪製實驗數據與遞歸。
參照
Philip Leitch, BSc (Applied Biology/Environmental Science Hons), MBA (Hons)