一文读懂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