老師好！手里資料給的都是雙門檻模型的畫圖命令，想請教一下單門檻值時是否可以畫LR圖，具體命令應(yīng)該如何寫？謝謝老師撥冗解答。

要靈活運用Stata，目前stata提供很多關(guān)于門檻回歸的命令，包括面板門檻，面板動態(tài)門檻，空間面板門檻等。

我們以官方給出的?xthreg?命令為例。下載這個命令自帶的“hansen1999”數(shù)據(jù)。

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有關(guān)的數(shù)據(jù)導(dǎo)入到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")

**你也可以用調(diào)用其他軟件來畫LR圖，例如R，畫出來要好看很多

往期回顧：

互助問答第186期：學(xué)員提問匯總（12）

互助問答第185期：分組回歸問題

互助問答第184期：分組回歸過程中的問題

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

如果您在計量學(xué)習(xí)和實證研究中遇到問題，請及時發(fā)到郵箱szlw58@126.com，專業(yè)委員會有30多名編輯都會看，您的問題會得到及時關(guān)注！請您將問題描述清楚，任何有助于把問題描述清楚的細(xì)節(jié)都能使我們更方便地回答您的問題，提問細(xì)則參見：實證研究互助平臺最新通知（點擊文末閱讀原文查看詳情）

如果您想成為問題解答者，在幫助他人過程中鞏固自己的知識，請發(fā)郵件至szlw58@126.com（優(yōu)先）或給本公眾號留言或加微信793481976給群主留言，我們誠摯歡迎熱心的學(xué)者和學(xué)生。具體招募信息請參見：實證研究互助平臺志愿者團(tuán)隊招募公告

鮮活的事例更有助于提高您的研究水平，呆板的教科書讓人生厭。如果您喜歡，請?zhí)岢瞿膯栴}，也請轉(zhuǎn)發(fā)推廣！

（歡迎轉(zhuǎn)發(fā)，歡迎分享；轉(zhuǎn)載請注明出處，引用和合作請留言。本文作者擁有所有版權(quán)，原創(chuàng)文章最早發(fā)表于“學(xué)術(shù)苑”。任何侵權(quán)行為將面臨追責(zé)?。?/p>

學(xué)術(shù)指導(dǎo)：張曉峒老師

本期解答人：李后建老師

編輯：統(tǒng)計小妹

統(tǒng)籌：易仰楠

技術(shù)：林毅