 |
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
Re: st: Multiple mean comparisons with complex survey design
At 02:28 PM 2/19/2008, Sven Klingemann wrote:
Hi,
I would like to compare means across six groups and adjust the significance
test accordingly. I noticed that the "mtest" command (and associated
Bonferroni option) is not allowed after the svy:mean command. Is there any
recommended alternative for multiple group comparisons using a complex
survey design?
Thank you!
At the risk of being accused of plagiarizing someone else's words -
Yes You Can :). Example:
. webuse nhanes2f, clear
. svy: mean health, over(race)
(running mean on estimation sample)
Survey: Mean estimation
Number of strata = 31 Number of obs = 10335
Number of PSUs = 62 Population size = 1.2e+08
Design df = 31
White: race = White
Black: race = Black
Other: race = Other
--------------------------------------------------------------
| Linearized
Over | Mean Std. Err. [95% Conf. Interval]
-------------+------------------------------------------------
health |
White | 3.663833 .0252884 3.612257 3.715409
Black | 3.117221 .0379095 3.039904 3.194538
Other | 3.510359 .080302 3.346582 3.674136
--------------------------------------------------------------
. test [#1]White = [#1]Black = [#1]Other, mtest(b)
Adjusted Wald test
( 1) [health]White - [health]Black = 0
( 2) [health]White - [health]Other = 0
---------------------------------------
| F(df,30) df p
-------+-------------------------------
(1) | 141.48 1 0.0000 #
(2) | 3.65 1 0.1304 #
-------+-------------------------------
all | 69.05 2 0.0000
---------------------------------------
# Bonferroni adjusted p-values
-------------------------------------------
Richard Williams, Notre Dame Dept of Sociology
OFFICE: (574)631-6668, (574)631-6463
HOME: (574)289-5227
EMAIL: [email protected]
WWW: http://www.nd.edu/~rwilliam
*
* 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/