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 |
Steve Samuels <sjsamuels@gmail.com> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re: st: Multistage sampling svyset |

Date |
Thu, 10 Jun 2010 16:31:49 -0400 |

You still haven't answered my question about the design, James: exactly how were HH selected? To answer your question: the sampling fraction within HH is not ignored; it should be reflected in the sampling weight. (Check to be sure that it is.) The -fpc- applies only to sampling units for which a contribution to variance can be estimated. As you have only one observation per family, there is no within-family variance component to estimate, hence no fpc. In fact even if several individuals had been selected within each hh, fpcs would usually not be computed, because they make no practical difference in the standard errors and because they would be invalid for many sub-population estimates. Steve On Thu, Jun 10, 2010 at 3:58 PM, <James.Valcour@med.mun.ca> wrote: > Thanks for the more detailed read of my original message ;) > > I tried as you suggested, and it works (thank you!). Does that not ignore the selection probability of the individuals within the household? > > I also tried Michael's suggestion of using the singleunit option. I get the same standard error estimates using that option and the more complex original syntax. > > Regardless, I am now getting the estimates I've been looking for. > > Thank you to you both. > > Cheers, > James > > -----Original Message----- > From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Steve Samuels > Sent: June 10, 2010 4:48 PM > To: statalist@hsphsun2.harvard.edu > Subject: Re: st: Multistage sampling svyset > > OK, I read your email in detail. Let pcomm = (N hh in the stratum)/n > hh in the sample from the stratum) > > _If_ the survey was done as you describe the -svyset- should go: > > svyset houseid [pw=pwgt], strata(pi4) fpc(pcomm) > > The proportion of people selected in the HH (the last stage) affects > only the sampling weight. > > Please verify how the HH were selected Was there a list of HH from > which HH were sampled? Or, was there an intermediate sample of areas > (such as blocks) drawn first, with HH sampled only in the selected > areas? > > Steve > > On Thu, Jun 10, 2010 at 3:04 PM, Steve Samuels <sjsamuels@gmail.com> wrote: >> James, what is "pi4"? What is the variable which identifies the >> community? How did you compute the sampling weight "pwgt"? (show an >> example, as well as code.) >> >> Steve >> >> >> On Thu, Jun 10, 2010 at 2:33 PM, <James.Valcour@med.mun.ca> wrote: >>> I'm having a problem using the survey commands to analyzing a survey >>> that used a complex sampling design. >>> >>> Eleven communities were sampled (these were all the communities in the >>> study area). Within each community households were randomly selected >>> from all available households. Within each household one person was >>> randomly select from all members of the household that were eligible to >>> be in the survey (we had a list of everybody eligible in the house). >>> >>> This is the code I used to set the data: >>> >>> svyset houseid [pw=pwgt], strata(pi4) fpc(p_household) || _n, >>> fpc(p_person) >>> >>> where >>> houseid is an identifier for the house in the community; a mix of the >>> community id number and the house number >>> pi4 is an identified for the community (1-11) >>> p_household is the probability of selection of the household in the >>> community >>> p_person is the probability of selection of the person in the household >>> >>> The problem I'm getting is when I'm trying to estimate the proportion of >>> study with certain ailments. I'm not getting confidence intervals >>> included in the estimates. The message that is given is " Note: missing >>> standard errors because of stratum with single sampling unit." >>> >>> I've tried other ways of specify the svyset command and tried >>> simplifying to the point of just assuming that household were sampled >>> but I get the same error. (for the record here's that command "svyset >>> _n [pw=pwgt], strata(pi4) fpc(p_household)") >>> >>> Any help would be appreciated, >>> >>> Thanks! >>> James Valcour >>> >>> -- >>> Dr. James Valcour >>> Division of Community Health and Humanities >>> Faculty of Medicine >>> The Health Sciences Centre >>> Memorial University of Newfoundland >>> St. John's, NL, Canada >>> A1B 3V6 >>> Tel: 709-777-2237 >>> Fax: 709-777-7382 >>> Email: james.valcour@med.mun.ca >>> >>> >>> >>> * >>> * 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/ >>> >> >> >> >> -- >> Steven Samuels >> sjsamuels@gmail.com >> 18 Cantine's Island >> Saugerties NY 12477 >> USA >> Voice: 845-246-0774 >> Fax: 206-202-4783 >> > > > > -- > Steven Samuels > sjsamuels@gmail.com > 18 Cantine's Island > Saugerties NY 12477 > USA > Voice: 845-246-0774 > Fax: 206-202-4783 > > * > * 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/ > -- Steven Samuels sjsamuels@gmail.com 18 Cantine's Island Saugerties NY 12477 USA Voice: 845-246-0774 Fax: 206-202-4783 * * 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: Multistage sampling svyset***From:*<James.Valcour@med.mun.ca>

**References**:**st: Multistage sampling svyset***From:*<James.Valcour@med.mun.ca>

**Re: st: Multistage sampling svyset***From:*Steve Samuels <sjsamuels@gmail.com>

**Re: st: Multistage sampling svyset***From:*Steve Samuels <sjsamuels@gmail.com>

**RE: st: Multistage sampling svyset***From:*<James.Valcour@med.mun.ca>

- Prev by Date:
**RE: st: Multistage sampling svyset** - Next by Date:
**st: bar graph axis color- frustrated** - Previous by thread:
**RE: st: Multistage sampling svyset** - Next by thread:
**RE: st: Multistage sampling svyset** - Index(es):