從 R 中同一圖表上的數據集創建多個箱線圖

Question

下表是我想用來制作三個箱線圖的數據集的一小部分：

##         No year month day hour PM2.5  PM10   SO2   NO2   CO    O3  TEMP   PRES
## 1     7345 2014     1   1    0  20.0  90.0  18.0  62.0   NA    NA  -1.5 1007.3
## 2     7346 2014     1   1    1  43.0 348.0  25.0  91.0 1100   1.0  -2.6 1006.9
## 3     7347 2014     1   1    2  79.0 423.0  41.0 103.0 1800   1.0  -3.0 1006.9
## 4     7348 2014     1   1    3  82.0 337.0  43.0 101.0 2100   1.0  -3.3 1006.4
## 5     7349 2014     1   1    4 124.0 594.0  59.0 130.0 2400   1.0  -2.7 1006.1
## 6     7350 2014     1   1    5  89.0 307.0  47.0 102.0 2500   1.0  -3.1 1006.6
## 7     7351 2014     1   1    6  59.0 161.0  45.0  91.0 1900   1.0  -2.6 1007.2
## 8     7352 2014     1   1    7  31.0  93.0  24.0  69.0  900   4.0  -2.9 1007.9
## 9     7353 2014     1   1    8  21.0  90.0  18.0  52.0  700  17.0   0.0 1008.8
## 10    7354 2014     1   1    9  38.0 142.0  25.0  87.0 1200   1.0   9.0 1009.6

我想制作變量 PM2.5 的三個箱線圖，它們都具有條件小時 == 12，但 2014、2015 和 2016 年各有一個。它們還需要在同一個圖表上。 這是我嘗試過的：

plots <- data.frame(box14 = subset(dat, hour == 12 & year == 2014)$PM2.5,
                    box15 = subset(dat, hour == 12 & year == 2015)$PM2.5,
                    box16 = subset(dat, hour == 12 & year == 2016)$PM2.5)
boxplot(plots)

我只是得到'data.frame 中的錯誤...... arguments 暗示不同的行數'。 我應該改變什么來獲得我想要的箱線圖？

Answer 1

這是使用dplyr和ggplot2的快速解決方案：

library(dplyr)
library(ggplot2)

dat %>% 
  as_tibble() %>% 
  filter(hour == 12) %>% 
  filter(year %in% c(2014, 2015, 2016)) %>% 
  ggplot(aes(x=as.factor(year), y=PM2.5)) + 
  geom_boxplot()

Answer 2

在基礎 r 中：

filtered_data <- original_data[original_data$hour == 12 & original_data$year %in% c(2014,2015,2016),]

接着

boxplot(PM2.5~year,data=filtered_data)