##Script for cleaning NHANES data
library(NHANES)
library(tidyverse)
data("NHANES")
##Build Diabetes dataset
##show there is an association between BMI and diabetes
NHANES %>% select(BMI, Diabetes, Gender) %>% ggplot(aes(x=Diabetes, y=BMI)) + geom_boxplot() + facet_wrap(~Gender)
##Bin Data into High (> 30) and Low BMI (<=30)
bmi_diabetes <- NHANES %>% dplyr::select(BMI,Diabetes, Gender, Age) %>% filter(Age > 20) %>%
mutate(BMIstatus=ifelse(BMI > 30, "High", "Low"), BMIstatus=as.factor(BMIstatus))
bmi_diabetes_raw <- bmi_diabetes
bmi_diabetes <- droplevels(bmi_diabetes)
bmi_diabetes %>% ggplot(aes(x=Diabetes, fill=BMIstatus)) +
geom_bar(position="fill", color="black") + facet_wrap(~Gender)
devtools::use_data(bmi_diabetes, overwrite=TRUE)
devtools::use_data(bmi_diabetes_raw, overwrite=TRUE)
##Whickham dataset
library(mosaic)
data("Whickham")
devtools::use_data(Whickham, overwrite=TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.