[R] Complex Surveys...Specifying Design

From: Zodet, Marc <MZodet_at_ahrq.gov>
Date: Sat 24 Jul 2004 - 02:23:35 EST


I need some guidance from someone who is familiar/has some experience with the survey package.  

The data that I am using is from the Medical Expenditure Panel Survey (www.meps.ahrq.gov <http://www.meps.ahrq.gov/> ). The STRATA and PSU variables are varstr01 and varpsu01 respectively. When I try to specify them with the svydesign function I get an error message. An excerpt of my session is as follows...  

library(foreign)

> h60 <- read.xport("h:\\meps\\temptransports\\h60.tpt")

> names(h60) <- casefold(names(h60))

> names(h60)

[1] "dupersid" "sex" "inscov01" "totexp01" "perwt01f" "varstr01" "varpsu01"

> library(survey)

> meps.design <- svydesign(id=~varpsu01, strata=~varstr01, weight=~perwt01f,
data=h60)

Error in svydesign(id = ~varpsu01, strata = ~varstr01, weight = ~perwt01f, :

        Clusters not nested in strata

>
 

I'm not sure what is causing this error message or how to get around it (i.e., remedy the problem). Any insight would be much appreciated.  

Thanks.  

Marc    

Marc W. Zodet, MS

Senior Health Statistician  

Agency for Healthcare Research and Quality

Center for Financing, Access, and Cost Trends

Division of Statistical Research and Methods

540 Gaither Road, Room 5058

Rockville, Maryland 20850  

Phone: 301-427-1563

FAX: 301-427-1276 E-mail: <mailto:mzodet@ahrq.gov> mzodet@ahrq.gov  

        [[alternative HTML version deleted]]



R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sat Jul 24 02:30:50 2004

This archive was generated by hypermail 2.1.8 : Fri 18 Mar 2005 - 02:37:41 EST