检验异方差性与调整异方差性

1、实验四异方差性【实验目的】掌握异方差性的检验及处理方法【实验内容】建立并检验我国制造业利润函数模型【实验步骤】【例1】表1列出了 1998年我国主要制造工业销售收入与销售利润的统计资料，请利用统 计软件Eviews建立我国制造业利润函数模型。表1我国制造工业1998年销售利润与销售收入情况行业名称销售利润销售收入行业名称销售利润销售收入食品加工业187.253180.44医药制造业238.711264.1食品制造业111.421119.88化学纤维制品81.57779.46饮料制造业205.421489.89橡胶制品业77.84692.08烟草加工业183.871328.59塑料制品业144.

2、341345纺织业316.793862.9非金属矿制品339.262866.14服装制品业157.71779.1黑色金属冶炼367.473868.28皮革羽绒制品81.71081.77有色金属冶炼144.291535.16木材加工业35.67443.74金属制品业201.421948.12家具制造业31.06226.78普通机械制造354.692351.68造纸及纸品业134.41124.94专用设备制造238.161714.73印刷业90.12499.83交通运输设备511.944011.53文教体育用品54.4504.44电子机械制造409.833286.15石油加工业194.452363

3、.8电子通讯设备508.154499.19化学原料纸品502.614195.22仪器仪表设备72.46663.68一、检验异方差性1. 图形分析检验观察销售利润（Y）与销售收入（X）的相关图（图1） : SCAT X Y图1我国制造工业销售利润与销售收入相关图从图中可以看出，随着销售收入的增加，销售利润的平均水平不断提高，但离散程度也逐步扩大。这说明变量之间可能存在递增的异方差性。残差分析首先将数据排序（命令格式为：SORT解释变量），然后建立回归方程。在方程窗口中点击Resids按钮就可以得到模型的残差分布图（或建立方程后在 Eviews工作文件窗口中点击resid对象来观察）。图2我国制造

4、业销售利润回归模型残差分布图2显示回归方程的残差分布有明显的扩大趋势，即表明存在异方差性。2. Goldfeld-Quant 检验将样本安解释变量排序（SORT X）并分成两部分（分别有1到10共11个样本合19到28共10个样本）利用样本1建立回归模型1 （回归结果如图3），其残差平方和为 2579.587。SMPL 1 10LS Y C XDopendsnl Variable- YMet hod. Leasta res'Date: 011/15/05 Time 20:47 Ssimplci 1 10Includedl10VariableCoefficiEnt StError Sta

5、tistic： Prob.C XIS 7G4ES n nfi5R!3414 ：02D22 1 0S9727 0 019182 4 477937a 3180.0021R-squaredAdjusted R-squared S.E. of regression Sum ftrw-iid Log likp|ih口口ri Durbin-Wgitsoni stt0.714B14 i i.szyitbI7.956S5 2C79.5C7 -Ji gsms2.26012Mean dapendsnt var 5.ID dependent var Akaiko Info ctrltericn Sc hwnirz

6、匚i il.vriicn F-statist I E Prob(F-statZ7.06400 di /0258.79067& 0G1193 20.051920.002051图3样本1回归结果利用样本2建立回归模型2 （回归结果如图4）,其残差平方和为 63769.67。SMPL 19 28LS Y C X图4样本2回归结果计算 F 统计量：F RSQ/RSS = 63769.67/2579.59=24.72, RSS禾口 RSQ 分别是模型3.44，而1和模型2的残差平方和。0.05 时，查 F 分布表得 Foo5 (10 1 1,10 1 1)F 24.72Fo.o53.44，所以存

7、在异方差性3. White 检验 建立回归模型：LS Y C X ，回归结果如图5。I * EVivLK4«ilI i ».1 lnrriTi.FninrTTTr.enjt Ob jcVi o许 Fte虫i ons 世人 rtdewVi W 1 Fr 口K?la. I OTbJ ><! 1 IL 1Friiit | Nunv | Fe* | E 11 rn m. 1> | F u>ir 口y<bt. | S t:ittab Ik dSb IDe pendent XZ-ariihls.YMethod: Least SquaresDate: 01/

8、'15/05 Tims: 21:05Sdnijjlv. 1 20lnelu£|*r| obRivtiansT 护VariableCoeffici&ntStd. Errort-StstiGtiProb.u12.C334919 610090 51653U115429X0.104394.06421 2.36653 .尺-squared口 .354694Medependent varZ1 3 4539Adjusted R-squared0.349105S.Dl dopondnt *ar146.1905S. IE - 口 f re grBE-sion55.90455Akcti

9、kt? infocri lerion1 0.99939Sum sqiLjareci reaid04191 .34Schwarz criterion1 1 -口B453Loq likelihood151 8513F eterietlc1 52.0322Durbiri-Walson el at2.497410Prob(F- statiffti c).oooooa图5我国制造业销售利润回归模型检验结果如图6。在方程窗口上点击 ViewResidualTestWhite Heteroskedastcity,辛，£丫“" - E讯.札业；。!二 imiTUl f orlfile： i

10、miTLEDjFile Edit Objects Hew Frees Qiuick Ogtioils Window Heljictn | Print | Name |FreezeEstinate|FWFiite Heterokedasticity Test:F-st atistic Ob s*R-squared3.607090 Probability6.270439 Probability0 0420400.043490图6 White检验结果其中F值为辅助回归模型的F统计量值。取显著水平0.05 ,由于20.0525.99 nR 6.2704,所以存在异方差性。 实际应用中可以直接观察相伴概

11、率p值的大小，若p值较小，则认为存在异方差性。反之，则认为不存在异方差性。4. Park检验建立回归模型(结果同图5所示)。生成新变量序列：GENR LNE2=log(RESID"2)GENR LNX=log建立新残差序列对解释变量的回归模型：LS LNE2 C LNX ，回归结果如图7所示。图7 Park检验回归模型从图7所示的回归结果中可以看出，LNX的系数估计值不为 0且能通过显著性检验，即随即误差项的方差与解释变量存在较强的相关关系，即认为存在异方差性。5. Gleiser检验(Gleiser检验与Park检验原理相同)建立回归模型(结果同图5所示)。生成新变量序列：GENR

12、 E=ABS(RESID)分别建立新残差序列(日 对各解释变量(X/XA2/XA(1/2)/XA( 1)/ XA( 2)/ XA( 1/2)的回归模型：LS E C X ，回归结果如图 8、9、10、11、12、13所示。IWVi.b l_<t 1irwTTTr-nm a. If 专 HL" ir-TTTT-tn I1 1 riJLs Mdi L ObjWi wwFt 口亡凰9*LL.iCjD L i ulThfii W L r& J 08WH-1Pr I Fi iTifr |C 1Hr 电Frit 事厂丹=t | 呂申fl. 11 | R* iiDitjfj um d

13、wi il匚Least SoilresEaate: 01/1505 Tiimm: 21 :S1Ssmple: i 2lncliud«d2Q/曰曰匕ImCaefflioienlStd. Error1-St al i sl idPirah.C12 Z393BID &1 BC1. 1 52B12O 26S6Xfl CH 267口 JOlCmS3.321123O 口02曰R- qu3rC!dl0 29242Muh ri dkpehdifrt vnrJIAdjusted R-sciuAred0.271251S.D. d«p0nd«M var35 26573S E o

14、f roq rsssi 口nAk alike info 匚ill 日ri口口S.Z71S47Sum s q ua)re d re std24Q-19. 7BSc lrrwa rz. izrit e nanLeg lik»lilicod-134.0073F-st4ftistieII.OJOCODu rb i n-Wn &|2.320361Prob(尸衍4目1H吕丰ic) .2&4黑上常$ 百剋 二inrTTTI-W!I>V«r IkFi TIJWTTTTJEllIlJ F i 1 星*ii «5J*_j as I kVi * F r q <

15、;r 勻畑i 口女3g t i a n,MW i riIfplp1 mW 1 Ft'XiCS | C&J w GXS ITri rut 1 &lOJTi>I 口丁中 QED |:E匸 El ar口 1 T旷囁匚直Et1 氏切D«p«ndQnl SHriQblLF. EK/lelhod Le ast S q ureDaie! O1 /1Q/n3 TIftie：21占弓S a mpls 1 卫曰Inc Budle d obse rva i i 0risVq ri dbloCcoffiCii 0111Std Errort-StariefiePro b

16、iC27.D5B3?0.23252B3.2BS7G3LJ.OD29x*x2.74E CB1 sQ2g O&2,日日日8吕2sy 123R-Siquia redl .Z21 7S93IIVIcam de pen de rrt "varA .E3ES-4Adjus 1 eihI F-f-E-qtii r>ciiiiO. 1B7B1 11S D d ei puncliarHl var3E. 26573S 冃 口旷 rnij rR-s ninnap mzpsAH講ik曰 infn nrito rinnA HHOfiSiZSum quaired resid5779 SBSchwr

17、 c-riterio n9.27&7=50Lng llkeliliood-1.36.3264F-«t«tl silo7 235303Durtain-WsitsonAtlasersProb(F-<t*tistic)0.012319I H 更i 工匚 Rd.il£14w £vqc e QuiOELiIfiEdww巴仝十 JgE* 爲口 盍jj 占- In> gjt亡| 空出1仝耳点 曰 上三黒严三1日話Dependfifem Variftbl®- EMetliiod: Least Sq u a r&aUatE :口 1

18、/I sms TlrnE NV:土心3 rn pin 1 20Inc-luicJd olbse»rvati d n s7 2SVoriQblE!-coi'irneii&ntstdl. Errort-SisiisiisProbcX( 1Z2) i 匚 G7GA31.3Q&17S17 0aG42-0.91 G9SS . 309207a SB 1弓趾吕0.3G7GCL 002IR- n q lj h= cdAdjueled R-&c|ijciiFed S 巨。riif r<»i©nSu m scuaii'e d resi d

19、 Log IIIk«=!Hui11 口口(1Ourbin-Wai = on si口. 32630 .30200g 2A7S723QS6.67-I £4 NOBN百二m日冲IM eain d e pe-nident var S O dipienriisnl gr Akni k#* infn c lit Arian S c： h *varz ciiterionRrobfF- laiiigii匚 J4LI 153(5&436.2BS73 a zaeTRfl 9.B239251 2 l=B4&U 口 DC I 45 E1图10图11i1塁iHFdi t-Qtj 旦*

20、書Vi "w畧尹-旦倉哺生利-i 百R jp ti ari-4W i n sri>nwW-l-FVi ow j| Fr 口 <?亘 | Dbj ad -a J Pr i !ii t W 口皿由 | Fraaz a |Re I irriQ t 口F心itau*= t 1 Stat = 1 Rusi d.cDnpfinrlHnl VarighlR £ hxlat hoHl_口nut Stu a口 =f 口： S'1 /1 ti/U51 lirn口 2 1 5SS-?tmpl® *1 26Iriujluidlucd nkiserwEai lams Z

21、?BVdrttiblQCMwfTiClwfiiSliJ. irark£1 aiLlGrProb.CX(-2J4G.93S99-2O227fc7. 1 9061116.5 1SEC4179*3010.-rOOliASG0.0000U.UUSR-sqijarecIAdjustod R-qusr-ed S E mf r=igrsi：saio>n Sum 各qum厂日H resid Ld g II Ik 目 liheczfdDi 1 rh i n-W =»1 = 口 n e! =ii0.1 1 09790 O7G7063-4 Q4SS9 31569.5B m 130. 110

22、71 SH1 了口卩Mi&a n d&'p e n dent vaiS D depeitisdent var Ak Aiik a info sritrion Sc bwsrz c-rit e ricnF- =i sa t is k ic|=口 h(F-st=>frisf ic)Ml .69654 3S 2G5731 OOB4B10i 103643.24SB41 口，口 曰 MANIVHT I T l-JKII Ik 一UHTJ-TIUSUMi.1 丄 . UMT丄ITUS川!11 Eile Kdi t Qbeet 3 ZiErocsSni ck Ob； t i 1

23、1*31 jtiinVkg rr qcu| Old1xu |u | Fvrw up.II Ru 匚 a 己匚Dependent Varisblej EMnl hindSqUrmDale: O1 /4S/DS Time: 2d EGSarYiple: 1 2CInclude J obGorvaliiont-Vania b le亡口 effi cie ntStd. ErTorT-Statistic尸 nzihC x(-1 Z2)0S QO7eS-1 S1 1 1 1 1IS D9I3345 Z5229749e 1931 -3 2490,4O.ODCD 口 口 IZT3 卫W-s-qijaredAd

24、jiusld R-aqMared S E of regreBsion Sum s quare d resi dl Lu gi li k EfiiihdudDurbi n-WnLi.2btJa7 i I.2SM32 31 173C9S I135 1HJOB2.320961Meein clE|n曰ndeint vsir S D de pendent w Akaike info c riteii a n BcHwsirz. c：iriierionF- Hl Lilli 9 IlL!Prob(F-*tirtlfftiic)41 _6«&643B.26B739.7057379.B3O0

25、ZS4ID. 42390.003206图13由上述各回归结果可知，各回归模型中解释变量的系数估计值显著不为0且均能通过显著性检验。所以认为存在异方差性。2由F值或R确定异方差类型2Gleiser检验中可以通过F值或R2值确定异方差的具体形式。本例中，图10所示的回2归方程F值（R ）最大，可以据次来确定异方差的形式。二、调整异方差性1.确定权数变量根据Park检验生成权数变量：GENR W1=1/XX.6743根据Gleiser 检验生成权数变量：GENR W2=1/XA0.5另外生成：GENR W3=1/ABS(RESIDGENR W4=1/ RESID 人22.利用加权最小二乘法估计模型在

26、Eviews命令窗口中依次键入命令:LS(W= ) Y C X或在方程窗口中点击 Estimate'Option 按钮，并在权数变量栏里依次输入W1 W2 W3W4回归结果图14、15、16、17所示。Equat i onZ IJMTTTLEB Voi-kf ile 2 inTTTLEDDependant Vanable: YMethod: Least SquaresDale： 01/154)5 Time: 22:41Sample: 1 26Included observations: 38Wsightirig series: 1/XA1.b734VariableCoefificien

27、tStd. Errort-St artisticProb.c5 9206553.8838581.5244260.1395X0 1085870.00994210 922230.0000Equation： UMTITLEB torkfile： 1MTITLEDniMiw山如和.B!倔丄¥i nr I From 口 bj.util Fr int | in. | Fth | Ext i凯t.|Fr 雹xt |st龟Ji.Dapendent Variable: YMethod: Ljeast SquaresDate: 01/1£«5 Tima: 22:44Sample: 1

28、 28Included observations; 2QWeightirig series:z UHTIT UDoxkfile： UNTITLEDVariableCoefficientStd. Errort-StatisticProb.C8 63927111.10760J. 772213D.447Q冥0.1061530.00774613704300.0000图15M Elation： UariTLEDtorkfilez ITHTITLEDVi ew |Frocs Objects | Print Hams |EsForecast IstatsDep endent Variable: Y Meth

29、od： Least SquaresDate; 01/15AJ5 Time； 22:49 Sample: 1 28Included observations: 20Weighting series. VASS(RESID)VariableCoefficientStd. Error卜 StatisticPrabC4.1689333.7797651.102964D2801X0 109408C.00363330.966530.0000图16Vi 舟|Ft1 ObjMtsI Frint| Ha亡| frw匕工EstimallForecBilstrtmlKEsi.d_sDependent Variable: YMethod: Least SqiaresDate： 01/15A15 Time: 22:50Sarnple： 1 28Included observatians: 26Weighting series: 1/RESID-*2VariableCoefficientStd. ErrorProb.5 1669381.6G03163.1132250.C0450 1113030.00205654.16207O.COOO图173对所估计的模型再进行