Weighted Data with Survey Package

[ella1]

Hey !

My dataset includes 5 variables:

gender - 1 (Male), 2 (Female)

car (has a car) - 1 (Yes), 2 (No), 3 (Don't know), 4 (Refused)

motorcycle (has a motorcycle) - 1 (Yes), 2 (No), 3 (Don't know), 4 (Refused)

helmet (if has a motorcycle: wears helmet when driving a motorcycle) - 1 (Yes), 2 (No), 3 (Don't know), 4 (Refused)

weights - sampling weights

I want to do some descriptive statistics on:

1: proportion of females having a car (among women)

2: proportion of people wearing a helmet when driving a motorcycle among those having a motorcycle

3: proportion of women wearing a helmet when driving a motorcycle among women having a motorcycle

4: proportion of women wearing a helmet when driving a motorcycle among women having both a car and a motorcycle

I am using the survey package. This is what I did:

weighteddata = svydesign(ids = ~1, data = data, weights = data$weights)

for 1: svyby(~car, ~gender, svymean)

With 2 - 4 I struggle.

Please send help!

Thanks, Ella

[williaml]

Can you provide a reproducible example for your dataset?

FAQ: What's a reproducible example (`reprex`) and how do I create one?

Why reprex? Getting unstuck is hard. Your first step here is usually to create a reprex, or reproducible example. The goal of a reprex is to package your code, and information about your problem so that others can run it and feel your pain. Then, hopefully, folks can more easily provide a solution. What's in a Reproducible Example? Parts of a reproducible example: background information - Describe what you are trying to do. What have you already done? complete set up - include any library() calls and data to reproduce your issue. data for a reprex: Here's a discussion on setting up data for a reprex make it run - include the minimal code required to reproduce your error on the data…