多元线性回归模型Stata

1、一、邹式检验（突变点检验、稳定性检验）1.突变点检验19852002年中国家用汽车拥有量（yt,万辆）与城镇居民家庭人均可支配收%）数据入（%,元），数据见表6.1年份y（万辆）Xt（元）年份y（万辆）Xt（元）198528.49739.11994205.423496.2198634.71899.61995249.964283198742.291002.21996289.674838.9198860.421181.41997358.365160.3198973.121375.71998423.655425.1199081.621510.21999533.885854199196.041700.

2、62000625.3362801992118.22026.62001770.786859.61993155.772577.42002968.987702.8表6.1中国家用汽车拥有量（yt）与城镇居民家庭人均可支配收入下图是关于yt和xt的散点图:Command：scattervar2var3从上图可以看出，1996年是一个突变点，当城镇居民家庭人均可支配收入突破4838.9元之后，城镇居民家庭购买家用汽车的能力大大提高。现在用邹突变点检验法检验1996年是不是一个突变点。H0：两个字样本（19851995年，19962002年）相对应的模型回归参数相等Hi:备择假设是两个子样本对应的回归参数

3、不等。在19852002年样本范围内做回归Commandregressvar?var3.regressvar-2var3工S?d;KSNueH”c=Qis=£7，电；Q.HOMedel11C90L7.2511169917.28Prob>F=0.0000Residual128203_43168013.03312a-3-quared=0.9027cQK3QUdrSU-U.Toifil1318D26.741777530.S845=tnotMSE-S551VIx<coer.Std.Ezr.tP>lt1S5*Conf.Inu«rv<lvarS.1138341.

4、OO3341S12.130.000C340303.133«3S_cons-113.60C33B,87504-2.53.oio-136,0183-31,19545在回归结果中作如下步骤（邹氏检验）：1、Chow模型稳定性检验（Irtest）用似然比作chow检验，chow检验的零假设：无结构变化，小概率发生结果变化*估计前阶段模型Commandquiregvar2var3in1/11eststoreA*估计后阶段模型；匚口mnn3ndquivar2var3in12/18eststoreC*整个区间上的估计结果保存为AllCcrnmandquiegvar2var3in1/18eststo

5、reAll*用似然比检验检验结构没有发生变化的约束CommandIrtcstAII)(AQstats得到结果如下；.Irteat(AllfAC】,目七瓶七Lilelih.ocdirat:j_e>testLHdii2(=107.74Prob>chi2=0,0000AdsaxLicn:4Al1)nestedin-：ArC)McdelOb911(nullJdfAICBLCAll18-12.3524-105.3303,22147607216.54L4*II-6216548-32.3329268,7657569.5L5SC7-47.78782-19.12435242.2-49942.1417

6、2Wcte:Hftimi_i曰mHincalculatingGIC;与色巨RBICnntfr（如何解释？）2.稳定性检验（邹氏稳定性检验）以表6.1为例，在用19851999年数据建立的模型基础上，检验当把20002002年数据加入样本后，模型的回归参数时候出现显著性变化。*用F-test作chow间断点检验检验模型稳定性* chow检验的零假设：无结构变化，小概率发生结果变化* 估计前阶段模型Commandquiregvar2var3in1/11scalarn1-eM）scalarrss!=erss）* 估计后阶段模型quiregvar2var3in12f18scalarn2-eN)scal

7、arrss2=erss)* 整个区间上的估计结果保存为AllCommandquiregvar2var3in1/18scalark=e(d1mscalarrssr=c(rss* 用F检验检验结构没有发生变化的约束* 计算和显示F检验统计量公式，零假设：无结构变化Cominandscalarf_t«5t=rs5rrss1+rss2)/k+l)/|r?&1盯3行2口1+n2-2xk+1J然后disf_test则得到结其;.disr_test270.716.* F统计量的临界概率CommanddisFtailk+1J4n1tn2-Z*(k+l)JJJesI)然后得到结果.dieFta

8、ll(,(nl+n2-3*(1+1),7.3306-15* F统计量的临界值ComnnancidisinvFtail(k+1)J(nl4n2-2x(l(4l),(.C5J然后得到结果dieinvFtiil(<k+L>r(nl+n2-2*lk41|0OS)3,（如何解释？）二、似然比（LR）检验有中国国债发行总量（DEBTt,亿元）模型如下：DEBT=P°+P1GDR+P2DEFt+P3REPAY+utREPAY表示REPAY)其中GDR表示国内生产总值（百亿元），DEFt表示年财政赤字额（亿元）年还本付息额（亿元）。19802001年数据见表6.2。表6.2国债发行总量D

9、EBT-GDP；、财政赤字额DEFt、年还本付息额（数据198043.0145.17868.928.581991461.4216.178237.14246.81981121.7448.624-37.3862.891992669.68266.381258.83438.57198283.8652.94717.6555.521993739.22346.344293.35336.22198379.4159.34542.5742.4719941175.25467.594574.52499.36198477.3471.7158.1628.919951549.76584.781581.52882.96198

10、589.8589.644-0.5739.5619961967.28678.846529.561355.031986138.25102.02282.950.1719972476.82744.626582.421918.371987223.55119.62562.8379.8319983310.93783.452922.232352.921988270.78149.283133.9776.7619993715.03820.67461743.591910.531989407.97169.092158.8872.3720004180.1894.4222491.271579.821990375.4518

11、5.479146.49190.0720014604959.3332516.542007.73对以上数据进行回归分析：Commandregressvar2va3vaMvar5得到以下结果：/国学后明viirSvar4varSSouxceSSitSUuEilieEo-bs=22ITJ个11。片Made!4323231.2315441077.1Peq1>>F=0,0000Residual46460.7754IB2632.23-squared-Q.Total46371692212208175.BlKn-QtMSE=51,887variCo#f.Std.Ezr.rP>111S5tCon

12、f.Interval.345Z018.15446962r230.0380206732.6697303vari9954D28.031C13I31.490.000,-9289911.0C1S2viifE0797595.04950751TF000071574919637,714.91400421.£72S0.200844-41.207249.S3S21对应的回归表达式为:DEBTt=4.310.35GDPt1.00DEFt0.88REPAY(0.2)(2.2)(31.5)(17.8)2R2=0.999,DW=2.1,F=5735.3现在用似然比（LR）统计量检验约束GDPt对应的回归系数已

13、等于零是否成立。（现在不会）三、Wald检验(以表6.2为例进行Wald检验，对输出结果进行检验。)检验过程如下：1.已知数据如表3.2YX1X211103298351541285-6(1)先根据表中数据估计以下回归模型的方程Y='0''2X2iu2iY=01Xii2X2iUi(2)回答下列问题：叫=质吗？为什么？入2=3吗？为什么?对上述3个方程进行回归分析，结果分别如下:Cemnnandregretvarlvar2得到结果如下:Tsgrtavarlv>r2Source3守d£MSNueib'erofobm=5=3口.820.0115=0.91

14、12=.9017=3.7594ModelPleaid-Lial435-£42.4131445.£1333233U1±pjPab¥F口二与毋1己工邑已AdjH-scjLisizedQgtMSETotal47S4t19.5vaxlCce£.Std.Err,t11白Ccn£.IntezrvH1vaz2_cons£.£-S.fl1.1989373.9429265.55-2.230.0120.1122eiesa-21.2461510.393413.748151Comnnandregressvarlvar3得到结果如下：Cen

15、sCoef.Std.直工工.一口69。皿P>lt|0.000二口口：.工1七口27/三1：-1.44983617.94075.49144419C.020DOOis.eoess19.07292Sonxu百SHdfMSETuusibexof6m=5TFf1专力一，/口勺"7。£FJJMedel7355百2H1475.55fi2Q«Frcb>E-0.DQ02Resitliial2.443753913.81455757R-sqnar&d=CL9945J±aJRsquar'ec.口.9-332Total4764119.5KSE=.3U2

16、55v*r3rer«AB从上述回归结果可知：*#"，月#胤。二元回归与分别对Xi与X2所作的一元回归，其对应的参数估计不相等，主要原因在于Xi与X2有很强的相关性。CommandcorrelatevarZvar3其相关分析结果如下:correlatevar?ttst3（心vax3ar2：1.ODOD-o.3675L,oaoo可见，两者的相关系数为0.9679。Y=%+4Xii+戛*21+UiCommandTegn&ssvar1var?，演3得到结果如下:represavarlvarZvbt3SaiLco«SSdfMSNuinberofob5=5M&

17、dsl47Gr423332223G,21G&G7Prab>F=Dr0033Residual1.5W3EEGN,763333333acruiEed=0.邃aw1q.L.a-nqqqj.Tot；al4764119.55laaE1£SE=Vif18aJ0.E工二*tF>|tI35t匚tin£.var2-1.17777Sl.L1302fi-10«0.401-53C67433.6111S7irar.3-1.344444.2693155-7.220.019-3.103216-.7356732_COTIS21.522224.35525S5.030,0373.

18、1930640，百613s3.表3.3列出了某地区家庭人均鸡肉年消费量Y与家庭月平均收入X,鸡肉价格P1,猪肉价格P2与牛肉价格P3的相关数据。年份Y/千X/P1/（元/P2/（元/P3/（元/年份Y/千P1/（元/P2/（元/P3/（元/克元千克）千克）千克）克X/元千克）千克）千克）19802.783974.225.077.8319924.189113.977.9111.4019812.994133.815.207.9219934.049315.219.5412.4119822.984394.035.407.9219944.0710214.899.4212.7619833.084593.9

19、55.537.9219954.0111655.8312.3514.2919843.124923.735.477.7419964.2713495.7912.9914.3619853.335283.816.378.0219974.4114495.6711.7613.9219863.565603.936.988.0419984.6715756.3713.0916.5519873.646243.786.598.3919995.0617596.1612.9820.3319883.676663.846.458.5520005.0119945.8912.8021.9619893.847174.017.009

20、.3720015.1722586.6414.1022.1619904.047683.867.3210.6120025.2924787.0416.8223.2619914.038433.986.7810.48(1)求出该地区关于家庭鸡肉消费需求的如下模型：lnY=：0ilnX2InR"为lnP2；，%lnP3u(2)请分析，鸡肉的家庭消费需求是否受猪肉及牛肉价格的影响依次生成变量lnvar2先做回归分析，过程如下:lnvar3lnvar4lnvar5lnvar6.gensr-ateLmra.T2=Log-(paHN,.generateLnvarS-Log-tTrirS)generate

21、Invar4Log-(iraT4Lnvar5=Log'tiFaTB)gener-ate：Ln.Ta.T-fi=log-CommandrepressInvar?InvarJEnvar4InvarSInvarG回归结果如下：regressLuvar2LmrarSInvarGSourceSSd.fKSHumberof=53上；吆.±qJ一-口1工一工1Model.T£L174644.19029366而Prob>F=a.QOOORftsxdual.013578211IQ.口口U754345R-aquaxad=0.9025ACLj5(JulaZSd-U.3CCTstal

22、774755fi7B25.04521504EicowMS"=.05747InvarZCcef.Err.P*lt|36*Ccnf.IntervalInvarS.3452563.08254?4.IB0.001.1717549匚4-.5DZ121B,1038305-4.57oooo-,752932-,37125041E1V«E&,1468672,0$?0D1.480.155-.0611368.354B711lrvarG.DS71D45.05985220S70.394-.1225972_een3-.731505-2.460.024-iassasa所以，回归方程为：lnY=-

23、0.7315+0.3463lnX0.5021lnR+0.14691n巳十0.08721nP3(-2.463)(4.182)(-4.569)(1.483)(0.873)由上述回归结果可以知道，鸡肉消费需求受家庭收入水平和鸡肉价格的影响,而牛肉价格和猪肉价格对鸡肉消费需求的影响并不显著。(AIC和SC准则不会算)去掉猪肉价格P2与牛肉价格P3重新进行回归分析Commandregress1rm2lnva3InvaM得出结果如下:片电目衣总日目Invar2InvarSLnvaHSauj：S3d£MSiJuznljeru£cfcsF(32Q)Pish>FsquaredAdjR-

24、aqusred三口ciGM53=52=4972启=0.DODO-0.3«Q3=0.5733=.02763ModelResidual.75-9480334.0152725412.379740167NO,000763627rotal,77475Z87522,03521604Coaf.Std.Eli.rP予1115邑看Conf.IntervalLnvarSLnvaf14u口口沟45154£5-.3727944.024554418.39.063104-5SI.0694201-150.0000.0000.00000327-E043fi71-1.21059S5C27CE241101S5

25、413£(AIC和SC准则不会算)2.某硫酸厂生产的硫酸的透明度指标一直达不到优质要求，经分析透明度低与硫酸中金属杂质的含量太高有关。影响透明度的主要金属杂质是铁、钙、铅、镁等。通过正交试验的方法发现铁是影响硫酸透明度的最主要原因。测量了47组样本值，数据见表3.4。表3.4硫；酸透明度y与铁杂质含量x数据序数XY序数XY1311902560502321902660413341802761524351402863345361502964406371203065257391103169308408132742094210033744010428034762511431103579301

26、243803685251348683787161449803889161550503999201652704076201752504110020185360421002019544443110152054544411015215648451222722565046154202358564721020245852硫酸透明度与铁杂质含量的散点图如下:Ccmnnandscattervar3var2得到以下结果:所以应该建立非线性回归模型。1 .通过线性化的方式估计非线性模型。生成变量：generatey=<l/var3generates父=(l/var21(1)建立倒数模型:Commandregrewyx得到以下结果：,regressyXSmu工ue52d=MSKiiTihezc£lr4£JP