全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 统计软件培训班VIP答疑区
2522 3
2011-12-26
连老师,我想问一下门槛回归命令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。
另外其他的结果我也不太明白,我的命令和参数的使用,以及测试数据做的有什么不对的地方还请连老师指点。



附件: 您需要登录才可以下载或查看附件。没有帐号?我要注册
二维码

扫码加我 拉你入群

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

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

全部回复
2011-12-27 11:04:05
你可以使用一个更大的样本量来执行测试;另外,测试过程中,尽量让两组的系数差别大一点,比如一组的系数是 0.8,另一组的是 -0.6,这样才比较容易看出问题。

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

扫码加我 拉你入群

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

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

2011-12-27 13:10:46
arlionn 发表于 2011-12-27 11:04
你可以使用一个更大的样本量来执行测试;另外,测试过程中,尽量让两组的系数差别大一点,比如一组的系数是 ...
感谢连老师的帮助,我确实设置的是:比如一组的系数是 0.8,另一组的是 -0.6。当然我的测试样本是小了点。我会再试试看大样本的。

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

扫码加我 拉你入群

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

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

2011-12-31 22:08:40
发邮件给我,用我写的程序试一下吧。
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群