[R] sampling design runs with no errors but returns empty data set

From: Simon Kiss <sjkiss_at_gmail.com>
Date: Wed, 30 Mar 2011 11:01:15 -0400


Dear colleagues,
I'm working with the 2008 Canada Election Studies (
http://www.queensu.ca/cora/_files/_CES/CES2008.sav.zip), trying to construct a weighted national sample using the survey package. Three weights are included in the national survey (a household weight, a provincial weight and a national weight which is a product of the first two). In the following code I removed variables with missing national weights and tried to construct the sample from advice I've gleaned from the documentation for the survey package and other help requests. There are no errors, but the data frame (weight_test) contains no What am I missing?
Yours, Simon Kiss
P.S. The code is only reproducible if the data set is downloadable. I'm nt sure

ces<-read.spss(file.choose(), to.data.frame=TRUE, use.value.labels=FALSE) missing_data<-subset(ces1, !is.na(ces08_NATWGT)) weight_test<-svydesign(id=~0, weights=~ces08_NATWGT, data=missing_data)

Note: this is some reproducible code that creates a data set that is a very stripped down version of what I'm working with, but with this, the surveydesign function appears to work properly.

mydat<-data.frame(ces08_HHWGT=runif(3000, 0.5, 5), ces08_PROVWGT=runif(3000, 0.6, 1.2), party=sample(c("NDP", "BQ", "Lib", "Con"), 3000, replace=TRUE), age=sample(seq(18, 72,1), 3000, replace=TRUE), income=sample(seq(21,121,1), 3000, replace=TRUE)) mydat$ces08_NATWGT<-mydat$ces08_HHWGT*mydat$ces08_PROVWGT weight_test<-svydesign(id=~1, weights=~ces08_NATWGT, data=mydat)



Simon J. Kiss, PhD
Assistant Professor, Wilfrid Laurier University 73 George Street
Brantford, Ontario, Canada
N3T 2C9
Cell: +1 519 761 7606

R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Wed 30 Mar 2011 - 15:11:40 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 30 Mar 2011 - 20:20:25 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive