好的，谢谢黄老师

_matplot e(LR), columns(1 2) yline(7.35, lpattern(dash)) connect(direct) msize(small) mlabp(0) mlabs(zero) ytitle("LR Statistics") xtitle("FIRST Threshold") recast(line) 单门槛画图代码
_matplot e(LR22), columns(1 2) yline(7.35, lpattern(dash)) connect(direct) msize(small) mlabp(0) mlabs(zero) ytitle("LR Statistics") xtitle("Second Threshold") recast(line) 双门槛画图代码
老师麻烦您帮我看下，为什么我的单门槛画出来的图为什么是这个样子，我刚刚问同学，她说建议按照单门槛模型来解释，希望您题出建议

黄老师您好，我有个问题想请教一下，我用xthreg加vce(robust)命令，变量结果都挺好，但是结果中的f检验显示是·（空白？）如图一所示。然后我又试了一下把vce(robust)选项删去，f检验结果就有了，变量结果也挺好，如图2所示。想请问一下老师，这是怎么回事啊？我应该怎么调整来使得加vce(robust)选项的情况下，f检验能正常出现呢（因为我论文结果里需要列出f检验的结果）。期待您的回复，谢谢！
图1：
xthreg y3 is m infr fd gov,rx(mot) qx(rd) thnum(1) grid(400) trim( 0.02 ) bs(300) vce(robust)
Estimating  the  threshold  parameters:   1st ......  Done
Boostrap for single threshold
.................................................. +   50
.................................................. +  100
.................................................. +  150
.................................................. +  200
.................................................. +  250
.................................................. +  300

Threshold estimator (level = 95):
-----------------------------------------------------
     model |    Threshold         Lower         Upper
-----------+-----------------------------------------
      Th-1 |       0.2100        0.2000        0.2200
-----------------------------------------------------

Threshold effect test (bootstrap = 300):
-------------------------------------------------------------------------------
Threshold |       RSS        MSE      Fstat    Prob   Crit10    Crit5    Crit1
-----------+-------------------------------------------------------------------
    Single |    0.0661     0.0004      49.16  0.0000  20.7404  23.2570  29.3277
-------------------------------------------------------------------------------
Fixed-effects (within) regression               Number of obs      =       195
Group variable: i                               Number of groups   =        13

R-sq:  within  = 0.8967                         Obs per group: min =        15
       between = 0.0026                                        avg =      15.0
       overall = 0.1409                                        max =        15

                                                F(6,12)            =         .
corr(u_i, Xb)  = -0.0308                        Prob > F           =         .

                                     (Std. Err. adjusted for 13 clusters in i)


图2：
xthreg y3 is m infr fd gov,rx(mot) qx(rd) thnum(1) grid(400) trim( 0.02 ) bs(300)
Estimating  the  threshold  parameters:   1st ......  Done
Boostrap for single threshold
.................................................. +   50
.................................................. +  100
.................................................. +  150
.................................................. +  200
.................................................. +  250
.................................................. +  300

Threshold estimator (level = 95):
-----------------------------------------------------
     model |    Threshold         Lower         Upper
-----------+-----------------------------------------
      Th-1 |       0.2100        0.2000        0.2200
-----------------------------------------------------

Threshold effect test (bootstrap = 300):
-------------------------------------------------------------------------------
Threshold |       RSS        MSE      Fstat    Prob   Crit10    Crit5    Crit1
-----------+-------------------------------------------------------------------
    Single |    0.0661     0.0004      49.16  0.0000  20.7404  23.2570  29.3277
-------------------------------------------------------------------------------
Fixed-effects (within) regression               Number of obs      =       195
Group variable: i                               Number of groups   =        13

R-sq:  within  = 0.8967                         Obs per group: min =        15
       between = 0.0026                                        avg =      15.0
       overall = 0.1409                                        max =        15

                                                F(7,175)           =    217.10
corr(u_i, Xb)  = -0.0308                        Prob > F           =    0.0000

我从来不 care 这个 F 值，我也不觉得它有什么重要性！

谢谢老师

老师好，我在检验三重门槛的时候，发现p值都是<0.1的但是下面图为什么少了一列2呢？还有老师，得出的三个门槛估计值的第三个和第一个相等这个是意味着什么呢，是否需要舍弃？希望老师可以指点一下

老师，做门槛回归时，发现单门槛检验的P值为0.09，是否要拒绝，直接报告无门槛效应存在？因为看到有些文章把10%也算显著，但不如0,05以下效果好

这由你决定 (我个人会认为有)。

黄老师您好，请问使用xthreg指令做固定效应门限回归，不想控制个体，只想控制行业该如何操作呢？

没办法！

好的，非常感谢！

黄老师您好，我看了您说的门槛值的计算方法，那如果在论文中简单提及一下门槛值的计算方法/原理应该怎么说啊？

黄老师您好，我看了您说的门槛值的计算方法，那如果在论文中简单提及一下门槛值的计算方法/原理应该怎么说啊？

你若是要在论文中叙述，请看一些相关"应用文章"之说明，这才是正确之道！

黄老师您好，stata 15 中的threshold 命令，只能用于做时间序列数据的门槛回归，您知道有什么命令可以做横截面数据的门槛回归吗？谢谢您！

请看看 https://www.ssc.wisc.edu/~bhansen/progs/progs_threshold.html。

非常感谢您的回复！我之前也尝试了Hansen的stata code，但是拿来run我自己的数据（大约2000个样本），发现速度非常慢，很久也得不到结果，不知是什么原因？如果用Hansen的matlab code，是否可能会好一些？

或許會快一點！

楼主真是个好人！非常好的资料！感谢楼主！

老师你好！我想用双向固定效应，代码是xi:xthreg y x1 x2 i.year,,rx(x3) qx(x3) thnum(1) grid（400） trim(0.01) bs(300),结果报错thest():  3200  conformability error
                thestm():     -  function returned error
                 <istmt>:     -  function returned error
请问是什么原因呢？双向固定效应和固定效应的门槛检验结果和回归结果会有很大出入吗？

老师，能否把第二个案例压缩包和第三个pdf给我发下一下，那个现在不了了，非常感谢！！！qq邮箱979534884@qq.com

黄老师，您好！我想请问一下面板门槛回归是只能用在固定效应模型中吗？那豪斯曼检验出基准回归应该使用随机效应模型，还可以进一步做门槛回归吗？

讲过很多次了，我只用 FE 相关方法。

好的，我还想再麻烦问下黄老师，如果要做分区域的面板门槛模型，stata命令是什么呢？我用的是全国30个省际面板数据，用xthreg进行门槛回归，结果显示存在单门槛值，现在还想进一步分区域做，就生成了区域虚拟变量，并用keep if group=1 (1表示东部）保留东部地区数据，然后用相同的门槛命令，但是结果出现了
（ thestm():  3301  subscript invalid
                 <istmt>:     -  function returned error
），不知道是怎么回事。

感谢楼主的分享，谢谢！

xthreg example下载不了了，能不能邮箱发一份，感谢楼主 646550590@qq.com

黄老师，我检验门槛存在，但过了门槛之后不显著，这样的话说明依旧不存在门槛么？谢谢黄老师~期待黄老师回复

存在，请将指令与节果发出来，以供判断！

黄老师，您好，我刚学门限模型，想请问一下，门限模型是否可以有两个核心解释变量？有些文献中有两个核心解释变量，请问是用什么程序做的?

老师，您好，我在下载xthreg的时候，输入search xthreg，出现这样的情况，请问这个怎么解决呢