Notice: On March 31, it was **announced** that Statalist is moving from an email list to a **forum**. The old list will shut down on April 23, and its replacement, **statalist.org** is already up and running.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
"Scholes, Shaun" <s.scholes@ucl.ac.uk> |

To |
"statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu> |

Subject |
st: RE: How to calculate standardized difference in means with survey weighted data? |

Date |
Mon, 5 Mar 2012 10:49:43 +0000 |

This was a lot trickier than I thought...and I cannot answer with 100% confidence (or maybe not even 95%). But I found a posting called "RE: st: svy: mean and descriptive tables " from July 2007. And I was able to replicate the following results using estat sd (using that posting: and I am using Stata SE v12.0): clear webuse nhanes2 svyset psu [pw=finalwgt], strata(strata) svy: mean smsa height weight bpsystol bpdiast tcresult tgresult estat sd // e(N) is in a scalar // unweighted number of observations di sqrt(e(N) * el(e(V_srs), 1, 1)) di sqrt(e(N) * el(e(V_srs), 2, 2)) ........ di sqrt(e(N) * el(e(V_srs), 7, 7)) So the quantities you may (I stress may) need are in the scalar e(N) and the diagonal of the variance-covariance matrix (which is where I normally obtain the Standard Errors from). Then I tried the same in the context of using -svy- to estimate mean BMI over two groups in this simple example: clear input id bmi group 1 25.1 1 2 26.2 2 3 26.4 2 4 25.4 1 5 25.7 2 6 25.3 1 7 26.3 2 8 26.3 1 9 25.9 2 10 24.1 2 11 25.6 1 12 25.4 1 13 28.3 2 14 27.2 1 15 26.2 2 16 26.2 1 17 28.2 1 18 27.1 2 19 28.1 1 20 28.4 2 end gen weight=1 svyset [pweight=weight] svy:mean bmi,over(group) estat sd mat A = e(_N) di sqrt(A[1,1]* el(e(V_srs), 1, 1)) di sqrt(A[1,2]* el(e(V_srs), 2, 2)) If this is on the right lines then you could adapt the posting from July 2007 and put all these into tables using the programmes in the July 2007 posting? Hope this helps Best wishes Shaun -----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Lok Wong Sent: 05 March 2012 00:46 To: statalist@hsphsun2.harvard.edu Cc: LOK WONG Subject: st: How to calculate standardized difference in means with survey weighted data? I need to calculate the standardized bias (the difference in means divided by the pooled standard deviation) with survey weighted data using STATA. I am comparing the means of 2 groups (Y: treatment and control) for a list of X predictor variables. The purpose is to evaluate differences before and after propensity score weighting (not matching so I cannot use PSMATCH2 or other similar packages). This is as far I got: svy: mean X, over (Y) estat sd lincom [X]1 - lincom [x]0 I calculated the means by treatment/control groups. Then obtained the standard deviations for each means (as the SE is reported by svy: means) I used lincom to obtain the difference in means from the svy post-estimation results. How do I now extract the stored standard deviations for the 2 means, so I can divide the difference in means by the pooled standard deviation? I need to do this twice for 20 variables, so I don't want to just read the output results and calculate by hand. Any suggestions? I did see an earlier posting (from 2005) on Standardized Response Mean, but the suggested code (diff in change score / sd of change score) does not address how to use survey weighted data. Thanks! Lok Wong Samson Doctoral Candidate JHSPH lokwongedu@gmail.com * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

**Follow-Ups**:**Re: st: RE: How to calculate standardized difference in means with survey weighted data?***From:*Steve Samuels <sjsamuels@gmail.com>

**Re: st: RE: How to calculate standardized difference in means with survey weighted data?***From:*Steve Samuels <sjsamuels@gmail.com>

**References**:**st: How to calculate standardized difference in means with survey weighted data?***From:*Lok Wong <lokwongedu@gmail.com>

- Prev by Date:
**st: Cannot drop a program** - Next by Date:
**Re: st: Cannot drop a program** - Previous by thread:
**st: How to calculate standardized difference in means with survey weighted data?** - Next by thread:
**Re: st: RE: How to calculate standardized difference in means with survey weighted data?** - Index(es):