请问这种模型的门限回归的命令要如何编呢

图片内容如下：
TFP = b0 + b1 TFP-1 + b2 in-1 × D(c < τ) + b3 in-1×D(c ≥ τ) + b4 im-1 × D(c < τ) + b5 im-1 × D(c ≥ τ)+b6 in-1 × system-1 + b7 im-1 × system-1+b8 in-1 × human-1 + b9 i -1 × human-1+b10 in-1 × transport-1 + b11 im-1 × transport-1+b12 lnagdp + b13 absorption-1+ ∑ bj ctrl + ε
其中，c 为门限变量，按照前述的理论模型，门限变量主要有 3 个：一个地区的经济发展水平、要素禀赋水平和制度环境，分别采用人均 GDP 、人力资本水平以及市场化水平予以代理。TFP 为全要素生产率，in 代表自主创新投入，im 代表模仿性创新投入，system 为制度变量,human 代表人力资本变量 ，transport 代表交通基础设施 ，agdp 代 表 人 均GDP，取对数形式，ctrl 为其他控制变量， -1 代表滞后一期。

这个不能用 Stata 做（除非你自己另外写程式）！

球大神帮忙看一下
xthreg lncs lninc, rx(lnhw) qx(lnyue) thnum(1) bs(300) trim(0.01) grid(10
> 0)
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.6553        0.6303        0.6643
-----------------------------------------------------

Threshold effect test (bootstrap = 300):
-----------------------------------------------------------------------
Threshold |       RSS        MSE      Fstat    Prob   Crit10    Crit5
-----------+-----------------------------------------------------------
    Single |    2.2006     0.0050      29.97  0.0033  17.6358  21.5110
-----------------------------------------------------------------------

--------------------
Threshold |   Crit1
-----------+--------
    Single | 25.4068
--------------------

Fixed-effects (within) regression               Number of obs      =      
> 455
Group variable: city                            Number of groups   =      
>  35

R-sq:  within  = 0.9420                         Obs per group: min =      
>  13
       between = 0.6695                                        avg =      1
> 3.0
       overall = 0.8208                                        max =      
>  13

                                                F(3,417)           =   2259
> .27
corr(u_i, Xb)  = -0.4265                        Prob > F           =    0.0
> 000

---------------------------------------------------------------------------
> ---
        lncs |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interv
> al]
-------------+-------------------------------------------------------------
> ---
       lninc |   .6587557   .0282209    23.34   0.000     .6032827    .7142
> 287
             |
_cat#c.lnhw |
          0  |   .1990422   .0249937     7.96   0.000     .1499129    .2481
> 714
          1  |   .2052844   .0252257     8.14   0.000     .1556989    .2548
> 699
             |
       _cons |   1.295394    .118624    10.92   0.000     1.062218    1.528
> 569
-------------+-------------------------------------------------------------
> ---
     sigma_u |  .14364899
     sigma_e |  .07254307
         rho |  .79679538   (fraction of variance due to u_i)
---------------------------------------------------------------------------
> ---
F test that all u_i=0: F(34, 417) = 35.40                    Prob > F = 0.0
> 000

球大神帮忙看一下
xthreg lncs lninc, rx(lnhw) qx(lnyue) thnum(1) bs(300) trim(0.01) grid(10
> 0)
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.6553        0.6303        0.6643
-----------------------------------------------------

Threshold effect test (bootstrap = 300):
-----------------------------------------------------------------------
Threshold |       RSS        MSE      Fstat    Prob   Crit10    Crit5
-----------+-----------------------------------------------------------
    Single |    2.2006     0.0050      29.97  0.0033  17.6358  21.5110
-----------------------------------------------------------------------

--------------------
Threshold |   Crit1
-----------+--------
    Single | 25.4068
--------------------

Fixed-effects (within) regression               Number of obs      =      
> 455
Group variable: city                            Number of groups   =      
>  35

R-sq:  within  = 0.9420                         Obs per group: min =      
>  13
       between = 0.6695                                        avg =      1
> 3.0
       overall = 0.8208                                        max =      
>  13

                                                F(3,417)           =   2259
> .27
corr(u_i, Xb)  = -0.4265                        Prob > F           =    0.0
> 000

---------------------------------------------------------------------------
> ---
        lncs |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interv
> al]
-------------+-------------------------------------------------------------
> ---
       lninc |   .6587557   .0282209    23.34   0.000     .6032827    .7142
> 287
             |
_cat#c.lnhw |
          0  |   .1990422   .0249937     7.96   0.000     .1499129    .2481
> 714
          1  |   .2052844   .0252257     8.14   0.000     .1556989    .2548
> 699
             |
       _cons |   1.295394    .118624    10.92   0.000     1.062218    1.528
> 569
-------------+-------------------------------------------------------------
> ---
     sigma_u |  .14364899
     sigma_e |  .07254307
         rho |  .79679538   (fraction of variance due to u_i)
---------------------------------------------------------------------------
> ---
F test that all u_i=0: F(34, 417) = 35.40                    Prob > F = 0.0
> 000

球大神帮忙看一下
xthreg lncs lninc, rx(lnhw) qx(lnyue) thnum(1) bs(300) trim(0.01) grid(10
> 0)
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.6553        0.6303        0.6643
-----------------------------------------------------

Threshold effect test (bootstrap = 300):
-----------------------------------------------------------------------
Threshold |       RSS        MSE      Fstat    Prob   Crit10    Crit5
-----------+-----------------------------------------------------------
    Single |    2.2006     0.0050      29.97  0.0033  17.6358  21.5110
-----------------------------------------------------------------------

--------------------
Threshold |   Crit1
-----------+--------
    Single | 25.4068
--------------------

Fixed-effects (within) regression               Number of obs      =      
> 455
Group variable: city                            Number of groups   =      
>  35

R-sq:  within  = 0.9420                         Obs per group: min =      
>  13
       between = 0.6695                                        avg =      1
> 3.0
       overall = 0.8208                                        max =      
>  13

                                                F(3,417)           =   2259
> .27
corr(u_i, Xb)  = -0.4265                        Prob > F           =    0.0
> 000

---------------------------------------------------------------------------
> ---
        lncs |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interv
> al]
-------------+-------------------------------------------------------------
> ---
       lninc |   .6587557   .0282209    23.34   0.000     .6032827    .7142
> 287
             |
_cat#c.lnhw |
          0  |   .1990422   .0249937     7.96   0.000     .1499129    .2481
> 714
          1  |   .2052844   .0252257     8.14   0.000     .1556989    .2548
> 699
             |
       _cons |   1.295394    .118624    10.92   0.000     1.062218    1.528
> 569
-------------+-------------------------------------------------------------
> ---
     sigma_u |  .14364899
     sigma_e |  .07254307
         rho |  .79679538   (fraction of variance due to u_i)
---------------------------------------------------------------------------
> ---
F test that all u_i=0: F(34, 417) = 35.40                    Prob > F = 0.0
> 000