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@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.
