. summarize ltthallpfc-rttcofc


    Variable |        Obs        Mean    Std. Dev.       Min        Max

-------------+---------------------------------------------------------

  ltthallpfc |        200      .37333    .0277996       .301        .47

  rtthallpfc |        200     .369105    .0246989       .311       .459

  ltthalmpfc |        200      .39878    .0205248       .349       .459

  rtthalmpfc |        200      .40587    .0211661       .342       .463

   ltthalofc |        200     .367645    .0311044       .274       .521

-------------+---------------------------------------------------------

   rtthalofc |        200      .36987    .0277549       .296       .465

    lttclpfc |        200     .352225    .0207193        .29       .411

    rttclpfc |        200      .36088     .020233       .264       .414

     lttcofc |        200      .30758    .0227173        .25       .365

     rttcofc |        200      .31088    .0186112       .252       .379


. tabulate group


      group |      Freq.     Percent        Cum.

------------+-----------------------------------

          1 |        100       50.00       50.00

          2 |         40       20.00       70.00

          3 |         60       30.00      100.00

------------+-----------------------------------

      Total |        200      100.00


. tabstat ltthallpfc-rttcofc, by(group)

Summary statistics: mean
  by categories of: group 

   group |  lt~llpfc  rt~llpfc  ltt~mpfc  rtt~mpfc  ltth~ofc  rtth~ofc  lttclpfc  rttclpfc   lttcofc   rttcofc
---------+----------------------------------------------------------------------------------------------------
       1 |    .37544    .37071     .3991    .40339    .37523     .3732    .35361    .35944    .31272    .31061
       2 |   .374875    .37045    .39555   .409275      .361     .3693    .35115   .362875    .30515   .315375
       3 |  .3687833  .3655333     .4004  .4077333  .3594333     .3647  .3506333    .36195  .3006333  .3083333
---------+----------------------------------------------------------------------------------------------------
   Total |    .37333   .369105    .39878    .40587   .367645    .36987   .352225    .36088    .30758    .31088
--------------------------------------------------------------------------------------------------------------


. pwcorr ltthallpfc-rttcofc, star(.05)

             | lt~llpfc rt~llpfc ltt~mpfc rtt~mpfc ltth~ofc rtth~ofc lttclpfc rttclpfc lttcofc rttcofc
-------------+----------------------------------------------------------------------------------------- 
  ltthallpfc |   1.0000 
  rtthallpfc |   0.4681*  1.0000 
  ltthalmpfc |   0.5358*  0.3849*  1.0000 
  rtthalmpfc |   0.4018*  0.5690*  0.6051*  1.0000 
   ltthalofc |   0.5068*  0.2856*  0.3209*  0.1968*  1.0000 
   rtthalofc |   0.3311*  0.5691*  0.1631*  0.2788*  0.3984*  1.0000 
    lttclpfc |   0.3995*  0.4462*  0.3287*  0.3409*  0.2640*  0.3323*  1.0000 
    rttclpfc |   0.2455*  0.3640*  0.2570*  0.3855*  0.1937*  0.2991*  0.4855*  1.0000
     lttcofc |   0.2774*  0.3128*  0.2063*  0.2530*  0.2914*  0.3179*  0.5015*  0.3517* . 1.0000
     rttcofc |   0.2053*  0.3370*  0.2149*  0.2690*  0.2937*  0.2905*  0.4205* . 0.6034*  0.4983* . 1.0000
 
. manova ltthallpfc-rttcofc = group

                       Number of obs =        200

                       W = Wilks' lambda      L = Lawley-Hotelling trace
                       P = Pillai's trace     R = Roy's largest root

              Source | Statistic        df    F(df1,     df2) =   F   Prob>F
          -----------+-------------------------------------------------------
               group |W   0.7790         2     20.0    376.0     2.50 0.0004 e
                     |P   0.2314               20.0    378.0     2.47 0.0005 a
                     |L   0.2704               20.0    374.0     2.53 0.0004 a
                     |R   0.2057               10.0    189.0     3.89 0.0001 u
                     |-------------------------------------------------------
            Residual |                 197
          -----------+-------------------------------------------------------
               Total |                 199
          -------------------------------------------------------------------
                       e = exact, a = approximate, u = upper bound on F

. manovatest, showorder

 Order of columns in the design matrix
      1: (group==1)
      2: (group==2)
      3: (group==3)
      4: _cons


. matrix c1=(1, -1, 0, 0)

. manovatest, test(c1)

 Test constraint
 (1)    1.group - 2.group = 0

                       W = Wilks' lambda      L = Lawley-Hotelling trace
                       P = Pillai's trace     R = Roy's largest root

              Source | Statistic        df    F(df1,     df2) =   F   Prob>F
          -----------+-------------------------------------------------------
          manovatest |W   0.8659         1     10.0    188.0     2.91 0.0021 e
                     |P   0.1341               10.0    188.0     2.91 0.0021 e
                     |L   0.1549               10.0    188.0     2.91 0.0021 e
                     |R   0.1549               10.0    188.0     2.91 0.0021 e
                     |-------------------------------------------------------
            Residual |                 197
          -------------------------------------------------------------------
                       e = exact, a = approximate, u = upper bound on F


. matrix c2=(1, 0, -1, 0)

. manovatest, test(c2)

 Test constraint
 (1)    1.group - 3.group = 0

                       W = Wilks' lambda      L = Lawley-Hotelling trace
                       P = Pillai's trace     R = Roy's largest root

              Source | Statistic        df    F(df1,     df2) =   F   Prob>F
          -----------+-------------------------------------------------------
          manovatest |W   0.8608         1     10.0    188.0     3.04 0.0014 e
                     |P   0.1392               10.0    188.0     3.04 0.0014 e
                     |L   0.1617               10.0    188.0     3.04 0.0014 e
                     |R   0.1617               10.0    188.0     3.04 0.0014 e
                     |-------------------------------------------------------
            Residual |                 197
          -------------------------------------------------------------------
                       e = exact, a = approximate, u = upper bound on F

. matrix c3=(0, 1, -1, 0)

. manovatest, test(c3)

 Test constraint
 (1)    2.group - 3.group = 0

                       W = Wilks' lambda      L = Lawley-Hotelling trace
                       P = Pillai's trace     R = Roy's largest root

              Source | Statistic        df    F(df1,     df2) =   F   Prob>F
          -----------+-------------------------------------------------------
          manovatest |W   0.9386         1     10.0    188.0     1.23 0.2747 e
                     |P   0.0614               10.0    188.0     1.23 0.2747 e
                     |L   0.0654               10.0    188.0     1.23 0.2747 e
                     |R   0.0654               10.0    188.0     1.23 0.2747 e
                     |-------------------------------------------------------
            Residual |                 197
          -------------------------------------------------------------------
                       e = exact, a = approximate, u = upper bound on F

. margins group, predict(equation(ltthallpfc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(ltthallpfc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |     .37544   .0027778   135.16   0.000     .3699619    .3809181

          2  |    .374875   .0043921    85.35   0.000     .3662134    .3835366

          3  |   .3687833   .0035861   102.84   0.000     .3617112    .3758555

------------------------------------------------------------------------------


. margins group, predict(equation(rtthallpfc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(rtthallpfc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |     .37071   .0024712   150.01   0.000     .3658367    .3755833

          2  |     .37045   .0039073    94.81   0.000     .3627446    .3781554

          3  |   .3655333   .0031903   114.58   0.000     .3592419    .3718248

------------------------------------------------------------------------------


. margins group, predict(equation(ltthalmpfc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(ltthalmpfc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |      .3991   .0020557   194.15   0.000     .3950461    .4031539

          2  |     .39555   .0032503   121.70   0.000     .3891402    .4019598

          3  |      .4004   .0026538   150.88   0.000     .3951664    .4056336

------------------------------------------------------------------------------


. margins group, predict(equation(rtthalmpfc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(rtthalmpfc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |     .40339   .0021119   191.01   0.000     .3992251    .4075549

          2  |    .409275   .0033392   122.57   0.000     .4026898    .4158602

          3  |   .4077333   .0027265   149.55   0.000     .4023565    .4131101

------------------------------------------------------------------------------


. margins group, predict(equation(ltthalofc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(ltthalofc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |     .37523   .0030308   123.80   0.000     .3692529    .3812071

          2  |       .361   .0047922    75.33   0.000     .3515494    .3704506

          3  |   .3594333   .0039128    91.86   0.000      .351717    .3671497

------------------------------------------------------------------------------


. margins group, predict(equation(rtthalofc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(rtthalofc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |      .3732   .0027646   134.99   0.000     .3677479    .3786521

          2  |      .3693   .0043713    84.48   0.000     .3606795    .3779205

          3  |      .3647   .0035691   102.18   0.000     .3576614    .3717386

------------------------------------------------------------------------------


. margins group, predict(equation(lttclpfc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(lttclpfc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |     .35361   .0020777   170.20   0.000     .3495127    .3577073

          2  |     .35115   .0032851   106.89   0.000     .3446716    .3576284

          3  |   .3506333   .0026823   130.72   0.000     .3453437    .3559229

------------------------------------------------------------------------------


. margins group, predict(equation(rttclpfc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(rttclpfc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |     .35944   .0020281   177.23   0.000     .3554404    .3634396

          2  |    .362875   .0032067   113.16   0.000     .3565511    .3691989

          3  |     .36195   .0026183   138.24   0.000     .3567866    .3671134

------------------------------------------------------------------------------


. margins group, predict(equation(lttcofc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(lttcofc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |     .31272   .0022181   140.98   0.000     .3083457    .3170943

          2  |     .30515   .0035072    87.01   0.000     .2982336    .3120664

          3  |   .3006333   .0028636   104.99   0.000     .2949861    .3062805

------------------------------------------------------------------------------


. margins group, predict(equation(rttcofc))


Adjusted predictions                            Number of obs     =        200


Expression   : Linear prediction, predict(equation(rttcofc))


------------------------------------------------------------------------------

             |            Delta-method

             |     Margin   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          1  |     .31061   .0018541   167.52   0.000     .3069535    .3142665

          2  |    .315375   .0029316   107.58   0.000     .3095936    .3211564

          3  |   .3083333   .0023937   128.81   0.000     .3036128    .3130538

------------------------------------------------------------------------------


. margins, dydx(group) predict(equation(ltthallpfc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(ltthallpfc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |   -.000565   .0051968    -0.11   0.914    -.0108135    .0096835

          3  |  -.0066567   .0045361    -1.47   0.144    -.0156023     .002289

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(rtthallpfc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(rtthallpfc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |    -.00026   .0046231    -0.06   0.955    -.0093772    .0088572

          3  |  -.0051767   .0040354    -1.28   0.201    -.0131348    .0027815

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(ltthalmpfc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(ltthalmpfc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |    -.00355   .0038458    -0.92   0.357    -.0111342    .0040342

          3  |      .0013   .0033569     0.39   0.699      -.00532      .00792

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(rtthalmpfc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(rtthalmpfc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |    .005885    .003951     1.49   0.138    -.0019067    .0136767

          3  |   .0043433   .0034487     1.26   0.209    -.0024579    .0111445

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(ltthalofc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(ltthalofc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |    -.01423   .0056702    -2.51   0.013    -.0254121   -.0030479

          3  |  -.0157967   .0049493    -3.19   0.002    -.0255572   -.0060362

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(rtthalofc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(rtthalofc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |     -.0039   .0051722    -0.75   0.452    -.0140999    .0062999

          3  |     -.0085   .0045146    -1.88   0.061    -.0174032    .0004032

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(lttclpfc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(lttclpfc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |    -.00246    .003887    -0.63   0.528    -.0101254    .0052054

          3  |  -.0029767   .0033928    -0.88   0.381    -.0096676    .0037142

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(rttclpfc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(rttclpfc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |    .003435   .0037942     0.91   0.366    -.0040475    .0109175

          3  |     .00251   .0033119     0.76   0.449    -.0040213    .0090413

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(lttcofc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(lttcofc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |    -.00757   .0041497    -1.82   0.070    -.0157536    .0006136

          3  |  -.0120867   .0036222    -3.34   0.001    -.0192299   -.0049434

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.


. margins, dydx(group) predict(equation(rttcofc))


Conditional marginal effects                    Number of obs     =        200


Expression   : Linear prediction, predict(equation(rttcofc))

dy/dx w.r.t. : 2.group 3.group


------------------------------------------------------------------------------

             |            Delta-method

             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]

-------------+----------------------------------------------------------------

       group |

          2  |    .004765   .0034688     1.37   0.171    -.0020757    .0116057

          3  |  -.0022767   .0030278    -0.75   0.453    -.0082477    .0036944

------------------------------------------------------------------------------

Note: dy/dx for factor levels is the discrete change from the base level.



저작자 표시 비영리 변경 금지
신고

'Statistics' 카테고리의 다른 글

Stata - MANOVA (Multivariate Analysis of Variance)  (0) 2017.11.17
Stata - bar graph with error bar  (0) 2017.11.01
Stata - incidence rate  (0) 2017.10.02

twoway (rcap lci uci sample, horizontal)(bar incidence sample, horizontal barwidth(0.7) ylabel(1(1)5, valuelabel angle(0)) ytitle("") legend(off))




저작자 표시 비영리 변경 금지
신고

'Statistics' 카테고리의 다른 글

Stata - MANOVA (Multivariate Analysis of Variance)  (0) 2017.11.17
Stata - bar graph with error bar  (0) 2017.11.01
Stata - incidence rate  (0) 2017.10.02

1. Incidence Rate  계산 

 

350 person-years에서 2명 발생한 경우


cii means 350 2 , poisson  

 

cii #exposure #events , poisson 

 

 

2. 메타 분석 관련

 

이렇게 구한 2개의 연구자료에 대한 Pooling 산출예제

 

ssc install metan

 

clear all
input str20 Study Rate LCL UCL 
"A" .0057143 .000692 .020642
"B" .0022624 .000274 .0081727
end

 

Rate에 1000곱하여 제시하는 경우 

 

clear all
input str20 Study Rate LCL UCL 
"A" 5.7143 .692 20.642
"B" 2.2624 .274 8.1727
end

 

gen ln_Rate=ln(Rate)
gen ln_LCL=ln(LCL)
gen ln_UCL=ln(UCL)

metan ln_Rate ln_LCL ln_UCL, eform random lcols(Study) 
 

저작자 표시 비영리 변경 금지
신고

'Statistics' 카테고리의 다른 글

Stata - MANOVA (Multivariate Analysis of Variance)  (0) 2017.11.17
Stata - bar graph with error bar  (0) 2017.11.01
Stata - incidence rate  (0) 2017.10.02

+ Recent posts