결측치가 들어있는 mpg 데이터를 활용해서 문제를 해결해보세요.
mpg <- as.data.frame(ggplot2::mpg)
mpg[c(65, 124, 131, 153, 212), "hwy"] <- NA
Q1. drv(구동방식)별로 hwy(고속도로 연비) 평균이 어떻게 다른지 알아보려고 합니다. 분석을 하기 전에 우선 두 변수에 결측치가 있는지 확인해야 합니다. drv 변수와 hwy 변수에 결측치가 몇 개 있는지 알아보세요.
table(is.na(mpg$drv))
##
## FALSE
## 234
table(is.na(mpg$hwy))
##
## FALSE TRUE
## 229 5
• Q2. filter()를 이용해 hwy 변수의 결측치를 제외하고, 어떤 구동방식의 hwy 평균이 높은지 알아보세요. 하나의 dplyr 구문으로 만들어야 합니다.
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
mpg %>%
filter(!is.na(hwy)) %>%
group_by(drv) %>%
summarise(mean_hwy = mean(hwy))
## # A tibble: 3 x 2
## drv mean_hwy
## <chr> <dbl>
## 1 4 19.2
## 2 f 28.2
## 3 r 21
이상치가 들어있는 mpg 데이터를 활용해서 문제를 해결해보세요. 구동방식별로 도시 연비가 다른지 알아보려고 합니다. 분석을 하려면 우선 두 변수에 이상치가 있는지 확인하려고 합니다.
mpg <- as.data.frame(ggplot2::mpg)
mpg[c(10, 14, 58, 93), "drv"] <- "k"
mpg[c(29, 43, 129, 203), "cty"] <- c(3, 4, 39, 42)
Q1. drv에 이상치가 있는지 확인하세요. 이상치를 결측 처리한 다음 이상치가 사라졌는지 확인하세요. 결측 처리 할 때는 %in% 기호를 활용하세요.
table(mpg$drv)
##
## 4 f k r
## 100 106 4 24
mpg$drv <- ifelse(mpg$drv ==c("4", "f", "r"), mpg$drv, NA)
table(mpg$drv)
##
## 4 f r
## 32 35 7
Q2. 상자 그림을 이용해서 cty에 이상치가 있는지 확인하세요. 상자 그림의 통계치를 이용해 정상 범위를 벗어난 값을 결측 처리한 후 다시 상자 그림을 만들어 이상치가 사라졌는지 확인하세요.
boxplot(mpg$cty)$stats
## [,1]
## [1,] 9
## [2,] 14
## [3,] 17
## [4,] 19
## [5,] 26
mpg$cty <- ifelse(mpg$cty < 9 | mpg$cty > 26, NA, mpg$cty)
boxplot(mpg$cty)
Q3. 두 변수의 이상치를 결측처리 했으니 이제 분석할 차례입니다. 이상치를 제외한 다음 drv별로 cty 평균이 어떻게 다른지 알아보세요. 하나의 dplyr구문으로 만들어야 합니다.
mpg %>% filter(!is.na(cty) & !is.na(drv)) %>%
group_by(drv) %>%
summarise(mean_cty = mean(cty))
## # A tibble: 3 x 2
## drv mean_cty
## <chr> <dbl>
## 1 4 14.1
## 2 f 19
## 3 r 13.7