計量經(jīng)濟學(xué)多元線性回歸模型(老師作業(yè)要求范本)范文

計量經(jīng)濟學(xué)?多元線性回歸模型應(yīng)用作業(yè)

一、概述

在當(dāng)今市場上，一國的原油產(chǎn)量與多個因素存在著緊密的聯(lián)系，例如民用汽車擁有量、

宏觀經(jīng)濟等都是影響一國原油產(chǎn)量的重要因素。本次將以中國1990-2006年原油產(chǎn)量

與國內(nèi)民用汽車擁有量、GDP等因素的數(shù)據(jù)，通過建計量經(jīng)濟模型來分析上述變量

之間的關(guān)系，強調(diào)的重要性，從而促進國內(nèi)原油產(chǎn)業(yè)的發(fā)展。

二、模型構(gòu)建過程

1.變量的定義

解釋變量：X1民用汽車擁有量，X2電力產(chǎn)量，X.國內(nèi)生產(chǎn)總值，X,能源消費總量。

被解釋變量：Y原油產(chǎn)量

建立計展經(jīng)濟模型：解釋原油產(chǎn)品與民用汽車擁有量、電力產(chǎn)量、國內(nèi)生產(chǎn)總值、

以及能源消費總量之間的關(guān)系。

2.模型的數(shù)學(xué)形式

設(shè)定原油產(chǎn)量與五個解釋變量相關(guān)美系模型，樣本回歸模型為：

/\AAAAA

y""X”x,"x「四X",

3.數(shù)據(jù)的收集

該模型的構(gòu)建過程中共有四個變量，分別是中國從1990—2006年民用汽車擁有量、電

力產(chǎn)量、國內(nèi)生產(chǎn)總值以及能源消費總量，因此為時間序列數(shù)據(jù)，最后一個即2006年的數(shù)

據(jù)作為預(yù)測對比數(shù)據(jù)，收集的數(shù)據(jù)如下所示：

年份YXIX2X3X4

199019745.18551.364988.25618667.898703

199120130.048606.114927.66821781.5103783

199220271.384691.745148.28826923.5109170

199320768.033817.585886.12735333.9115993

199420896.304941.957005.01148197.9122737

199521419.64410408000.10860793.7131176

199622544.721100.087691.72871176.6138948

199722906.931219.098606.6578973137798

199822986.251319.38821.7584402.3132214

199922920.171452.948311.7189677.1133831

200023345.0181608.919286.41699214.6138553

200123365.651802.0411270.49109655.2143199

200223872.462053.1711648.61120332.7151797

200324248.6162382.9311960.466135822.8174990

200425103.6942693.7114425.257159878.3203227

200525940.3763159.6615852.452183867.9224682

200626305.6643697.3517463.424210871246270

Y原油產(chǎn)量（萬噸標(biāo)準(zhǔn)煤）

XI民用汽車擁有量（萬輛）

X2電力產(chǎn)量（萬噸標(biāo)準(zhǔn)煤）

X3國內(nèi)生產(chǎn)總值（億元）

X4能源消費總量（萬噸標(biāo)準(zhǔn)煤）

4.用OLS法估計模型

回歸結(jié)果，散點圖分別如下：

DependentVariable:Y

Method:LeastSquares

Date:05/04/09Time:18:45

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

20425.46531.159238.454500.0000

-21872140.487949-4.4824650.0007

-01981180.112342-1.7635190.1032

00822710.00821810,010620.0000

00011450.0057330.1997650.8450

R-squared0993314Meandependentvar22751.18

AdjustedR-squared0991085S.D.dependentvar1998.786

S.E.ofregression188.7229Akaikeinfocriterion13,55836

Sumsquaredresid427395.8Schwarzcriterion13.80343

Loglikelihood-110.2461F-statistic445.6871

Durbin-Watsonstat1951792Prob(F-statistic)0.000000

』「()?、

Yi=20425.46-2872X1981X+0.0823X3+0.(X)llX4

d.f=12,R2=0,9933,

Se=(531.1592)(0.4879)(0.1123)(0.0082)(0.0057)

t=(38.4545)(-4.4825)(-1.7635)(10.0106)(0.1998)

。C

?X1

eX2

*X3

*X4

180002000022000240002600028000

Y

三、模型的檢驗及結(jié)果的解釋、評價

2.擬合優(yōu)度檢驗及統(tǒng)計檢驗

R2=0.9933,可以看到模型的擬合優(yōu)度非常高，說明原油產(chǎn)量與上述四個解釋變量之間

總體線性關(guān)系顯著。

?模型總體性檢驗(F檢驗)：給定顯著水平a=0.05,查自由度為(4,12)的F分布表，得

F(4,I2)=3.26,可見該模型的F值遠大于臨界值，因此該回歸方程很明顯是顯著的。但由

于X1與X2系數(shù)不顯著且符號為負，與經(jīng)濟意義不符，因此我們認為解釋變量之間存

在多重共線性。

?變量的顯著性檢驗(t檢驗)：給定顯著水平a=0.05,查自由度為12的t分布表，得

ta/212=2.179,大于該臨界值的的顯著變量為X3：其余的解釋變量未通過檢驗，說明

這些變量與被解釋變量之間不存在顯著的線性相關(guān)關(guān)系。

3.多重共線性的檢驗

⑴相關(guān)系數(shù)檢驗法

CorrelationMa”僅

YX1X2X3X4

Y1.0000000.9576960.9676140.9852250.933689

X10.9576961.0000000.9895570.9913360.979697

X20.9676140.9895571.0000000.9928200.971923

X30.9852250.9913360.9928201.0000000.969240

X40.9336890.9796970.9719230.9692401.000000

DependentVariable:Y

Method:LeastSquares

Date:05/04/09Time:21:23

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

C19443.20294.334666.058150.0000

X12.0722190.16077712.888770.0000

R-squared0.917182Meandependentvar22751.18

AdjustedR-squared0.911661S.D.dependentvar1998.786

S.E.ofregression594.0769Akaikeinfocriterion15.72203

Sumsquaredresid5293910.Schwarzcriterion15.82005

Loglikelihood-131.6372F-statistic166.1204

Durbin-Watsonstat0.428702Prob(F-statistic)0.000000

DependentVariable:Y

Method:LeastSquares

Date:05/04/09Time:21:26

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

C17906.56349.955151.168180.0000

X20.5106100.03439514.845550.0000

R-squared0.936276Meandependentvar22751.18

AdjustedR-squared0.932028S.D.dependentvar1998.786

S.E.ofregression621.1132Akaikeinfocriterion16.45994

Sumsquaredresid4073385.Schwarzcriterion15.55797

Loglikelihood-129.4095F-statistic220.3902

Durbin-Watsonstat0.883430Prob(F-statistic)0.000000

DependentVariable:Y

Method:LeastSquares

Date:05/04/09Time:21:27

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

C19577.31166.2736117.74150.0000

X30.0346860.00155722.279530.0000

R-squared0.970667Meandependentvar22751.18

AdjustedR-squared0.968712S.D.dependentvar1998.786

S.E.ofregression353.5539Akaikeinfocriterion14.68408

Sumsquaredresid1875005.Schwarzcriterion14.78210

Loglikelihood-122.8147F-statistic496.3772

Durbin-Watsonstat0.742903Prob(F-statistic)0.000000

DependentVariable:Y

Method:LeastSquares

Date:05/04/09Time:21:27

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

C16161.76676.688323.883610.0000

X40.0446820.00442510.098640.0000

R-squared0.871776Meandependentvar22751.18

AdjustedR-squared0.863228S.D.dependentvar1998.786

S.E.ofregression739.2064Akaikeinfocriterion16,15916

Sumsquaredresid8'96391.Schwarzcriterion16.25719

Loglikelihood-135.3529F-statistic101.9826

Durbin-Watsonstat0.255859Prob(F-statistic)0.000000

可以看出，Y與X3擬合優(yōu)度R:最大，因此將這個方程作為基本方程，然后往里加入其他

變量。

2.引入第二個變量

DependentVariable:Y

Method:LeastSquares

Date:05/04A)9Time:22:31

Sample:19902006

Includedobservations:17

VariableCoeiicientStd.Errort-StatisticProb.

C19864.14104.2852190.478900000

X30.0731060.00657511.118750.0000

X1-2.3819730.404103-5.8944690.0000

R-squared0.991575Meandependentvar22751.18

AdjustedR-squared0.990372S.D.dependentvar1998.786

S.E.ofregression196.1269Akaikeinfocriterion13.55419

Sumsquaredresid538520.9Schwarzcriterion13.70122

Loglikelihood-112.2106F-statistic823.8987

Durbin-Watsonstat1.979233Prob(F-statistic)0.000000

如上圖所示，引入變量X1后，X1的系數(shù)通不過顯著性檢驗。

DependentVariable:Y

Method:LeastSquares

Date:05/04/09Time:22:35

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

C20909.10611.544334.190650.0000

X30.0604260.0115535.2301850.0001

X2-0.3886140.173173-2.2440800.0415

R-squared0.978427Meandependentvar22751.18

AdjustedR-squared0.975345S.D.dependentvar1998.786

S.E.ofregression313.8446Akaikeinfocriterion14.49446

Sumsquaredresid1378978.Schwarzcriterion14.64150

Loglikelihood-120.2029F-statistic317.4841

Durbin-Watsonstat1.134309Prob(F-statistic)0.000000

如上圖所示，引入變量X2后，其系數(shù)也通不過顯著性檢驗。

DependentVariable:Y

Method:LeastSquares

Date:05/04/D9Time:22:37

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

C2C956.43649.493332.265810.0000

X30.0466450.0056578.2461420.0000

X4-0.0167720.007689-2.1812940.0467

R-squared0.978108Meandependentvar22751.18

AdjustedR-squared0.974980S.D.dependentvar1998.786

S.E.ofregression316.1604Akaikeinfocriterion14,50916

Sumsquaredresid1399403.Schwarzcriterion14.65620

Loglikelihood-120.3279F-statistic312.7479

Durbin-Watsonstat1.105811Prob(F-statistic)0.000000

引入變量X4后，其系數(shù)同樣通不過顯著性檢驗。

綜.上所述，本次模型只引入變量X,,其最終輸出結(jié)果如下:

DependentVariable:Y

Method:LeastSquares

Date:05/04/09Time:21:27

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

C19577.31166.2736117.74150.0000

X30.0346860.00155722.279530.0000

R-squared0.970667Meandependentvar22751.18

AdjustedR-squared0.968712S.D.dependentvar1998.786

S.E.ofregression353.5539Akaikeinfocriterion14.68408

Sumsquaredresid1875005.Schwarzcriterion14,78210

Loglikelihood-122.8147F-statistic496.3772

Durbin-Watsonstat0.742903Prob(F-statistic)0.000000

模型的最終結(jié)果為

Y=19577.31+0.0347

(117.7415)(22.2795)

R2=0.9707,*=0.9687,F=496.3772,DW=0.7429

五、異方差檢驗（懷特檢驗）

WhiteHeteroskedasticityTest:

F-statistic0.313118Probability0.736153

Obs*R-squared0.727872Probability0.694936

TestEquation:

DependentVariable:RESIDA2

Method:LeastSquares

Date:05/04/09Time:22:53

Sample:19902006

Includedobservations:17

VariableCoefficientStd.Errort-StatisticProb.

C172673.899625.641.7332260.1050

X3-1.6309602.141675-0.7615350.4590

X3A2761E-069.62E-060.7913240.4420

R-squared0.042816Meandependentvar110294.4

AdjustedR-squared-0.093925S.D.dependentvar122707.1

S.E.ofregression128340.4Akaikeinfocriterion26.52155

Sumsquaredresid2.31E+11Schwarzcriterion26,66858

Loglikelihood-222.4331F-statistic0.313118

Durbin-Watsonstat1.106513Prob(F-statistic)0.736153

n*R2=0.7279<⑵=5.991,不存在異方差。

六、自相關(guān)檢驗及修正

Breusch-GodfreySerialCorrelationLMTest:

F-statistic6.034154Probability0.027698

Obs*R-squared5.120287Probability0.023648

TestEquation:

DependentVariable:RESID

Method:LeastSquares

Date:05/04/09Time:23:06

Presamplemissingvaluelaggedresidualssettozero.

VariableCoefficientStd.Errort-StatisticProb.

C60.06185145.93720.4115590.6869

X3-0.0008950.001395-0.6412350.5317

RESID(-1)0.6331800.2577622.4564520.0277

R-squared0.301193Meandependentvar-1.10E-12

AdjustedR-squared0.201364S.D.dependentvar342.3271

S.E.ofregression305.9255Akaikeinfocriterion14,44335

Sumsquaredresid1310266.Schwarzcriterion14.59038

Loglikelihood-119.7684F-statistic3.017077

Durbin-Watsonstat1.658208Prob(F-statistic)0.081377

LM=n*R2=5.1203>(1)=3.841,模型存在一階自相關(guān)。

同理，可通過LM檢驗法檢驗是否存在二階自相關(guān)，具體如下:

Breusch-GodfreySerialCorrelationLMTest:

F-statistic2.802700Probability0.097270

Obs*R-squared5.121728Probability0.077238

TestEquation:

DependentVariable:RESID

Method:LeastSquares

Date:05/04/09Time:23:12

Presamplemissingvaluelaggedresidualssettozero.

VariableCoefficientStd.Errort-StatisticProb.

C59,38060152.40560.3896220.7031

X3-0.0008850.001469-0.6025180.5572

RESID(-1)0.6399420.3170722.0182880.0647

RESID(-2)-0.0126740.319160-0.0397110.9689

R-squared0.301278Meandependentvar-1.10E-12

AdjustedR-squared0.140035S.D.dependentvar342.3271

SE.ofregression317.4547Akaikeinfocriterion14.56087

Sumsquaredresid1310107.Schwarzcriterion14.75692

Loglikelihood-r9.7674F-statistic1.868467

Durbin-Watsonstat1.670097Prob(F-statistic)0.184761

。2

LMBRF⑵7</。,。產(chǎn)991，模型不存在二階自相關(guān),

七、科一奧迭代法修正

Breusch-GodfreySerialCorrelationLMTest:

F-statistic0.243740Probability0.787814

Obs*R-squared0.678973Probability0.712136

TestEquation:

DependentVariable:RESID

Method:LeastSquares

Date:05/0W9Time:23:28

Presamplemissingvaluelaggedresidualssettozero.

VariableCoefficientStd.Errort-StatisticProb.

C