从R中同一图形上的数据集创建多个箱线图

下表是我想用来绘制三个箱线图的数据集的一小部分：

##         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中的错误…”。。。表示不同行数的参数'。我应该更改什么以获得所需的箱线图？

这里有一个使用dplyr和ggplot2的快速解决方案：

---

下面是一个使用dplyr和ggplot2的快速解决方案：

在r基中：

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

在r基中：

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

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