Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- lm(DR1TKCAL ~ population + sex + age, data = combined.data)
- diet_design <- svydesign(id = ~ SDMVPSU, strata = ~ SDMVSTRA, nest = TRUE, weight = ~ WTDRD1, data = diet_data)
- svyby(~DR1TKCAL, ~sex, diet_design, svymean, na.rm=TRUE, keep.names=FALSE)
Add Comment
Please, Sign In to add comment