I have survey data from a complex survey with stratification, weights and clustering.
I'm using the survey
package in R to run regressions:
svyglm(y ~ x, my_design_object)
Now, I need to add more weights - that is, I have a separate need to weight the respondents, independent of the survey design. Normally, I'd just run:
lm(y ~ x, my_data, weights = my_weights)
Statistically, is it OK to just multiply the old weights by the new weights?
E.g. to do something like
new_design <- svydesign(id = ~id, strata = ~strata,
weights = ~ survey_weights * my_weights,
data = my_design_object$variables)
Or is it (sigh) more complicated than that? And if so, is there an R function that does it?