请教：如何在stata中用reg实现与ttest结果一致(两期数据)

问：ttest ZBMI_10==ZBMI_09 if eatstuchange==1，这是ttest的命令
（意思是eatstuchange变量的取值由0变成1时，10年的ZBMI与09年的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后，结果如附件。

不知道我的问题有没有表述清楚，求解答！(数据也在附件中)

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

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

. 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
------------------------------------------------------------------------------