李子奈计量经济学课堂笔记.doc

上传人:文库蛋蛋多 文档编号:3824523 上传时间:2023-03-23 格式:DOC 页数:26 大小:589KB
返回 下载 相关 举报
李子奈计量经济学课堂笔记.doc_第1页
第1页 / 共26页
李子奈计量经济学课堂笔记.doc_第2页
第2页 / 共26页
李子奈计量经济学课堂笔记.doc_第3页
第3页 / 共26页
李子奈计量经济学课堂笔记.doc_第4页
第4页 / 共26页
李子奈计量经济学课堂笔记.doc_第5页
第5页 / 共26页
点击查看更多>>
资源描述

《李子奈计量经济学课堂笔记.doc》由会员分享,可在线阅读,更多相关《李子奈计量经济学课堂笔记.doc(26页珍藏版)》请在三一办公上搜索。

1、计量经济学模型的各种检验经济意义检验:检验求得的经济模型的参数估计值的符号和大小是否符合经验和经济理论。统计检验:检验参数估计值的可靠性,包括拟合优度检验、回归效果的检验。拟合优度检验反映了解释变量的变化可以解释被解释变量R2%的变动;回归效果的检验主要检验方程显著性、解释变量是否对被解释变量具有显著性影响。计量经济学检验:包括随机干扰项的异方差性、序列相关性,解释变量的多重共线性,模型设定的偏误性检验。模型的预测检验:检验模型参数估计量的稳定性及其在样本量变化时的灵敏度。一元线性回归模型普通最小二乘估计(ordinary least squares,OLS):使用OLS方法需要满足的假设有:

2、对于随机干扰项,零均值、同方差、无序列相关性、服从正态分布;对于解释变量,应具有非随机性、如果是随机的不能与随机干扰项相关、各解释变量间不存在线性相关性。最大似然估计=最大或然估计(maximum likelihood,ML) 拟合优度检验:检验模型对样本观测值的拟合程度,其统计量是可决系数R2,R2直接由回归结果给出。但是R2只是比较模糊的推测,不能给出严格的统计结论。变量的显著性检验:考察解释变量是否对被解释变量有显著的线性影响,假设参数为0。使用t统计量检验,可以直接依据回归结果给出的P值判断(Pa 时显著)。在一元回归中,变量的显著性检验与方程总体线性的显著性检验一致,因为只有一个解释

3、变量。参数的置信区间:Eviews5中,函数C(1)返回系数参数的估计值,C(2)返回截距参数的估计值,分别等于Eviews5输出结果中的Coefficient。置信度95%的置信区间的上下限为C(i)stderrs(i)*qtdist(0.975, n-2)(i=1,2),其中stderrs(i)(standard errors)为第i个解释变量的标准误差,等于Eviews5输出结果中的Std. Error;stderrs(i)*qtdist(0.975, n-2)(i=1,2)为估计误差。参数置信区间Eviews5操作:把置信区间保存为矩阵,建一个2行3列的矩阵,步骤:(1)使用matri

4、x(2,3) para_conf_intval,创建一个存放置信区间的2行3列的矩阵(2)打开回归方程,使用scalar tinv=qtdist(0.975, regobs-2)创建并保存t分布临界值的标量tinv,因为下面的方法禁止使用二维数组形式的括号,所以上式无法合并到下式(3)打开回归方程,使用para_conf_intval.fill 1,2,C(1)-stderrs(1)*tinv,C(2)-stderrs(2)*tinv,C(1)+stderrs(1)* tinv,C(2)+stderrs(2)* tinv,按先列后行的顺序填充矩阵。矩阵中的结果即为参数的置信区间,形式为(第几个

5、参数,估计值下限,估计值上限),其中参数显示的顺序与回归结果中的变量顺序一致。均值或个别值的估计值或预测值Eviews5操作:(1)使用scalar X0=解释变量的新观测值,将自变量的取值存放到标量X0中(2)由Y0=C(2)+C(1)*X0(X0为解释变量的取值),打开回归方程,使用scalar forecast_value=C(1)+C(2)*X0(注意:C(1)为截距系数),将Y0预测值存放到标量forecast_value中。 均值预测值的置信区间=预测值均值预测值的估计误差,均值预测值的估计误差=t临界值均值预测值的估计误差,均值预测值的估计误差=残差标准误差sqrt(1/n+X0

6、的离差方X的总离差方和)。Eviews5中的均值预测值95%置信区间公式为:Y0qtdist(0.975, n2)sesqrt(1/n+(X0mean(X)2/(nvar(X),其中,X0为解释变量X的样本观测值;函数se返回残差标准误差(也称回归标准误差或随机干扰项标准差,其值等于Eviews5结果中的S.E. of regression,即s2的估计值);(X0mean(X) 2是X0的离差方;nvar(X) 返回X的总离差方和。此问题必须在方程被显示的情况下进行。均值预测值置信区间Eviews5操作:把均值预测值的置信区间保存为矩阵,建一个1行2列的矩阵,步骤:(1)使用matrix(1

7、,2) mean_forecast_intval,创建一个存放置信区间的2行3列的矩阵(2)打开回归方程,使用mean_forecast_intval.fill forecast_value-tinv*se*sqrt(1/regobs+(X0-mean(X)2/(regobs*var(X), forecast_value+tinv*se*sqrt(1/regobs+(X0-mean(X)2/(regobs*var(X),按先列后行的顺序填充矩阵。矩阵中的结果即为参数的置信区间,形式为(第几个参数,估计值下限,估计值上限),其中参数显示的顺序与回归结果中的变量顺序一致。个别值预测值的置信区间=预

8、测值个别值预测值的估计误差,个别值预测值的估计误差=t临界值个别值预测值的估计误差,个别值预测值的估计误差=残差标准误差sqrt(1+1/n+X0的离差方X的总离差方和)。Eviews5中的个别值预测值95%置信区间公式为:Y0qtdist(0.975, n2)sesqrt(1+1/n+(X0mean(X)2/(nvar(X) (注意:sqrt中要加1)。其中,X0为解释变量X的新观测值;函数se返回残差标准误差(也称回归标准误差,其值等于Eviews5结果中的S.E. of regression);(X0mean(X) 2是X0的离差方;nvar(X) 返回X的总离差方和。此问题必须在方程被

9、显示的情况下进行。个值预测值置信区间Eviews5操作:把个别值预测值的置信区间保存为矩阵,建一个1行2列的矩阵,步骤:(1)使用matrix(1,2) single_forecast_intval,创建一个存放置信区间的2行3列的矩阵(2)打开回归方程,使用single_forecast_intval.fill forecast_value-tinv*se*sqrt(1+1/regobs+(X0-mean(X)2/(regobs*var(X), forecast_value+tinv*se*sqrt(1+1/regobs+(X0-mean(X)2/(regobs*var(X),按先列后行的顺

10、序填充矩阵。矩阵中的结果即为参数的置信区间,形式为(第几个参数,估计值下限,估计值上限),其中参数显示的顺序与回归结果中的变量顺序一致。多元线性回归模型矩估计(moment method,MM)工具变量方法(intstrumental variables,IV)广义矩估计法(generalized moment method,GMM)拟合优度检验:统计量是可决系数R2,常使用调整的可决系数(Adjusted R-squared)。赤池信息准则(Akaike information criterion,AIC):用于比较所含解释变量个数不同的多元回归的拟合优度。如果增加变量后使得AIC减少,说明

11、该增加的变量具有解释能力,则在模型中增加该解释变量。AIC由Eviews5直接输出。施瓦茨准则(Schwarz criterion,SC):用于比较所含解释变量个数不同的多元回归的拟合优度。如果增加变量后使得SC减少,说明该增加的变量具有解释能力,则在模型中增加该解释变量。SC由Eviews5直接输出。方程总体线性显著性检验:简称方程显著性检验,假设方程总体不呈线性。使用F统计量(F-statistic)检验,可以直接依据回归结果给出的P值(Prob(F-statistic))判断(Pa 时显著)。变量的显著性检验:考察解释变量是否对被解释变量有显著的线性影响,以决定相应的解释变量能否被保留在

12、模型中,当某解释变量的影响不显著时,应当将该变量从模型中剔除。假设参数为0。使用t统计量(t-Statistic)检验,可以直接依据回归结果给出的P值判断(Pa 时显著)。参数的置信区间:Eviews5中,使用函数C(i)获得第i个参数的估计值(i=1,2,k+1),等于Eviews5输出结果中的Coefficient,C(k+1)为截距的估计值,k为解释变量的个数,不含常数项,k+1为参数的个数。置信度95%的置信区间的上下限等于C(i)stderrs(i)qtdist(0.975, n-k)(i=1,2,k+1),其中stderrs(i)(standard errors)为第i个解释变量或

13、截距的标准误差,等于Eviews5输出结果中的Std. Error;stderrs(i)qtdist(0.975,n-k)为估计误差。参数的置信区间Eviews5操作:把置信区间保存为矩阵,建一个k行3列的矩阵,步骤:(1)打开回归方程,使用scalar n=regobs将样本容量存放到标量n中(2)使用scalar k=ncoef将参数个数存放到标量k中(3)使用scalar tinv=qtdist(0.975, n-k)创建并保存t分布临界值的标量tinv,因为下面的方法禁止使用二维数组形式的括号(4)使用matrix(k,3) para_conf_intval,创建一个存放置信区间的k行

14、3列的矩阵(5)回归方程具有3个参数时,可使用para_conf_intval.fill 1,2,3,C(1)-stderrs(1)*tinv,C(2)-stderrs(2)*tinv,C(3)-stderrs(3)*tinv,C(1)+stderrs(1)*tinv,C(2)+stderrs(2)*tinv,C(3)+stderrs(3)*tinv,按先列后行的顺序填充矩阵。矩阵中的结果即为参数的置信区间,形式为(第几个参数,估计值下限,估计值上限),其中参数显示的顺序与回归结果中的变量顺序一致。均值或个别值的估计值或预测值:由Y0=C(k+1)C(1)X10C(2)X20C(k)Xk0,其

15、中,X10,X20,Xk0为解释变量X1,X2,Xk的观测值;C(k+1)为常数项;k为解释变量的个数,不含常数项,故参数个数为k+1。均值与个值的估计值Eviews5操作:(1)使用scalar x10=X1的观测值,scalar x20=X2的观测值scalar xk0=Xk的观测值,将k个解释变量的观测值存放到标量中(2)使用matrix forecast_value=C(1)+C(2)*X10+C(3)*X20+C(k)*Xm0(注意:k=ncoef为参数个数,m=k-1),将预测值Y0存放到矩阵forecast_value中。均值预测值的置信区间:多元回归的置信区间较麻烦,需使用矩阵

16、。均值的置信区间Eviews5操作:(1)使用命令:group gx 1 x1 x2 (注意:只用空格间隔,开头没有等号,1为常数列向量,X1,X2为两个解释变量序列),将常数项和解释变量的列向量合并为n行k列的组(不是矩阵,此处的k是参数个数)(2):使用命令:stom(gx, x)将组gx转化为矩阵X(3)使用命令:matrix(1, k) x0(此处的k=ncoef是参数个数)创建一个1行k列的矩阵x0,表示新的观测值矩阵(4)使用命令:x0.fill 1, X10, X20,(注意有逗号,有常数项列向量1,Xk0表示第k个解释变量的观测值)把样本观测值X0用行向量矩阵表示为(1 X10

17、 X20 )(5)使用matrix(1,2) mean_conf_intval,创建一个1行2列的0矩阵(6)使用colplace(mean_conf_intval, forecast_value-tinv*se*sqrt(X0*inverse(transpose(X)*X)*transpose(X0),1)(注意:forecast_value是矩阵,不是上面的预测值)。将均值预测值置信度95%的置信区间的下限写入矩阵mean_conf_intval入的第1个单元格(7)使用colplace(mean_conf_intval, forecast_value+tinv*se*sqrt(X0*in

18、verse(transpose(X)*X)*transpose(X0),2)(注意:forecast_value是矩阵,不是上面的预测值;最后一个逗号可能会出问题)。将均值预测值置信度95%的置信区间的上限写入矩阵mean_conf_intval入的第2个单元格。个别值预测值的置信区间:多元回归的置信区间较麻烦,需使用矩阵。个值的置信区间Eviews5操作:(1)使用matrix(1,2) single_conf_intval,创建一个1行2列的0矩阵(2)使用matrix one=1创建一个1行1列的单位矩阵(3)使用colplace(single_conf_intval, forecast

19、_value-tinv*se*sqrt(one+X0*inverse(transpose(X)*X)*transpose(X0),1)(注意:forecast_value是矩阵,不是上面的预测值;one不能改为1)。将均值预测值置信度95%的置信区间的下限写入矩阵mean_conf_intval入的第1个单元格(4)使用colplace(single_conf_intval, forecast_value+tinv*se*sqrt(one+X0*inverse(transpose(X)*X)*transpose(X0),2)(注意:forecast_value是矩阵,不是上面的预测值)。将均值

20、预测值置信度95%的置信区间的上限写入矩阵mean_conf_intval入的第2个单元格。可化为线性的多元非线性回归模型倒数模型:1/Q=a+b/P+m,用Y=1/Q,X=1/P置换,得到Y=a+bX+m多项式模型:s=a+br+cr2+m,用X1=r,X2=r2置换,得到Y=a+bX1+cX2+m幂函数模型:Q=AKaLem,取对数后,lnQ=lnA+a lnK+lnL+m,Q=f(X1,X2,Xk)的幂函数形式为Q=AX11X22Xkk,注意使用全体样本取对数。无约束回归模型:Q=AX1P12 P03,取对数后lnQ=0+1lnX+2lnP1+3lnP0+m,其中0=lnA,然后进行多元

21、线性回归。注意使用全体样本取对数。有约束回归模型:Q=AX1P12 P03,s.t.1+2+3=0,取对数后同上,将3换做-1-2,得到lnQ=0+1ln(X/ P0)+2ln(P1/ P0)+ m。施加的约束为真的受约束回归模型与无约束回归模型具有相同的解释能力,但约束是否为真需要进行检验。线性约束检验:即线性约束的零阶齐次性假设检验,假设方程受线性约束。如对模型URR:Y=0+1X1+2X2+kXk+m,施加1+2=1,k-1=k的约束变为RR模型:Y=0+1X1+(1-1)X2+k-1Xk-1+kXk+m(样本量未变)。检验约束条件1+2=1,k-1=k是否为真的统计量为F=(RSSR-

22、RSSu)/(ku-kR)/(RSSU/(n-ku-1)F(ku-kR,n-ku-1),其中,RSSR为有约束条件(Restrict)的残差平方和,RSSU为无约束回归模型(Unrestrict)的残差平方和,kR为有约束回归模型中的解释变量的个数,kU为无约束回归模型中的解释变量的个数,不含常数项,kU-kR为约束条件的个数,n为受约束模型的样本量。如果RSSR-RSSU值较大,则F较大,说明约束条件无效。计算置信度为95%的F临界值的公式是qfdist(0.95,1,10)。线性约束检验稍复杂的Eviews5操作:注意样本量不变,实质的解释变量个数不变。(1)进行无约束回归,求scalar

23、 n=regobs,scalar rssu=ssr,scalar ku=ncoef-1;(2)按零阶齐次性表达式进行有约束回归,求scalar rssr=ssr,scalar kr=ncoef-1;(3)计算F统计量scalar linear_restrict_F_stat=(rssr-rssu)/(ku-kr)/(rssu/(n-ku-1);(4)使用scalar linear_restrict_P_val=1-cfdist(linear_restrict_F_stat, ku-kr, rssu/(n-ku-1)求P值;(5)如果linear_restrict_P_val0.05,不能拒绝约

24、束不成立的假设,即不能拒绝解释变量不具有零阶齐次特性,亦即不能拒绝线性约束成立。线性约束检验Eviews5菜单操作:进行无约束回归viewCoefficient TestsWald-Coefficient Restrictions输入约束条件根据结果中的F统计量值对应的P值判定,若P0.05则不能拒绝约束成立的假设。解释变量取舍:属于线性约束问题。假设变量没有解释能力。如在模型RR:Y=0+1X1+2X2+kXk+m中添加变量后变成模型URR:Y=0+1X1+1X2+kXk +k+1Xk+1+k+pXk+p+m,RR即是在URR中施加了k+1=k+2=k+p=0的(线性)约束条件(解释变量变化

25、,但样本量未变)。检验约束条件k+1=k+2=k+p=0是否为真的统计量为F=(RSSR-RSSU)/(kU-kR)/(RSSU/(n-k-kU+kR-1)F(kU-kR,n-k-kU+kR-1),一般使用其等价式:F=(R2u-R2R)/(ku-kR)/(R2u/(n-k-ku+kR-1)F(ku-kR,n-ku-1),其中,R2u为无约束回归模型的可决系数,R2R为有约束回归模型的可决系数,kU为无约束回归模型中的解释变量的个数, kR为有约束回归模型中的解释变量的个数,n为受约束模型的样本量。如果F计算值比临界值大,说明约束条件无效,即添加的解释变量有解释能力,或者说P0.05,增加解释

26、变量才是合适的。解释变量取舍Eviews5操作:注意样本量不变,解释变量个数有变化。(1)进行解释变量个数较少的有约束回归,求scalar n=regobs,scalar rr_square=r2,scalar kr=ncoef-1;(2)进行解释变量个数较多的无约束回归,求scalar ru_square=r2,scalar ku=ncoef-1;(3)根据统计量F=(R2u-R2R)/(ku-kR)/(1-R2u/(n-ku-1),计算F统计量scalar add_x_F_stat=(ru_square-rr_square)/(ku-kr)/(1-ru_square)/(n-ku-1);(

27、4)scalar add_x_P_val=1-cfdist(add_x_F_stat, ku-kr, n-ku-1)计算P值;(5)如果add_x_P_val0.05,说明新增的解释变量需要引入模型。线性约束中的(邹氏)参数的稳定性检验(Chow test for parameter stablity):在使用不同样本量的连续的序列(无重合的多期)中,如果结构发生变化,会导致参数不稳定,从而模型的预测能力会下降。假设参数稳定。无约束模型URR:前期Y=0+1X1+2X2+kXk+m1及后期Y=a0+a1X1+a2X2+akXk+m2。如果序列结构没有发生变化,则=a(矩阵式),故受约束模型RR

28、:前期Y=0+1X1+2X2+kXk+m1及后期Y=0+1X1+2X2+kXk+m2。检验=a是否为真的统计量为F=(RSSR-RSSU)/(k+1)/(RSSU/(n1+n2-2(k+1)F(k+1,n1+n2-2(k+1),常用其等价式:F=(RSSR-(RSS1+RSS2)/(k+1)/(RSS1+RSS2)/(n1+n2-2(k+1)F(k+1,n1+n2-2(k+1)。线性约束中的邹氏参数的稳定性检验稍麻烦的Eviews5操作:注意样本量有变化,解释变量个数不变。(1)使用前期小样本进行无约束回归,求scalar n1=regobs,scalar rss1=ssr,scalar k=

29、ncoef-1;(2)使用后期小样本进行无约束回归,scalar n2=regobs,scalar rss2=ssr;(3)使用假设前后期参数一致的合并大样本进行有约束回归,scalar rssr_all=ssr;(4)计算F统计量scalar chow_para_stabl_F_stat=(rssr_all-rss1-rss2)/(k+1)/(rss1+rss2)/(n1+n2-2*(k+1);(5)scalar chow_para_stabl_P_val=1-cfdist(chow_para_stabl_F_stat,k+1, n1+n2-2*(k+1),计算P值;(6)如果chow_pa

30、ra_stabl_P_valk,即后期样本量不小于参数个数。线性约束中的邹氏预测检验(Chow test for predictive failure):其思想是用第一时间段的估计参数对第二时间段的样本进行预测,若预测误差较大,说明参数发生了变化。假设参数稳定。假设参数不同的无约束模型(向量形式)URR:前期Y1=X1+m1及后期Y2=X2a+m2=X2a+X2-X2+m2=X2+g+m2(记g=X2(a-),如果序列结构没有发生变化,即=a,则假设参数不变的受约束模型RR:前期Y1=X1+m1及后期Y2=X2+m2。检验=a是否为真的统计量为F=(RSSR-RSSU)/n2)/(RSSU/(

31、n1-k-1)F(k+1,n1+n2-2(k+1)。线性约束中的邹氏预测检验稍麻烦的Eviews5操作:注意样本量有变化,解释变量个数不变。(1)使用前期小样本进行无约束回归,求scalar n1=regobs,scalar rss1=ssr;(2)使用前后期参数一致的合并大样本进行有约束回归,scalar rssr_all=ssr;(3)计算F统计量scalar chow_predict_F_stat=(rssr_all-rss1)/n2)/(rss1/(n1-k-1)(可手工录入后期样本量n2);(4)scalar chow_predict_P_val=1-cfdist(chow_pred

32、ict_F_stat,n2,n1-k-1),计算P值;(5)如果chow_predict_P_val0.05,则不能拒绝约束成立。放宽基本假定的模型计量经济学模型基本假定违背:随机干扰项序列存在异方差性、随机干扰项序列存在序列相关性、解释变量之间存在多重共线性、解释变量是随机变量且与随机干扰项相关、模型设定有偏误、解释变量的方差随样本容量的增加而增加。使用OLS方法进行回归时,不能有对基本假定的违背,即不能出现上述情况之一。计量经济学检验:对回归模型是否满足OLS的基本假定进行检验。异方差性异方差性(Heteroskedasticity):总体中每个样本单位的随机干扰项的方差不相同。异方差意味

33、着未进入模型的随机干扰项的误差会随着解释变量的变化而变化;异方差存在时;因变量取值的差异会随解释变量的不同而比较大,导致样本整体看来缺乏规律性。截面数据常存在异方差性。一般使用加权最小二乘法进行修正。图示异方差检验:观察是否存在明显的散点扩大或缩小或复杂趋势即不在一个带状区域的Y-X散点图,或观察是否是一条水平直线的残差平方-X的散点图进行检验。Park(帕克)检验和Gleiser(戈里瑟)检验:假设回归模型不存在异方差性。基本思想是,以残差绝对值或残差平方为被解释变量,某一解释变量Xj为解释变量,进行回归,对回归方程进行显著性检验,若检验显著,则存在异方差性。Park检验:假设回归模型不存在

34、异方差性。(1)使用OLS方法建立回归方程,用series resid_ols=resid保存残差序列;(2)使用ls log(resid_ols2) c log(x)求回归方程,若回归方程的检验与参数检验均显著,则存在异方差性。Gleiser检验:假设回归模型不存在异方差性。选择X不同的函数形式,以| resid |为被解释变量,使用OLS方法依次建立回归方程,用ls abs(resid_ols) c x,ls abs(resid_ols) c sqrt(x),ls abs(resid_ols) c 1/x,若某个回归方程的检验与参数检验均显著,则存在异方差性;有的方程检验也许不能拒绝同方差

35、性。GQ(Goldfeld-Quandt=戈德菲尔德-匡特)检验:假设回归模型不存在异方差性。(1)使用OLS方法建立回归方程,用scalar n=regobs保存总样本量,scalar k=ncoef-1保存解释变量个数,使用series resid_ols=resid保存残差序列。找出最大的系数对应的解释变量,它是差异的主要来源,可能是引起异方差的主要因素;(2)使用sort(d) 主要解释变量,对所有序列降序排序(默认升序);(3)打开回归方程,使用scalar n_sub=ceiling(3*n/8)计算子样本容量(n是样本总量,ceiling为上取整函数);(4)对排序后数据较大的子

36、样本进行回归(Quick-Estimate Equation,在Sample文本框中输入1 n_sub的值),并求残差平方和scalar RSS1=ssr;(5)使用sort 主要解释变量,对所有序列升序排序,对排序后数据较小的子样本进行回归(Quick-Estimate Equation,在Sample文本框中输入1 n_sub的值),并求残差平方和scalar RSS2=ssr;(6)利用两个子样的残差平方和之比构造统计量F=RSS1/RSS2F(n_sub-k-1,n_sub-k-1)(k是解释变量个数)并进行异方差检验:使用scalar G_Q_F_stat=RSS1/RSS2计算F统

37、计量,使用scalar G_Q_P_val=1-cfdist(g_q_f_stat,n_sub-k-1,n_sub-k-1)计算F统计量的P值;(7)如果G_Q_P_val0.05,则说明总体随机干扰项存在异方差。White检验:假设回归模型不存在异方差性。无需排序,使用样本总量进行含交叉项辅助回归(即以残差平方和为被解释变量,各原解释变量的一次幂项、二次幂项、两变量的交叉乘积项为解释变量),根据辅助回归的结果,nR2c2(k),其中n为样本容量,R2为辅助回归的可决系数,k为解释变量的个数,不含常数项。(1)含交叉项的White检验:其思想是先作含交叉项的辅助回归,即使用(各项)残差平方re

38、sid_square(因变量)、Xi、Xi2及交叉项Xi*Xj(或resid_square、lnXi、(lnXi)2、lnXi*lnXj)进行回归,回归的结果使用怀特检验,统计量为辅助回归的结果中的Obs*R-squaredc2(k),根据其对应的P值判断随机干扰项是否存在异方差。Eviews5菜单操作:选取原变量为Y、Xi(或lnY、lnXi)作OLS回归viewResidualTestsWhite Heteroskedasticity(cross terms)LM=Obs*R-squared即为怀特统计量c2(k),其后即为对应的P值,若P0.05,说明总体随机干扰项存在异方差;含交叉项怀

39、特检验结果中还给出了辅助回归结果;(2)不含交叉项的White检验:其思想是先作去掉交叉项的辅助回归(样本量为总样本量),即使用(各项)残差平方resid_square(因变量)、Xi、Xi2(或resid_square、lnXi、(lnXi)2)进行回归,回归的结果使用怀特检验,统计量为Obs*R-squaredc2(k),根据其对应的P值,可判断总体随机干扰项是否存在异方差。Eviews5菜单操作:选取原变量为Y、Xi(或lnY、lnXi)作OLS回归viewResidual TestsWhite Heteroskedasticity(no cross terms)LM=Obs*R-squ

40、ared即为怀特统计量c2(k),其后即为对应的P值,若P0.05,则说明已不存在异方差性。在实际使用中,通常不进行异方差检验,而是直接选择加权最小二乘法。若确实存在异方差性,则会被消除;若不存在,则等价于最小二乘法。序列相关性序列相关性(Serial Correlation):随机干扰项不相互独立,总体中每个样本单位的随机干扰项相关。如果一样本单位只与上一个样本单位相关,称为一阶序列相关或自相关,其形式为mi=rmi-1+ei。很多情况下,时间序列数据中解释变量以外的因素在时间上存在连续性,带来对被解释变量的影响的连续性,从而导致序列相关性。未进入回归模型的随机干扰因素受到历史的影响使不考虑

41、历史影响的回归误差加大。时间序列数据常存在序列相关性。常使用广义最小二乘法、广义差分法修正。图示序列相关性检验:观察残差与时间的散点图,或观察本期残差和上期残差的散点图。回归检验法:假设回归模型不存在序列相关性。使用et=ret-1+et、et=r1et-1+r2et-2+et一个一个的建立方程,并分别进行检验。如果某个函数形式存在,则原模型存在序列相关性。DW=杜宾-瓦森(Durbin-Watson)检验:假设不存在一阶自回归。只能检验一阶自相关,不能使用滞后被解释变量。根据样本量n、参数个数k查DW上下界表,若DW(0,dl),则存在正相关;若DW(4-du,4-dl),则存在负相关;若DW(du,4-du),则无自相关;若DW( dl,du),则不能确定。若DW不能判定,可借助LM判定。一般,DW在2附近时不存在序列相关性,偏离2越远,序列相关性越严重。拉格朗日乘数(Lagrange multiplier,LM)检验:也称Breusch-Godfrey检验或GB检验。假设回归模型不存在序列相关性。用于检验高阶序列相关,以及存在滞后被解释变量的情形。作OLS回归viewResidual TestsSerial Correlatio

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > 办公文档 > 其他范文


备案号:宁ICP备20000045号-2

经营许可证:宁B2-20210002

宁公网安备 64010402000987号