问一个有关tobit的问题

clear
set seed 13579
set obs 1000
gen x1=invnorm(uniform())
gen x2=2*invnorm(uniform())
gen e=invnorm(uniform())
gen y=3+x1+3*x2+e

gen z=y if y<10
replace z=10 if y>=10


*drop in 1/100
reg y x1 x2
tobit z x1 x2, ul(10)  //eq1
tobit y x1 x2, ul(10) //eq2

sum y
/*

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
           y |      1000    2.792235    5.905514   -14.7648   24.75783
*/


Q:为什么eq1得到和eq2一样的结果？

因为你两个都用了ul(10)，对y用ul(10)得到的就是z，当然会一样。
如果你这样，就不会有这个问题：

1. . tobit y x1 x2, ul(1000000)
   
2. 
   
3. Tobit regression                                  Number of obs   =       1000
   
4.                                                   LR chi2(2)      =    3547.54
   
5.                                                   Prob > chi2     =     0.0000
   
6. Log likelihood = -1420.5558                       Pseudo R2       =     0.5553
   
7. 
   
8. ------------------------------------------------------------------------------
   
9.            y |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
   
10. -------------+----------------------------------------------------------------
    
11.           x1 |   .9534411   .0315907    30.18   0.000     .8914492    1.015433
    
12.           x2 |    2.99618   .0163966   182.73   0.000     2.964005    3.028356
    
13.        _cons |   3.013828   .0317027    95.07   0.000     2.951617     3.07604
    
14. -------------+----------------------------------------------------------------
    
15.       /sigma |   1.001619   .0223969                      .9576682    1.045569
    
16. ------------------------------------------------------------------------------
    
17.   Obs. summary:          0  left-censored observations
    
18.                       1000     uncensored observations
    
19.                          0 right-censored observations

复制代码

夏目贵志 发表于 2012-9-23 02:44
因为你两个都用了ul(10)，对y用ul(10)得到的就是z，当然会一样。
如果你这样，就不会有这个问题：

感谢夏目贵志:)