全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 统计软件培训班VIP答疑区
4068 8
2012-03-23
连老师,您的程序我估计了一下,有几个简单的问题还需要请教
1. 门槛的F检验显著,但95%的置信区间范围很大,即几乎全部在7.35的虚线以下,这种门槛值是否应判为不合理。
2. 估计结果中有x_1,x_2,x_3,这些变量应该是按从小到大的分断顺序列出的吧?比如一门槛情况下,x_1就是小于门槛值的估计结果吧?
3. 用这个程序经常出现如下错误:
file C:\Users\ESTABL~1\AppData\Local\Temp\ST_0100000d.tmp could not be opened
matrix has missing values

4. 假如,经检验只有一个门槛值10,估计的x_1为0.5,x_2为-0.2。
在这种情况下,我如果用xtreg y x if q<=10,fe为什么得不到0.5,或者接近0.5的估计结果?
同样,如果用xtreg y x if q>10,fe为什么也得不到-0.2?

再次感谢连老师,以及您写的这么好的程序。

二维码

扫码加我 拉你入群

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

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

全部回复
2012-3-31 21:42:01
1. 这样的门槛值非常不准确,不可用。
2. 是从小到大排序的,系数的含义参见 Hansen (1999);
3. 可能是因为缺漏值没有处理干净;
4. 你需要仔细看看 Hansen(1999) 的论文,搞清楚原理后自然就理解了。
二维码

扫码加我 拉你入群

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

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

2012-4-7 00:04:44
arlionn 发表于 2012-3-31 21:42
1. 这样的门槛值非常不准确,不可用。
2. 是从小到大排序的,系数的含义参见 Hansen (1999);
3. 可能是因 ...
连老师,感谢你的回答,长时间没有回论坛看,今天看到回复,对于第2个问题,我有自己的看法。
如果按照您的回答,那我觉得可能是你的程序出现了BUG,还请连老师查看。
我自己生成一组测试数据(见附件,stata可拿来直接测试,excel可拿来查看生成公式)
其中Y的生成公式为=IF(q>0.5,a+0.8*x+e,a-0.6*x+e)
按您的说法,估计的门槛值应该为0.5,当小于0.5时,估计系数应该为-0.6,当大于0.5时,估计系数应该为0.8
但是,实际的估计系数如下(命令为xtthres y, th(q) d(x)):
——————————————————————————————————————
第一个门槛估计值:.49774173
Fixed-effects (within) regression               Number of obs      =       750
Group variable: id                              Number of groups   =        30
R-sq:  within  = 0.9983                         Obs per group: min =        25
       between = 0.0039                                        avg =      25.0
       overall = 0.2199                                        max =        25
                                                F(2,718)           = 211246.22
corr(u_i, Xb)  = -0.0130                        Prob > F           =    0.0000
------------------------------------------------------------------------------
           y |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         x_1 |   .8036123   .0027283   294.55   0.000      .798256    .8089686
         x_2 |  -.5952335   .0026473  -224.84   0.000     -.600431   -.5900361
       _cons |    18.5634   .0165484  1121.77   0.000     18.53091    18.59589
-------------+----------------------------------------------------------------
     sigma_u |  8.7985837
     sigma_e |  .19331623
         rho |   .9995175   (fraction of variance due to u_i)
------------------------------------------------------------------------------
F test that all u_i=0:     F(29, 718) = 51779.26             Prob > F = 0.0000

——————————————————————————————————————————————
可以看到门槛值的估计为0.497,非常接近0.5,很完美。
x_1,按你设计的应该是小于门槛值0.5的系数,结果为0.803,而按模拟数据应该为-0.6,
x_2,按你的设计应该是大于门槛值0.5的系数,结果为-0.595,而按模拟数据应该为0.8。
这二者正好相反,而且如果用王群勇老师的xtptm可以得到接近模拟数据的结果。
如果把=IF(q>0.5,a+0.8*x+e,a-0.6*x+e) 中的q>0.5改为小于,那将得到相反的估计系数。
综合,我猜想,可能是你的程序的BUG。
对于双重或者三重门槛的测试,由于相对复杂一点,我没有做,所以也不知道那里面的变量是如何排序的
而且我的论文的数据也用不到三重门槛。
附件: 您需要登录才可以下载或查看附件。没有帐号?我要注册
二维码

扫码加我 拉你入群

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

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

2012-4-13 11:50:07
多谢你的提醒,这个的确是我解释有误。我查验了源程序:
viewsource xtthres.ado

adoedit xtthres.ado
主要问题在于我对结果的解释不够清晰,这个我会想办法修改原程序,让结果看起来更明了。

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

扫码加我 拉你入群

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

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

2012-4-13 12:14:28
arlionn 发表于 2012-4-13 11:50
多谢你的提醒,这个的确是我解释有误。我查验了源程序:
viewsource xtthres.ado
十分感谢连老师,以及您辛苦编写的程序,
这为我完成了毕业论文提供了很大的帮助。

关于估计系数的这点问题,我在论文的结果中做一个小说明,连老师不介意吧。
二维码

扫码加我 拉你入群

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

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

2012-4-13 12:16:25
arlionn 发表于 2012-4-13 11:50
多谢你的提醒,这个的确是我解释有误。我查验了源程序:
viewsource xtthres.ado
啊,还有一个小请求,如果程序修改好了,能再发给我一个最新版的吗,establown@gmail.com
多谢连老师。
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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