響應面分析

響應面分析

在迴歸分析中，觀察值y可以表述為：

其中 是自變量 的函數， 是誤差項。

在響應面分析中，首先要得到迴歸方程 ，然後通過對自變量 的合理取值，求得使 最優的值，這就是響應面分析的目的。

[例13.15] 有一個大麥氮磷肥配比試驗，施氮肥量為每畝尿素0，3，6，9，12，15，18kg 7個水平，施磷肥量為每畝過磷酸鈣0，7，14，21，28，35，42kg 7個水平，共49個處理組合，試驗結果列於表13.66，試作產量對於氮、磷施肥量的響應面分析。

表13.66 大麥氮磷肥配比

對於表13.66的數據可以採用二元二次多項式擬合，那麼產量可表示為：

其中Ni、Pj、ij分別表示N、P施用量和誤差，按此模型的方差分析見表13.67。結果表明b2和b3這兩個偏回歸係數不顯著，應該將模型縮減，逐步去掉不顯著的迴歸係數，得到

的模型為： 。使用該模型分析的結果為表13.68，從

中可以看出b1，b4，b5是顯著的，b2達到顯著，該模型的迴歸變異佔總變異的98%，因此可以較好地説明施用N、P對產量的影響。對此資料作多項式迴歸分析的方法可參見第11章和附錄的SAS程序LT13-15.sas。

表13.67 二元二次多項式迴歸分析的方差分析(全模型)

表13.68 二元二次多項式迴歸的方差分析(縮減模型)

表13.69 二元二次多項式迴歸的迴歸係數及其顯著性測驗(縮減模型)

由表13.69，可以列出產量對N、P施用量的迴歸方程為：

由迴歸方程，可以作出產量對N、P施用量的響應曲面圖。

分別對迴歸方程求對N和P的偏導數，並令偏導數等於0，可以求得極值：

，N=13.87(kg)

，P=21.61(kg)

因而由迴歸方程估計得尿素施用量為13.87kg，過磷酸鈣施用量為21.61kg時產量最高。

響應面分析中通過迴歸方程進行預測時一般不能超過自變量的取值範圍，例如氮肥的取值範圍為0至18kg/畝，而磷肥的取值範圍為0至42kg/畝。推論合理的處理組合時，也應該這樣。