setwd("E:/STAT 50")getwd()## [1] "E:/STAT 50"
library(ggplot2)
library(dplyr)load("brfss2013.RData")brfss2013 %>%
group_by(sleptim1) %>%
filter(is.na(sleptim1)) %>%
summarise(count=n())## # A tibble: 1 × 2
## sleptim1 count
## <int> <int>
## 1 NA 7387
sleepA<-brfss2013 %>%
group_by(sleptim1) %>%
filter(sleptim1 <=5) %>%
summarise(count=n())
sum(sleepA$count)## [1] 52498
sleepF<-brfss2013 %>%
group_by(sleptim1) %>%
filter(sleptim1 >5, sleptim1 <11) %>%
summarise(count=n())
sum(sleepF$count)## [1] 425670
sleepB<-brfss2013 %>%
group_by(sleptim1) %>%
filter(sleptim1 >=11) %>%
summarise(count=n())
sum(sleepB$count)## [1] 6220
sleepM<-brfss2013 %>%
group_by(sleptim1) %>%
filter(sleptim1 <=5, marital == "Married") %>%
summarise(count=n())
sum(sleepM$count)## [1] 21417