互助问答第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，画出来要好看很多

往期回顾：

互助问答第186期：学员提问汇总（12）

互助问答第185期：分组回归问题

互助问答第184期：分组回归过程中的问题

互助问答第183期：三重差分模型

如果您在计量学习和实证研究中遇到问题，请及时发到邮箱szlw58@126.com，专业委员会有30多名编辑都会看，您的问题会得到及时关注！请您将问题描述清楚，任何有助于把问题描述清楚的细节都能使我们更方便地回答您的问题，提问细则参见：实证研究互助平台最新通知（点击文末阅读原文查看详情）

如果您想成为问题解答者，在帮助他人过程中巩固自己的知识，请发邮件至szlw58@126.com（优先）或给本公众号留言或加微信793481976给群主留言，我们诚挚欢迎热心的学者和学生。具体招募信息请参见：实证研究互助平台志愿者团队招募公告

鲜活的事例更有助于提高您的研究水平，呆板的教科书让人生厌。如果您喜欢，请提出您的问题，也请转发推广！

（欢迎转发，欢迎分享；转载请注明出处，引用和合作请留言。本文作者拥有所有版权，原创文章最早发表于“学术苑”。任何侵权行为将面临追责！）

学术指导：张晓峒老师

本期解答人：李后建老师

编辑：统计小妹

统筹：易仰楠

技术：林毅