I have survey data that needs to be weighted, and to help me with this task, I have access to the full joint distributions of the variables I want to use. As I understand it, I should use the postStratify()
function in the survey
package for this task, but I must admit that I am quite lost when it comes to the required syntax.
I have tried to use marginal distributions and the rake()
function, this is quite straight-forward:
http://www.r-bloggers.com/survey-computing-your-own-post-stratification-weights-in-r/
But how do I incorporate all my joint distributions into the postStratify()
function? The manual only lists a simple example and is a bit vague, to me at least. I have full joint distributions for 5 different variables.
Moreover, if want to use my weighted result with a function that does not take weight as an input, how should I go about and "duplicate" my rows in the final data? When I tried the rake()
function I could only trim weights to an interval with a min value of 0.87. I was thinking I could round all weights to integers.
postStratify()
syntax here. $\endgroup$