公司財(cái)務(wù)管理回歸分析(參考版)
2025-01-13 08:07本頁(yè)面
【正文】 2022/2/6 47 逐步回歸： （ 1）先在初始模型中取全部自變量： stepwise(x,y) 得圖 Stepwise Plot 和表 Stepwise Table 圖 Stepwise Plot中四條直線都是虛線，說(shuō)明模型的顯著性不好 從表 Stepwise Table中看出變量 x3和 x4的顯著性最差 . 2022/2/6 48 （ 2）在圖 Stepwise Plot中點(diǎn)擊直線 3和直線 4，移去變量 x3和 x4 移去變量 x3和 x4后模型具有顯著性 . 雖然剩余標(biāo)準(zhǔn)差（ RMSE）沒(méi)有太大的變化，但是統(tǒng)計(jì)量 F的 值明顯增大，因此新的回歸模型更好 . To MATLAB(liti51） 2022/2/6 49 （ 3）對(duì)變量 y和 x x2作線性回歸： X=[ones(13,1) x1 x2]。 y=[ ]39。 x3=[6 15 8 8 6 9 17 22 18 4 23 9 8]39。,beta0) To MATLAB(liti6） 2022/2/6 44 betafit = 即 y= ++ 結(jié)果為 : 返 回 2022/2/6 45 逐 步 回 歸 逐步回歸的命令是： stepwise（ x， y， inmodel， alpha） 運(yùn)行 stepwise命令時(shí)產(chǎn)生三個(gè)圖形窗口： Stepwise Plot，Stepwise Table， Stepwise History. 在 Stepwise Plot窗口，顯示出各項(xiàng)的回歸系數(shù)及其置信區(qū)間 . Stepwise Table 窗口中列出了一個(gè)統(tǒng)計(jì)表，包括回歸系數(shù)及其置信區(qū)間，以及模型的統(tǒng)計(jì)量剩余標(biāo)準(zhǔn)差（ RMSE）、相關(guān)系數(shù)（ Rsquare）、 F值、與 F對(duì)應(yīng)的概率 P. 矩陣的列數(shù)的指標(biāo)，給出初始模型中包括的子集（缺省時(shí)設(shè)定為全部自變量） 顯著性水平（缺省時(shí)為 ） 自變量數(shù)據(jù) , 階矩陣 mn? 因變量數(shù)據(jù)， 階矩陣 1?n2022/2/6 46 例 6 水泥凝固時(shí)放出的熱量 y與水泥中 4種化學(xué)成分 x x x x4 有關(guān)，今測(cè)得一組數(shù)據(jù)如下，試用逐步回歸法確定一個(gè) 線性模 型 . 序號(hào) 1 2 3 4 5 6 7 8 9 10 11 12 13x 1 7 1 11 11 7 11 3 1 2 21 1 11 10x 2 26 29 56 31 52 55 71 31 54 47 40 66 68x 3 6 15 8 8 6 9 17 22 18 4 23 9 8x 4 60 52 20 47 33 22 6 44 22 26 34 12 12y 78 . 5 74 . 3 10 4. 3 87 . 6 95 . 9 10 9. 2 10 2. 7 72 . 5 93 . 1 11 5. 9 83 . 8 11 3. 3 10 9. 4數(shù)據(jù)輸入： x1=[7 1 11 11 7 11 3 1 2 21 1 11 10]39。 betafit = nlinfit(X,y,39。 y=[ ... ... ... ]39。 yy=a*x1+b*x2+c*x3+d*x4+e*x5+f*x6。 x5=X(:,5)。 x3=X(:,3)。 x1=X(:,1)。 e=beta0(5)。 c=beta0(3)。 2022/2/6 42 1． 對(duì)回歸模型建立 M文件 : function yy=model(beta0,X) a=beta0(1)。下表列出了 19521981年的原始數(shù)據(jù) ，試構(gòu)造預(yù)測(cè)模型。r39。k+39。,x39。,beta0)； beta 得結(jié)果： beta = 即得回歸模型為： xey1 0 6 ??To MATLAB(liti41) 題目 2022/2/6 40 預(yù)測(cè)及作圖： [YY,delta]=nlpredci(39。,39。 求回歸系數(shù)： [beta,r ,J]=nlinfit(x39。 y=[ 10 ]。 mn?預(yù)測(cè)和預(yù)測(cè)誤差估計(jì)： [Y， DELTA]=nlpredci（ ’ model’, x， beta， r， J） 求 nlinfit 或 nlintool所得的回歸函數(shù)在 x處的預(yù)測(cè)值 Y及預(yù)測(cè)值的顯著性為 1alpha的置信區(qū)間 Y DELTA. ?2022/2/6 39 例 4 對(duì)第一節(jié)例 2，求解如下： 1 、對(duì)將要擬合的非線性模型 y = a xbe / ，建立 m 文件 v o l u m . m 如下： f u n c t i o n y h a t = v o l u m ( b e t a , x ) yhat = b e t a ( 1 ) * e x p ( b e t a ( 2 ) . / x ) 。[ b , bi nt , r , r i nt , s t a t s ] = r e gr e s s ( y , X ) 。 ( x2 . ^ 2) 39。 x2 39。) 2022/2/6 35 在畫(huà)面左下方的下拉式菜單中選” all”, 則 beta、 rmse和 residuals都傳送到 Matlab工作區(qū)中 . 在左邊圖形下方的方框中輸入 1000，右邊圖形下方的方框中輸入 6。 rstool(x,y,39。 x239。 y=[100 75 80 70 50 65 90 100 110 60]39。) 預(yù)測(cè)及作圖 To MATLAB(liti23) 2022/2/6 33 （二）多元二項(xiàng)式回歸 命令： rstool（ x， y，’ model’, alpha） n?m矩陣 顯著性水平 （缺省時(shí)為 ） n維列向量 由下列 4 個(gè)模型中選擇 1 個(gè)（用字符串輸入，缺省時(shí)為線性模型）： l i n e a r （ 線性 ）：mmxxy ??? ???? ?110 p u r e q u a dr a t i c （ 純二次 ）： ???????njjjjmmxxxy12110???? ? i n t e r a c t i o n （ 交叉 ）： ????????? mkjkjjkmmxxxxy1110???? ? q u a dr a t i c （ 完全二次 ）： ???????? mkjkjjkmmxxxxy,1110???? ? 2022/2/6 34 例 3 設(shè)某商品的需求量與消費(fèi)者的平均收入、商品價(jià)格的統(tǒng)計(jì)數(shù) 據(jù)如下，建立回歸模型，預(yù)測(cè)平均收入為 1000、價(jià)格為 6時(shí) 的商品需求量 . 需求量 100 75 80 70 50 65 90 100 1 1 0 60收入 1000 600 1200 500 300 400 1300 1 1 0 0 1300 300價(jià)格 5 7 6 6 8 7 5 4 3 9選擇純二次模型，即 2222211122110 xxxxy ????? ?????法一 直接用多
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1