R语言与Stata结果不一致

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

直接上代码

1. library(car)
   
2. library(plm)
   
3. library(estimatr)
   
4. Paneldata<-pdata.frame(Data,index=c("city","year"))
   
5. 
   
6. ols1=lm_robust(lce~letp+lur,data=Paneldata)
   
7. ols2=lm_robust(lce~letp+lur+ecs,data=Paneldata)
   
8. ols3=lm_robust(lce~letp+lur+lecs,data=Paneldata)
   
9. 
   
10. summary(ols1) #结果和stata相同
    
11. summary(ols2) #加入非对数项后fdi in2 ecs后和stata出现差异
    
12. summary(ols3) #加上对数后，结果与stata相同

复制代码

R的回归结果

lm_robust(formula = lce ~ letp + lur + ecs, data = Paneldata)

Standard error type:  HC2

Coefficients:
            Estimate Std. Error t value   Pr(>|t|) CI Lower CI Upper   DF
(Intercept) 4.054283   0.124363  32.600 9.140e-194 3.810415 4.298151 2435
letp        0.377539   0.012566  30.045 5.743e-169 0.352898 0.402180 2435
lur         0.249592   0.018549  13.456  7.261e-40 0.213219 0.285965 2435
ecs         0.005104   0.001295   3.942  8.309e-05 0.002565 0.007642 2435

Multiple R-squared:  0.5817 ,        Adjusted R-squared:  0.5812
F-statistic:  1348 on 3 and 2435 DF,  p-value: < 2.2e-16

stata的回归结果

. reg lce letp lur ecs, r

Linear regression                               Number of obs     =      2,439
                                                F(3, 2435)        =    1379.17
                                                Prob > F          =     0.0000
                                                R-squared         =     0.5811
                                                Root MSE          =     .73403

------------------------------------------------------------------------------
             |               Robust
         lce |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        letp |   .3830517   .0125671    30.48   0.000     .3584084     .407695
        lur  |   .2470695   .0187324    13.19   0.000     .2103364    .2838026
         ecs |  -.0044635   .0011671    -3.82   0.000    -.0067521   -.0021749
       _cons |   4.859845   .1195206    40.66   0.000     4.625473    5.094218
------------------------------------------------------------------------------

其中ecs的正负都不一样，但在取对数后lecs的结果与Stata相同请问是什么问题导致？

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Stata tata R语言 regression regressio