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。
对于双重或者三重门槛的测试,由于相对复杂一点,我没有做,所以也不知道那里面的变量是如何排序的
而且我的论文的数据也用不到三重门槛。