互助问答第187期：双门槛模型解惑

老师好！手里资料给的都是双门槛模型的画图命令，想请教一下单门槛值时是否可以画LR图，具体命令应该如何写？谢谢老师拨冗解答。

要灵活运用Stata，目前stata提供很多关于门槛回归的命令，包括面板门槛，面板动态门槛，空间面板门槛等。

我们以官方给出的 xthreg 命令为例。下载这个命令自带的“hansen1999”数据。

use hansen1999, clear

cap log close

set more off

local indep "q1 q2 q3 d1 qd1"

local condition "grid(400) trim(0.01) bs(300)"

xthreg i `indep',  rx(c1) qx(d1) thnum(1) `condition'

mat d=e(LR)

**提前在c盘里面建立一个名为output的文件夹

cd "C:\output"

putexcel set "C:\output\LR.xls", modify

putexcel A1=mat(d)

mat drop d

**把LR有关的数据导入到excel里面，你想用什么软件画图都可以了

import excel "C:\output\LR.xls", sheet("Sheet1")

rename A LR

rename B threshold

twoway (connected LR threshold,  ms(O oh) connect(direct)  msize(small) mlabp(0) mlabs(zero)), yline(7.35,lc(pink) lp(dash)) xla(,nogrid) yla(,nogrid)  title("") xtitle("Threshold") ytitle("LR Statistics")

**你也可以用调用其他软件来画LR图，例如R，画出来要好看很多