Statalist The Stata Listserver


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: Re: geometric mean with svy (version 9.0)


From   "Ahmed Arif" <arifaa@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   st: Re: geometric mean with svy (version 9.0)
Date   Sun, 19 Mar 2006 16:58:17 -0600

I was able to figure out the answer. after running svy:mean on log
transformed variables, I typed ereturn display,eform(gm) and it gives
me geometric mean and sd. I tried using foreach syntax to automate the
process but was not able to run the query.

foreach x of varlist lnbenz lnchol {
      quietly svy:mean`x'
      gen gm`x'=exp(r(mean))
}

The problem here was r(mean) . Since I am running estimation i wasnt
able to find any equivalence of r(mean). Here is list of scalar I got
when I typed ereturn list

. svy:mean  lnbenz
(running mean on estimation sample)

Survey: Mean estimation

Number of strata =      13          Number of obs    =     647
Number of PSUs   =      27          Population size  = 1.5e+08
                                    Design df        =      14

--------------------------------------------------------------
             |             Linearized
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
      lnbenz |   .2441459    .221858     -.2316921    .7199839
--------------------------------------------------------------

. ereturn list

scalars:
              e(df_r) =  14
         e(singleton) =  0
            e(census) =  0
             e(N_pop) =  145459402.6002451
             e(N_psu) =  27
          e(N_strata) =  13
            e(N_over) =  1
                 e(N) =  647
     e(N_strata_omit) =  0
            e(stages) =  1
              e(k_eq) =  1
           e(k_eform) =  0

macros:
               e(cmd) : "mean"
            e(prefix) : "svy"
           e(cmdname) : "mean"
           e(command) : "mean lnbenz,"
             e(title) : "Survey: Mean estimation"
           e(vcetype) : "Linearized"
         e(estat_cmd) : "svy_estat"
           e(varlist) : "lnbenz"
           e(predict) : "_no_predict"
             e(wtype) : "pweight"
              e(wvar) : "wtsvoc2y"
              e(wexp) : "= wtsvoc2y"
               e(vce) : "linearized"
               e(su1) : "sdmvpsu"
           e(strata1) : "sdmvstra"
        e(properties) : "b V"

matrices:
                 e(b) :  1 x 1
                 e(V) :  1 x 1
           e(_N_subp) :  1 x 1
             e(V_srs) :  1 x 1
                e(_N) :  1 x 1
             e(error) :  1 x 1

functions:
            e(sample)

thanks for the help
Ahmed.


On 3/17/06, Ahmed Arif <arifaa@gmail.com> wrote:
> Dear Stata Guru's,
>
> I want to compute geometric mean and SD from survey data. Here is an example:
>
> . svy:mean  lnbenz,over (sex)
> (running mean on estimation sample)
>
> Survey: Mean estimation
>
> Number of strata =      13          Number of obs    =     647 Number
> of PSUs   =      27          Population size  = 1.5e+08
>                                    Design df        =      14
>
>    _subpop_1: sex = 0. Male
>    _subpop_2: sex = 1. Female
>
> --------------------------------------------------------------
>             |             Linearized
>        Over |       Mean   Std. Err.     [95% Conf. Interval]
> -------------+------------------------------------------------
> lnbenz       |
>   _subpop_1 |   .3884244   .2993531     -.2536241    1.030473
>   _subpop_2 |   .1070038   .1907388     -.3020902    .5160979
> --------------------------------------------------------------
>
> I can manually type
> Di exp(0.3884244)
>
> For each of the category but I have to type hundreds of time since I
> have several variables. I tried egenmore option but I couldn't find
> how to use it with svy command. Can any help point me to the resource
> to find more info? Is there program I can use to automate this
> process? Ideally I would like Stata to run the survey command and then
> use the mean value to compute geometric mean and SD. Any help will be
> greatly appreciated
>
> Thank you
> Ahmed
> Ahmed A. Arif, MD., PhD
> Assistant Professor & Director
> Health Services Research Program
> Department of Family and Community Medicine Texas Tech University
> Health Sciences Center
> 3601 4th St, STOP 8161
> Lubbock, TX 79430
> Tel: 806-743-3048
> Fax: 806-743-1292
> email: ahmed.arif@ttuhsc.edu
>

*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2020 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index