一文读懂stata合并多个图示(附合并面板门槛LR统计量对应的图形)
本文简单介绍在Stata中如何将多个图示合并为一个整体图示。具体操作如下:
语法格式为:
graph combine name [name ...] [, options]
案例操作命令如下(以系统自带数据为例):
sysuse uslifeexp
line le_male year, saving(male)
line le_female year, saving(female)
gr combine male.gph female.gph
结果显示如下:
接下来本文以门槛回归模型为例,为大家展示结果展示面板门槛LR统计量对应的图形,首先回顾一下门槛回归。
门槛回归的案例
use hansen1999
Estimate a single-threshold model
xthreg i q1 q2 q3 d1 qd1, rx(c1) qx(d1) thnum(1) trim(0.01) grid(400) bs(300)
Estimate a triple-threshold model given the estimated result above
xthreg i q1 q2 q3 d1 qd1, rx(c1) qx(d1) thnum(3) trim(0.01 0.01 0.05) grid(400) bs(300 300 300)
输出结果包括四个部分。第一部分输出门限估计值和自举法的结果。第二部分列表输出门限值及置信区间,Th-1代表单一门限估计值,Th-21 和Th-22代表双门限回归的两个估计值,有时Th-21和Th-1相同。第三部分列出了门限检验,包括RSS、MSE、F统计量及概率值,以及10%、5%、1%的置信水平。第四部分是固定效应回归结果。
在3个门限值情形下,画LR统计量对应的图形
_matplot e(LR21), columns(1 2) yline(7.3523, lpattern(dash)) connect(direct) recast(line) ytitle('LR Statistics') xtitle('First Threshold') name(LR21)
_matplot e(LR22), columns(1 2) yline(7.3523, lpattern(dash)) connect(direct) recast(line) ytitle('LR Statistics') xtitle('2nd Threshold Parameter') name(LR22)
_matplot e(LR3), columns(1 2) yline(7.3523, lpattern(dash)) connect(direct) recast(line) ytitle('LR Statistics') xtitle('3rd Threshold Parameter') name(LR3)
graph combine LR21 LR22 LR3, cols(1)
·
·
·
END
l 计量经济学服务中心 l