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.