全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
3853 3
2011-07-06
悬赏 2 个论坛币 已解决


问:ttest ZBMI_10==ZBMI_09 if eatstuchange==1,这是ttest的命令
(意思是eatstuchange变量的取值由0变成1时,10年的ZBMI09年的ZBMI有没有显著的差异)。
我要在此基础上加入一些控制变量,如性别,能用其它命令来实现么,比如reg,请教如何实现与这个ttest的结果一致?


eatstuchange的定义为:
g eatstuchange=.
replace eatstuchange=0 ifeatstu_09==0&eatstu_10==0
replace eatstuchange=1 ifeatstu_09==0&eatstu_10==1)


我的数据如下:
  

stucode

  
  

eatstuchange

  
  

ZBMI_10

  
  

ZBMI_09

  
  

331061118

  
  

1

  
  

-1.843575

  
  

-2.750096

  
  

331051101

  
  

1

  
  

-0.7376546

  
  

-2.688892

  
  

331031116

  
  

0

  
  

-1.878961

  
  

-2.667522

  
  

331021212

  
  

1

  
  

-1.848347

  
  

-2.523079

  
  

325041118

  
  

0

  
  

-1.917948

  
  

-2.497282

  
  

331041118

  
  

1

  
  

-2.216679

  
  

-2.220008

  
  

323021117

  
  

0

  
  

-2.651529

  
  

-2.019347

  
  

331081132

  
  

1

  
  

-1.427066

  
  

-1.977728

  
  

332041137

  
  

0

  
  

-2.151803

  
  

-1.954701

  
  

326071134

  
  

1

  
  

0.372839

  
  

-1.873647

  
  

325041210

  
  

0

  
  

-2.045145

  
  

-1.860512

  
  

331031102

  
  

0

  
  

-1.517159

  
  

-1.814215

  
  

326011118

  
  

0

  
  

-1.355653

  
  

-1.791252

  
  

326061108

  
  

0

  
  

-0.8850324

  
  

-1.752629

  
  

325032213

  
  

1

  
  

-1.286155

  
  

-1.747417

  
  

331051109

  
  

1

  
  

-2.092448

  
  

-1.738412

  
  

327041227

  
  

1

  
  

-0.7699418

  
  

-1.733157

  
  

322031213

  
  

1

  
  

-1.04974

  
  

-1.720897

  
  

332061117

  
  

1

  
  

-0.8741232

  
  

-1.687178

  
  

332071108

  
  

0

  
  

-1.935525

  
  

-1.663775

  



ttest ZBMI_10==ZBMI_09 if eatstuchange==1后,结果如附件。

不知道我的问题有没有表述清楚,求解答!(数据也在附件中)
截图00.jpg

原图尺寸 45.27 KB

截图00.jpg

数据.xlsx

大小:9.47 KB

 马上下载

最佳答案

ermutuxia 查看完整内容

关于两个变量均值相等的检验是不能用回归分析得到的,需要用ttest命令
二维码

扫码加我 拉你入群

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

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

全部回复
2011-7-6 17:31:35
关于两个变量均值相等的检验是不能用回归分析得到的,需要用ttest命令
二维码

扫码加我 拉你入群

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

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

2014-10-10 21:17:24
clear

input ///
double stucode        eatstuchange        double ZBMI_10        double ZBMI_09
331061118        1        -1.843575        -2.750096
331051101        1        -0.7376546        -2.688892
331031116        0        -1.878961        -2.667522
331021212        1        -1.848347        -2.523079
325041118        0        -1.917948        -2.497282
331041118        1        -2.216679        -2.220008
323021117        0        -2.651529        -2.019347
331081132        1        -1.427066        -1.977728
332041137        0        -2.151803        -1.954701
326071134        1        0.372839        -1.873647
325041210        0        -2.045145        -1.860512
331031102        0        -1.517159        -1.814215
326011118        0        -1.355653        -1.791252
326061108        0        -0.8850324        -1.752629
325032213        1        -1.286155        -1.747417
331051109        1        -2.092448        -1.738412
327041227        1        -0.7699418        -1.733157
322031213        1        -1.04974        -1.720897
332061117        1        -0.8741232        -1.687178
332071108        0        -1.935525        -1.663775
end
format stucode        %10.0f
sort stucode        eatstuchange
reshape long ZBMI, i(stucode) j(group) str
replace group=substr(group,2,3)
destring group,replace


bysort eatstuchange: ttest ZBMI,by(group)
bysort eatstuchange: reg ZBMI group
二维码

扫码加我 拉你入群

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

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

2014-10-10 21:19:42
. clear

.
. input ///
> double stucode  eatstuchange    double ZBMI_10  double ZBMI_09

        stucode  eatstuc~e     ZBMI_10     ZBMI_09
  1. 331061118       1       -1.843575       -2.750096
  2. 331051101       1       -0.7376546      -2.688892
  3. 331031116       0       -1.878961       -2.667522
  4. 331021212       1       -1.848347       -2.523079
  5. 325041118       0       -1.917948       -2.497282
  6. 331041118       1       -2.216679       -2.220008
  7. 323021117       0       -2.651529       -2.019347
  8. 331081132       1       -1.427066       -1.977728
  9. 332041137       0       -2.151803       -1.954701
10. 326071134       1       0.372839        -1.873647
11. 325041210       0       -2.045145       -1.860512
12. 331031102       0       -1.517159       -1.814215
13. 326011118       0       -1.355653       -1.791252
14. 326061108       0       -0.8850324      -1.752629
15. 325032213       1       -1.286155       -1.747417
16. 331051109       1       -2.092448       -1.738412
17. 327041227       1       -0.7699418      -1.733157
18. 322031213       1       -1.04974        -1.720897
19. 332061117       1       -0.8741232      -1.687178
20. 332071108       0       -1.935525       -1.663775
21. end

. format stucode  %10.0f

. sort stucode    eatstuchange

. reshape long ZBMI, i(stucode) j(group) str
(note: j = _09 _10)

Data                               wide   ->   long
-----------------------------------------------------------------------------
Number of obs.                       20   ->      40
Number of variables                   4   ->       4
j variable (2 values)                     ->   group
xij variables:
                        ZBMI_09 ZBMI_10   ->   ZBMI
-----------------------------------------------------------------------------

. replace group=substr(group,2,3)
(40 real changes made)

. destring group,replace
group has all characters numeric; replaced as byte

.
.
. bysort eatstuchange: ttest ZBMI,by(group)

-------------------------------------------------------------------------------------------------------------
-> eatstuchange = 0

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       9 |       9   -2.002359    .1158926    .3476779   -2.269608   -1.735111
      10 |       9   -1.815417     .169146     .507438   -2.205469   -1.425366
---------+--------------------------------------------------------------------
combined |      18   -1.908888      .10201    .4327918   -2.124111   -1.693666
---------+--------------------------------------------------------------------
    diff |           -.1869422    .2050402               -.6216079    .2477236
------------------------------------------------------------------------------
    diff = mean(9) - mean(10)                                     t =  -0.9117
Ho: diff = 0                                     degrees of freedom =       16

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
Pr(T < t) = 0.1877         Pr(|T| > |t|) = 0.3754          Pr(T > t) = 0.8123

-------------------------------------------------------------------------------------------------------------
-> eatstuchange = 1

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       9 |      11   -2.060046    .1248715     .414152   -2.338278   -1.781815
      10 |      11   -1.252081    .2275381    .7546587   -1.759068   -.7450944
---------+--------------------------------------------------------------------
combined |      22   -1.656064    .1543087    .7237719   -1.976966   -1.335161
---------+--------------------------------------------------------------------
    diff |           -.8079655    .2595506               -1.349379   -.2665525
------------------------------------------------------------------------------
    diff = mean(9) - mean(10)                                     t =  -3.1129
Ho: diff = 0                                     degrees of freedom =       20

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
Pr(T < t) = 0.0027         Pr(|T| > |t|) = 0.0055          Pr(T > t) = 0.9973

. bysort eatstuchange: reg ZBMI group

-------------------------------------------------------------------------------------------------------------
-> eatstuchange = 0

      Source |       SS       df       MS              Number of obs =      18
-------------+------------------------------           F(  1,    16) =    0.83
       Model |    .1572632     1    .1572632           Prob > F      =  0.3754
    Residual |  3.02698599    16  .189186625           R-squared     =  0.0494
-------------+------------------------------           Adj R-squared = -0.0100
       Total |  3.18424919    17  .187308776           Root MSE      =  .43496

------------------------------------------------------------------------------
        ZBMI |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
       group |   .1869422   .2050402     0.91   0.375    -.2477236    .6216079
       _cons |  -3.684839   1.950578    -1.89   0.077    -7.819879    .4502009
------------------------------------------------------------------------------

-------------------------------------------------------------------------------------------------------------
-> eatstuchange = 1

      Source |       SS       df       MS              Number of obs =      22
-------------+------------------------------           F(  1,    20) =    9.69
       Model |  3.59044529     1  3.59044529           Prob > F      =  0.0055
    Residual |  7.41031566    20  .370515783           R-squared     =  0.3264
-------------+------------------------------           Adj R-squared =  0.2927
       Total |  11.0007609    21  .523845759           Root MSE      =   .6087

------------------------------------------------------------------------------
        ZBMI |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
       group |   .8079655   .2595506    3.11   0.005     .2665525    1.349379
       _cons |  -9.331736   2.469143    -3.78   0.001    -14.48228   -4.181193
------------------------------------------------------------------------------
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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