《异方差检验的套路(共5页).docx》由会员分享,可在线阅读,更多相关《异方差检验的套路(共5页).docx(5页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上异方差检验的套路:(1)Breusch-pagan法:reg y x1 x2 x3predict u,residgen usq=u2reg usq x1 x2 x3求F值display R/(1-R)*n2/n1(n1表示分子除数,n2表示分母除数)display Ftail()求LM值display R*n(n表示总样本量)display chi2tail()(2)white法:reg y x1 x2 x3predict u,residgen usq=u2predict ygen ysq=y2reg usq y ysq求F值display R/(1-R)*n2/n1
2、(n1表示分子除数,n2表示分母除数)display Ftail()求LM值display R*n(n表示总样本量)display chi2tail()(3)必要补充F值和LM值转换为P值的命令:display Ftail(n1,n2,a)(利用F值求p值,n1表示分子除数,n2表示分母除数,a为F值)display chi2tail(n3,b)(利用LM值求p值,n3表示自由度的损失量,一般等于n1,b为LM值) 异方差的纠正WLS(weighted least square estimator)(1)基本思路:reg y x1 x2 x3 aw=x1(将x1作为异方差的来源,对方程进行修正
3、)上式相当于:reg y/(x10.5) 1/(x10.5) x1/(x10.5) x2/(x10.5) x3/(x10.5),noconstant(2)纠正异方差的常用套路(构造h值)reg y x1 x2 x3predict u,residgen usq=u2gen logusq=log(usq)reg logusq x1 x2 x3predict ggen h=exp(g)reg y x1 x2 x3 aw=1/h异方差hausman检验:reg y x1 x2 x3est store A(将上述回归结果储存到A中)reg y x1 x2 x3 aw=1/hest store Bhaus
4、man A B当因变量为对数形式时(log(y))如何预测yreg logy x1 x2 x3predict kgen m=exp(k)reg y m,noconstantm的系数为iy的预测值=iexp(k)方差分析:一元方差分析anova y g1 / g1|g2 /(g*表示不同分类变量,计算g1和交互项/ g1|g2 /这两种分类的y值是否存在组内差异)anova y d1 d2 d1*d2(d*表示虚拟变量,计算d1、d2和d1*d2的这三种分类的y值是否有组内差异) anova y d1 d2 x1 d2*x1, continuous(x1)(x*表示连续的控制变量)多元方差分析w
5、ebuse jawmanova y1 y2 y3 = gender fracture gender*fracture(按性别、是否骨折及二者的交互项对y1、y2和y3进行方差分析)manova y1 = gender fracture gender*fracture(相当于一元方差分析,以y1为因变量)webuse nobetweengen mycons = 1manova test1 test2 test3 = mycons, noconstantmat c = (1,0,-1 0,1,-1)manovatest mycons, ytransform(c)进行多元回归的方法:多元回归分析:(
6、与mvreg相同)foreach vname in y1 y2 y3 (确定y变量组vname)reg vname x1 x2 x3 (将y变量组中的各个变量与诸x变量进行回归分析,注意vname的标点符号) 上式等价于:mvreg y1 y2 y3 = x1 x2 x3reg3命令:(1)简单用法:reg3 (y1 = x1 x2 x3) (y2 = x1 x3 x4) (y3 = x1 x2 x5)测试y1 coefs = 0test y1测试不同回归中相同变量的系数:test y1=y2=y3, commontest (y1=y2) (y1=y3), common constant(constant表示包含截距项)(2)用reg3进行2SLSreg3 (y1 = y2 x1 x2) (y2 = y1 x4),2sls(2)用reg3进行OLSreg3 (y1 = y2 x1 x2) (y2 = y1 x4),ols对两个回归结果进行hausman检验:专心-专注-专业