sas多元線性回歸

1、應(yīng)用數(shù)理統(tǒng)計(jì)報(bào)告應(yīng)用數(shù)理統(tǒng)計(jì)報(bào)告報(bào)告人：宋玲報(bào)告人：宋玲地地 點(diǎn)：計(jì)算機(jī)院軟工實(shí)訓(xùn)室點(diǎn)：計(jì)算機(jī)院軟工實(shí)訓(xùn)室時(shí)時(shí) 間：間：20132013年年1212月月2525日日主要內(nèi)容主要內(nèi)容報(bào)告題目報(bào)告題目1輸輸出出結(jié)結(jié)果果3編寫程序編寫程序2分析及總結(jié)分析及總結(jié)4報(bào)告題目報(bào)告題目 在林木生物量生產(chǎn)率研究中，為了了解林地施肥量（x1，kg）、灌水量（x2，10）與生物量（Y，kg）的關(guān)系，在同一林區(qū)共進(jìn)行了20次試驗(yàn)，觀察值見下表，試建立Y關(guān)于x1,x2的線性回歸方程。 從數(shù)據(jù)中可知，通過給定數(shù)據(jù)之間的關(guān)系，采用 線性回歸的方法分析過程REG.編寫程序編寫程序DATA ct;INPUT x1 x2 y

2、;XSQ=x1*x2;CARDS;54 29 5061 39 5152 26 5270 48 5463 42 5379 64 6068 45 5965 30 6579 51 6776 44 7071 36 7082 50 7375 39 7492 60 7896 62 8292 61 8091 50 8785 47 84106 72 8890 52 92PROC REG;MODEL y=x1 x2/P CLI;MODEL y=x1 x2 xsq/P CLI;RUN;定義變量數(shù)據(jù)數(shù)據(jù)過程步輸出結(jié)果輸出結(jié)果(model1)圖1圖2輸出結(jié)果（輸出結(jié)果（model2） 圖3圖4結(jié)果分析結(jié)果分析 （1）

3、回歸模型是否顯著，顯著水平是多少？復(fù)相關(guān)系數(shù)是多少？ 答：回歸方程顯著，顯著水平是0.0001。復(fù)相關(guān)系數(shù)是0.9659。（2）回歸系數(shù)的估計(jì)值是多少？顯著性如何？ 答：Intercept -4.94048 0.1711 X1 1.53952 0.0001 X2 -0.94385 0.0001 X1與X2的系數(shù)對(duì)于表達(dá)式極顯著，截距項(xiàng)對(duì)應(yīng)的系數(shù)對(duì)表達(dá)式在 0.01下不顯著。（3）寫出回歸方程的表達(dá)式。 y=1.53952x1-0.94385x2 - 4.94048（4）利用殘差（實(shí)測(cè)值與預(yù)測(cè)值之差）、95%置信取間的上下限討論預(yù)測(cè)預(yù)報(bào)效果及預(yù)報(bào)的穩(wěn)定性。 答：根據(jù)上面結(jié)果可知?dú)埐詈?5%置信區(qū)間的上下限的差異很大，最大的達(dá)到7.4640.最小的達(dá)到0.2868.幅度比較大。所以穩(wěn)定性也很差。結(jié)果分析結(jié)果分析（5）對(duì)本問題再求出Y關(guān)于X1、X2的二次多項(xiàng)式回歸方程，并與線性回歸方程比較，說明優(yōu)缺點(diǎn)。 Intercept -14.55333 0.2210 X1 1.66857 |t|的值小于0.0001，遠(yuǎn)小于0.05，說明截距項(xiàng)（即常數(shù)項(xiàng)Intercept）通過檢驗(yàn)，估計(jì)值為5.82331.對(duì)自變量x1分析同樣可以得知，x1系數(shù)通過檢驗(yàn)，估計(jì)值為-0.00493. 所以回歸