计量经济学：实验一、线性回归模型

实验一多元线性回归

［实验目的］使学生掌握数据的输入，掌握线性回归分析的操作方法，进行相关检验，并

能对结果加以分析和解释，

［实验内容］①数据输入、输出及处理；②建立线性回归模型；③相关检验：④对结果加

以分析。

［实验步骤］数据输入处理一建立模型一假设检验一分析预测

问题：对女性劳动参工率的决定因素进行计量分析。数据dataLcsv

方法一：利用gretl软件，gretl界面为:

«F)WMtnML3

点击后会出现下面窗口:

点击OK,找到数据datal.csv,点右下角Open

点击No,得到:

filelootsQauyiewAddSampleVariableModd

dau4-5.gdt

10»，Variablename<Descriptivelabel

0constauto-generatedcomtant

Undated:Fullrange1•50

m2曰™:”惜0LiB&

这里各变量含义:

const-自动生成的固定项

wlfp-16岁以上女性参工率

yfT5岁以上女性收入中位数

ymT5岁以上男性收入中位数

educ-25岁以上女性初中以上学历占比

ue-失业率

mr-15岁以上女性结婚占比

drT5岁以上女性离婚占比

urb-市区人口占比

whT6岁以上白人女性占比

估计的模型如下：

wljp=/?()+pxyf+yin+P3educ+/34ue+ft5mr+(36dr+(^urh+凤wh+e

并给出自变量yf等分别为：17,30,80,6,55,10,60,65时,wlfp的预测值及95%的

置信区间。

菜单操作：

模型估计：点击Model选OrdinaryLeastSquares,打开对话窗:

选中wlfp作为因变量，yf、ym等作为自变量，再点0K,得:

命令操作：

点击File->scriptfiles->Newscript->gretlscript,

E*leloolsgy”-Add，・y，/-卜Help

Appenddm

日$41"dMA

G3Savedataas

£xportd(*t.

QS«ndTo...

0dauMtCtrl3

▲d*AS@t

dorkfngdirectory...

Scriptfiles4UMrBe...

Sessionftlcs4Practkcfile.

RMab.^esgretlscript

functionfll^tgnuplot“rffX

CTH*XRscript

a4riar，vmt»Ya.-

打开程序输入窗口

。启GG9<W»，>U。入〜*o=0X

输入程序

opendata4-5

#导入数据，也可以用open"J:\ProgramFiles\gretl\data\data4-5.gdt*

model1<-olswlfpconstyfymeducuemrdrurbwh

#建模,const也可以写成0

mode11.show

#查看模型结果，或者在grell:scriptoulpul窗口查看

Run

既可以执行单条命令（选中某条命令，点Run按钮即可），也可以同时执行多条命令，直接

点击Run按钮。得到上面同样的结果。

由于ym和mr不显著，P值很大，因此我们设其前面的参数为0,命令窗口中输入

restrict

b[3]=0

b[6]=0

endrestrict

执行后得到：

?restrict

?b[3]=0

?b[6]-0

?endrestrict

Restrictionset

1：b[ym]-0

2：b[mr]=0

Teststatistic：F(2,41)-0.318535,withp-value-0.728997

Restrictedestimates:

coefficientstd.errort-ratiop-value

const41.83365.475287.6401.54e-09***

yf0.8492640.1581525.3702.99e-06***

ym0.00000D0.000000NANA

educ0.2491520.06909873.6060.0008***

ue-1.677580.276859-6.0593.00e-07***

mr0.0000000.000000NANA

dr0.4341040.2250801.9290.0604*

urb-0.09421720.0293363-3.2120.0025***

Wh-0.09608610.0352037-2.7290.0092***

StandarderrorofCheregression-2.14036

P值为0.728997,因此我们可认为这两个参数为0。这样最终模型为:

卬加=41.834+0.849弁+0.249—1.678+0.434dr—0.094〃法一0.096iM

结果分析：

在其它因素不变的条件下，妇女工资每增长1000美元，会带来妇女的参工率增加0.84%；

男性收入影响不显著；

中学毕业妇女人数每增加魏，妇女参工率增加0.249%；

ue前为符号，表明失业率增加，妇女会放弃找工作，退出劳动力市场，失业率每增加战,

妇女参工率降低1.678%

离婚率增加会带来参工率增加，5%显著性水平上不显著；

urb表明农村人口较高地区妇女参工率高，农村妇女要干很多农活;

白种人口增加1%会使妇女参工率降低0.096%。

预测:在主窗口，点击Data->Addobservations***

P|Addobservations

Numberofobservationstoadd:1

CancelOK

点击OK。

回到主窗口，点击Data->Selectall,再点击Data->Editvalues,输入第51个观测值

以便预测。

点Applyo

打开modcll窗口，点Analysis->Forecasts,

改为0,不画

0K,得

这样,我们可以得到妇女参工率的预测值为58.0486,95%的置信区间信2.747,63.350］。

预测的命令操作：

addobs1

scalarn=$nobs

scalark=$nvars

yf[n]=17

ym[n]=30

educ-nl=8()

uc[n]=6

mr[n]=55

dr[n]=10

urb[n]=60

wh[n]=65

feastfwlfp-out-of-sample

scalarwlfpP=fwlfp[n]

genrtval=critical(t,n-k,0.025)

genrwlfpL=wlfpP-tval*$fccrr

genrwlfpU=wlfpP+tval*$fcerr

printwlfpPwlfpLwlfpU

#预测值、预测下限和上限

ff=$yhat

printff

方法二：利用EViws软件，EViews软件界面如下(菜单操作省略):

点击file-New-Program,Hl现窗口:

在窗口里输入下列程序，再点击run：

createu50

cdK:\

readdatal.csv9

equationeql.Iswlfpcyfymeducucmrdrurbwh

freezeeql.results

freezeeql.waldc(3)=0,c(6)=0

scalarn=eql.©regobs

scalark=eql.©ncoef

pagestruct(end=@last+l)*

yf(n+l)=17

ym(n+l)=30

educ(n+l)=80

ue(n+l)=6

mr(n+l)=55

dr(n+l)=10

urb(n+l)=60

wh(n+1)=65

cql.forecastwlfpfs_se

genrypl=wlfpf(n)-@qtdist(0.975,n_k)*s_se(n)

genrypu=wlfpf(n)+Sqtdist(0.975,n-k)*sse(n)

groupexlwlfpf(n)yplypu

showexl

注：文件路彳仝在桌面时cdMC:\DocumentsandSettings\Administrator\桌面"中需加入双引号!

共产生下列四个窗口，其中包括工作表，里面为所生成的对象:

U)Workfil©:UNTITLED[=M目

[View||Proe][object]|Print^||s^ye|Qegils+/-][show||Fet^i|^StoreJ|Delete|[GerwJ[sampl^]

Range-151—51obsDisplayFilter:*

Sample.151—51obs

yrr>

K3ypl

ypu

wlfpf

\Untitled人ZewPage/

估计结果:

A|口o

OependentVanat>leiWUFR

Method:LeastSquares

3Dat❸04/09/10Time1213

4Sample*1no

5Iriczlt-iclcclobservations:50

8

7Vuri&iUleOoefTicaentSld.errort-StiitisticProt>.

Q

9C44.SO9S789749G449593030.0000

ioYFO987983O4075832424007O0198

1T丫M-O1743^5O306207-O5GS37OO£722

-1二,OOO931GG3OGO478O0039

13UE-161OS85O3136175135SO9OOOOO

14MR-O.OT82-14O.1Z31390451743O.6S38

15L>M0.43/3/10238336*1.6930360.0980

1GURB-0.0926340.033335-2.7783400.0082

17WH-0.087492O.O3984S-21958200.0338

-18

-1»R-wq・iparadO700710Maandapacdactv^rG747400

NCJ/XHjusteelR-HquinrudO737922SOdoporidaritvsar4248784

21S.Eofregression217S104人knikuinfocriterion4553579

22Sumnqjxrudrunid193.9742Suhwxr/criterion4.89ZZ43

23LogliKehhooxd-1O4.839S—statistic18.24S9。

24Ourt>in-Watsonstat1.636G47Rrot>(F-statistic)0.000000

OU

2G■I_I>

系数限制检验结果:

□□Table:UNTITLEDWorkfile:UZTITLED::LlntiTed\

[view]|Proc；[Object|[Pont||Name〔Edit♦/•]|CellFmt][1nsOel)〔Grid+/•[Titie](Comme

A1B1C1D1

1WaldTest♦

2Equation:EQ1

3

4FestStatisticValuedfProbability

C

6F-statistic0318535(2,41)07290

7Chi-square0.63707020.7272

Q

9

10NullHypothesisSummary

*11

12NormalizedRestriction(=0)ValueStd.Err.

Ta

14C(3)-0.1743450.306207

15C(6)-00782140173139

16

17Restrictionsarelinearincoefficients.

.O

1O

.c4(