Hello,
I just want to briefly mention that I’m fairly new to R, and this is the first analysis I try to run in the context of a complex survey design. Below a description of my problem along with a few questions.
I want to output a one-way frequency analysis. The data set contains survey weights (WTS_S), along with 500 bootstrap weights (BSW1-BSW500). The variable of interest is PROVINCE. The method to estimate variance is BRR. I’ using the following code to capture the survey design.
CCHSDesign = svrepdesign(data = myData,
weights = ~WTS_S,
repweights = "BSW[1-500]",
combined.weights = F,
type = "BRR")
I use combined.weights = F since repweights does not contain the sampling weights. Neverthelss, I get the following warning message:
Warning message:
In svrepdesign.default(data = myData, weights = ~WTS_S, repweights = "BSW[1-500]", :
Data look like combined weights: mean replication weight is 1303.63701983038 and mean sampling weight is 1302.86018957346
The previous message does no appear when I used combined.weights = F, which seems not to be appropriate in my situation.
Question 1: Is the previous warning message of any concern to me?
Question 2: Why do I get “Balanced Repeated Replicates with 456 replicates.” instead of 500 when I run CCHSDesign?
In addition to the frequencies by category of PROVINCE, I also want the standard errors (SEs). Running the following code provides the frequencies, but not the SEs.
prop.table(svytable(~PROVINCE,design=CCHSDesign))*100
Question 3: How can I get the SEs in addition to the frequencies by categories of PROVINCE?
Thanks a lot for your support,
A.G.