實(shí)驗(yàn)二一元回歸模型

1、實(shí)驗(yàn)二 一元回歸模型【實(shí)驗(yàn)?zāi)康摹空莆找辉€性、非線性回歸模型的建模方法【實(shí)驗(yàn)內(nèi)容】建立我國(guó)稅收預(yù)測(cè)模型【實(shí)驗(yàn)步驟】【例1】建立我國(guó)稅收預(yù)測(cè)模型。表1列出了我國(guó)19851998年間稅收收入Y和國(guó)內(nèi)生產(chǎn)總值（GDP）x的時(shí)間序列數(shù)據(jù)，請(qǐng)利用統(tǒng)計(jì)軟件Eviews建立一元線性回歸模型。表1 我國(guó)稅收與GDP統(tǒng)計(jì)資料年份稅收GDP年份稅收GDP198520418964199232972663819862091102021993425534634198721401196319945127467591988239114928199560385847819892727169091996691067885199

2、0282218548199782347446319912990216181998926379396一、 建立工作文件菜單方式在錄入和分析數(shù)據(jù)之前，應(yīng)先創(chuàng)建一個(gè)工作文件（Workfile）。啟動(dòng)Eviews軟件之后，在主菜單上依次點(diǎn)擊FileNewWorkfile（菜單選擇方式如圖1所示），將彈出一個(gè)對(duì)話框（如圖2所示）。用戶可以選擇數(shù)據(jù)的時(shí)間頻率（Frequency）、起始期和終止期。圖1 Eviews菜單方式創(chuàng)建工作文件示意圖圖2 工作文件定義對(duì)話框本例中選擇時(shí)間頻率為Annual（年度數(shù)據(jù)），在起始欄和終止欄分別輸入相應(yīng)的日期85和98。然后點(diǎn)擊OK，在Eviews軟件的主顯示窗口將顯示相

3、應(yīng)的工作文件窗口（如圖3所示）。圖3 Eviews工作文件窗口一個(gè)新建的工作文件窗口內(nèi)只有2個(gè)對(duì)象(Object）,分別為c（系數(shù)向量）和resid（殘差）。它們當(dāng)前的取值分別是0和NA（空值）?？梢酝ㄟ^(guò)鼠標(biāo)左鍵雙擊對(duì)象名打開(kāi)該對(duì)象查看其數(shù)據(jù)，也可以用相同的方法查看工作文件窗口中其它對(duì)象的數(shù)值。命令方式還可以用輸入命令的方式建立工作文件。在Eviews軟件的命令窗口中直接鍵入CREATE命令，其格式為：CREATE 時(shí)間頻率類型 起始期 終止期本例應(yīng)為：CREATE A 85 98二、 輸入數(shù)據(jù)在Eviews軟件的命令窗口中鍵入數(shù)據(jù)輸入/編輯命令：DATA Y X此時(shí)將顯示一個(gè)數(shù)組窗口（如圖4

4、所示），即可以輸入每個(gè)變量的數(shù)值 圖4 Eviews數(shù)組窗口三、圖形分析借助圖形分析可以直觀地觀察經(jīng)濟(jì)變量的變動(dòng)規(guī)律和相關(guān)關(guān)系，以便合理地確定模型的數(shù)學(xué)形式。趨勢(shì)圖分析命令格式：PLOT 變量1 變量2 變量K作用：分析經(jīng)濟(jì)變量的發(fā)展變化趨勢(shì)觀察是否存在異常值本例為：PLOT Y X相關(guān)圖分析命令格式：SCAT 變量1 變量2 作用：觀察變量之間的相關(guān)程度觀察變量之間的相關(guān)類型，即為線性相關(guān)還是曲線相關(guān)，曲線相關(guān)時(shí)大致是哪種類型的曲線說(shuō)明：SCAT命令中，第一個(gè)變量為橫軸變量，一般取為解釋變量；第二個(gè)變量為縱軸變量，一般取為被解釋變量SCAT命令每次只能顯示兩個(gè)變量之間的相關(guān)圖，若模型中含有多

5、個(gè)解釋變量，可以逐個(gè)進(jìn)行分析通過(guò)改變圖形的類型，可以將趨勢(shì)圖轉(zhuǎn)變?yōu)橄嚓P(guān)圖本例為：SCAT Y X圖5 稅收與GDP趨勢(shì)圖圖5、圖6分別是我國(guó)稅收與GDP時(shí)間序列趨勢(shì)圖和相關(guān)圖分析結(jié)果。兩變量趨勢(shì)圖分析結(jié)果顯示，我國(guó)稅收收入與GDP二者存在差距逐漸增大的增長(zhǎng)趨勢(shì)。相關(guān)圖分析顯示，我國(guó)稅收收入增長(zhǎng)與GDP密切相關(guān)，二者為非線性的曲線相關(guān)關(guān)系。圖6 稅收與GDP相關(guān)圖三、 估計(jì)線性回歸模型在數(shù)組窗口中點(diǎn)擊ProcMake Equation，如果不需要重新確定方程中的變量或調(diào)整樣本區(qū)間，可以直接點(diǎn)擊OK進(jìn)行估計(jì)。也可以在Eviews主窗口中點(diǎn)擊QuickEstimate Equation，在彈出的方程

6、設(shè)定框（圖7）內(nèi)輸入模型：Y C X 或 圖7 方程設(shè)定對(duì)話框還可以通過(guò)在Eviews命令窗口中鍵入LS命令來(lái)估計(jì)模型，其命令格式為：LS 被解釋變量 C 解釋變量系統(tǒng)將彈出一個(gè)窗口來(lái)顯示有關(guān)估計(jì)結(jié)果（如圖8所示）。因此，我國(guó)稅收模型的估計(jì)式為：這個(gè)估計(jì)結(jié)果表明，GDP每增長(zhǎng)1億元，我國(guó)稅收收入將增加0.09646億元。圖8 我國(guó)稅收預(yù)測(cè)模型的輸出結(jié)果五、 估計(jì)非線性回歸模型由相關(guān)圖分析可知，變量之間是非線性的曲線相關(guān)關(guān)系。因此，可初步將模型設(shè)定為指數(shù)函數(shù)模型、對(duì)數(shù)模型和二次函數(shù)模型并分別進(jìn)行估計(jì)。在Eviews命令窗口中分別鍵入以下命令命令來(lái)估計(jì)模型：雙對(duì)數(shù)函數(shù)模型：LS log(Y) C

7、log(X)對(duì)數(shù)函數(shù)模型：LS Y C log(X)指數(shù)函數(shù)模型：LS log(Y) C X二次函數(shù)模型：LS Y C X X2還可以采取菜單方式，在上述已經(jīng)估計(jì)過(guò)的線性方程窗口中點(diǎn)擊Estimate項(xiàng)，然后在彈出的方程定義窗口中依次輸入上述模型（方法同線性方程的估計(jì)），其估計(jì)結(jié)果顯示如圖9、圖10、圖11圖、12所示。雙對(duì)數(shù)模型：(3.8305) (21.0487) 對(duì)數(shù)模型：(-8.3066) (9.6999) 指數(shù)模型：(231.7463) (27.2685) 二次函數(shù)模型：(7.4918) (3.3422) (3.4806) 圖9 雙對(duì)數(shù)模型回歸結(jié)果圖10 對(duì)數(shù)模型回歸結(jié)果圖11 指數(shù)模型回歸結(jié)果圖12 二次函數(shù)模型回歸結(jié)果六、 模型比較四個(gè)模型的經(jīng)濟(jì)意義豆比較合理，解釋變量也都通過(guò)了T檢驗(yàn)。但是從模型的擬合優(yōu)度來(lái)看，二次函數(shù)模型的值最大，其次為指數(shù)函數(shù)模型。因此，對(duì)這兩個(gè)模型再做進(jìn)一步比較。在回歸方程（以二次函數(shù)模型為例）窗口中點(diǎn)擊ViewActual,Fitted,Residual Actual,Fitted,Residual Table（如圖13）,可以得到相應(yīng)的殘差分布表。圖13 回歸方程殘差分析菜單上述兩個(gè)回歸模型的殘差