悬赏 20 个论坛币 未解决
求助门限回归命令xtptm的命令详解与结果详解,以及比较好的测试数据
一直没有人帮我解释,我自己设置了一个微型的测试数据在excel中的sheet2里面。
数据共含6个个体,29期的观测值。变量共两个,即自变量x与因变量y。
当x<=26时,y=0.8x+c1;当x>26,y=-0.6x+c2。
我用命令xtptm y,rx(x) thrvar(x) iters(3000) trim(0.05),得到结果如下:
================ Fundamental information: =================== 
Number of Regime independent variables:    0
Number of Regime dependent variables:    1
Number of individuals in panel:    6
Number of periods in panel:   29
================ Ordinary Fixed effect regression: ========= 
Sum of Squared Residuls:     739.0820
Stardard error of regression:       2.1037
Regression Result(ordinary standard error):
----------------------------------------------------
    |       Coef         Std           t        Prob
----+-----------------------------------------------
  1 |     0.4216      0.0233     18.0876      0.0000
----------------------------------------------------
Regression Result(Robust standard error):
-----------------------------------------------------
    |       Coef   Std_Robust           t        Prob
----+------------------------------------------------
  1 |     0.4216       0.0273     15.4616      0.0000
-----------------------------------------------------
================ Single threshold regession: =================== 
Minimized Sum of Squared Residuals:  169.4322
Standard error of residuals:    1.0103
Threshold estimator:   29.0000
95% conf. intv. of threshold:   29.0000   29.0000
LR Critical value to test gamma=gamma0:    7.3523
Threshold regression(Ordinary Std. Error): 
----------------------------------------------------
    |       Coef         Std           t        prob
----+-----------------------------------------------
  1 |     0.7012      0.0163     43.0440      0.0000
  2 |     0.5004      0.0117     42.8418      0.0000
----------------------------------------------------
Threshold regression(Robust Std. Error): 
-----------------------------------------------------
    |       Coef   Std_Robust           t        prob
----+------------------------------------------------
  1 |     0.7012       0.0149     47.1453      0.0000
  2 |     0.5004       0.0116     43.1624      0.0000
-----------------------------------------------------
Note: Critcal (Inverse CDF of LR stat): -2*ln(1-sqrt(1-alpha))
Ho: No threshold; Ha: Single threshold
Number of bootstrap:3000
F-stat & Prob:  558.1105   0.0000
F-critical value of 90% 95% 99%: 
                 1
    +---------------+
  1 |  2.802225033  |
  2 |  3.890695915  |
  3 |  6.932982111  |
    +---------------+
Thresholds in single model:
  29
============== Descrpitive statistic in each regime: =========== 
Descriptive statistics of y-X-THR at regime :  1
----------------------------------------------------------------
    |       Mean         Std         Min         Max       Count
----+-----------------------------------------------------------
  1 |    30.3895      7.6025     14.0369     45.2109         134
  2 |    19.8582      5.4458     10.0000     28.0000         134
  3 |    19.8582      5.4458     10.0000     28.0000         134
----------------------------------------------------------------
Descriptive statistics of y-X-THR at regime :  2
----------------------------------------------------------------
    |       Mean         Std         Min         Max       Count
----+-----------------------------------------------------------
  1 |    32.4134      6.4932     21.8680     42.6598          40
  2 |    31.5000      1.5359     29.0000     34.0000          40
  3 |    31.5000      1.5359     29.0000     34.0000          40
----------------------------------------------------------------
LR and Thresholds series are stored in Stata matrix: LR#
You can observe the scatter plot by: _matplot LR#, columns(1 2)
----------------------------------------------------------------------------------------
我最不能明白的为什么Threshold值是29,不是26或者27。
另外其他的结果我也不太明白,数据在附件中,我的命令和参数的使用,以及测试数据做的有什么不对的地方也请高手指点。
~~~~~~~~~~~~~~~~~~帮同学问的~~~~~~~~~~~~~~~~~~~~~~~~~~