请教连老师xtptm的问题，或者应该怎么处理

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

连老师，我想问一下门槛回归命令xtptm的命令详解与结果详解，我自己设置了一个微型的测试数据在附件中
数据设置如下：
数据共含6个个体，29期的观测值。变量共两个，即自变量x与因变量y，x也是门槛变量。
当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。
另外其他的结果我也不太明白，我的命令和参数的使用，以及测试数据做的有什么不对的地方还请连老师指点。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：xtptm 怎么处理 xtp PTM 连老师 测试 periods excel 因变量

你可以使用一个更大的样本量来执行测试；另外，测试过程中，尽量让两组的系数差别大一点，比如一组的系数是 0.8，另一组的是 -0.6，这样才比较容易看出问题。

至于其它结果的解释，需要仔细看 ado 文件才能明白，在这个文件中，王群勇老师有一段文字说明。
其他的细节问题还需要咨询程序的写作者——南开大学的王群勇老师。

arlionn 发表于 2011-12-27 11:04
你可以使用一个更大的样本量来执行测试；另外，测试过程中，尽量让两组的系数差别大一点，比如一组的系数是 ...

感谢连老师的帮助，我确实设置的是：比如一组的系数是 0.8，另一组的是 -0.6。当然我的测试样本是小了点。我会再试试看大样本的。

不知道门槛回归还有什么别的方法吗？虽然我编程还行，但数学理论太差，所以还是自己写不了程序。
还望连老师指点。

发邮件给我，用我写的程序试一下吧。