資料為4/14~4/28日各地測站收集每小時空污指標濃度後換算為每日平均濃度
PM10指粒徑在10微米以下之粒子，又稱浮游塵。 主要來源包括道路揚塵、車輛排放廢氣、露天燃燒、營建施工及農地耕作等，或由原生性空氣污染物轉化成之二次污染物，由於粒徑小於10微米以下，能深入人體肺部深處，如該粒子附著其他污染物，則將加深對呼吸系統之危害。
rm(list = ls(all=TRUE))
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
library(knitr)
rawdata = read.csv(
'airpollutiontest.csv',
encoding = 'UTF-8')
#各地測站懸浮微粒日平均濃度
df <- rawdata %>% group_by(SiteName)
result1 = summarise(df, mean(PM10SubIndex))
write.table(result1, "result1")
data = read.table("result1")
kable(data)
| 安南 |
48.00000 |
| 板橋 |
NA |
| 埔里 |
36.71429 |
| 菜寮 |
42.92308 |
| 潮州 |
51.61538 |
| 大里 |
45.61538 |
| 大寮 |
NA |
| 大同 |
52.84615 |
| 大園 |
46.92308 |
| 淡水 |
44.69231 |
| 冬山 |
39.00000 |
| 斗六 |
47.46154 |
| 二林 |
47.53846 |
| 豐原 |
34.61538 |
| 鳳山 |
51.69231 |
| 復興 |
40.07143 |
| 古亭 |
38.46154 |
| 關山 |
25.50000 |
| 觀音 |
50.30769 |
| 恆春 |
25.46154 |
| 湖口 |
41.00000 |
| 花蓮 |
32.46154 |
| 基隆 |
29.84615 |
| 嘉義 |
NA |
| 金門 |
NA |
| 林口 |
55.15385 |
| 林園 |
26.92308 |
| 龍潭 |
47.76923 |
| 崙背 |
46.46154 |
| 馬公 |
33.92857 |
| 馬祖 |
52.35714 |
| 麥寮 |
55.92857 |
| 美濃 |
38.38462 |
| 苗栗 |
44.07692 |
| 南投 |
46.61538 |
| 楠梓 |
55.53846 |
| 平鎮 |
47.69231 |
| 屏東 |
43.84615 |
| 朴子 |
53.61538 |
| 前金 |
45.69231 |
| 前鎮 |
42.00000 |
| 橋頭 |
45.53846 |
| 仁武 |
49.30769 |
| 三義 |
35.30769 |
| 三重 |
56.07143 |
| 沙鹿 |
41.00000 |
| 善化 |
44.30769 |
| 士林 |
39.38462 |
| 松山 |
39.69231 |
| 臺東 |
28.84615 |
| 臺南 |
41.23077 |
| 臺西 |
NA |
| 桃園 |
44.84615 |
| 頭份 |
NA |
| 土城 |
39.15385 |
| 萬華 |
40.30769 |
| 萬里 |
45.07692 |
| 汐止 |
38.38462 |
| 西屯 |
44.61538 |
| 線西 |
49.46154 |
| 小港 |
45.30769 |
| 新店 |
36.92308 |
| 新港 |
49.38462 |
| 新營 |
52.30769 |
| 新竹 |
NA |
| 新莊 |
45.23077 |
| 陽明 |
18.00000 |
| 宜蘭 |
35.92857 |
| 永和 |
42.35714 |
| 彰化 |
45.53846 |
| 中壢 |
NA |
| 中山 |
41.46154 |
| 忠明 |
39.23077 |
| 竹東 |
40.30769 |
| 竹山 |
48.28571 |
| 左營 |
46.84615 |
#4/14~28所有測站整體平均
mean1=mean(rawdata$PM10SubIndex,na.rm = T)
mean1
## [1] 43.40142
#每日懸浮微粒平均濃度
df <- rawdata %>% group_by(MonitorDate)
result11 = summarise(df, mean(PM10SubIndex))
write.table(result11, "result11")
data = read.table("result11")
kable(data)
| 2017-04-14 |
46.53846 |
| 2017-04-15 |
NA |
| 2017-04-16 |
NA |
| 2017-04-17 |
NA |
| 2017-04-18 |
NA |
| 2017-04-19 |
NA |
| 2017-04-20 |
NA |
| 2017-04-21 |
32.01316 |
| 2017-04-22 |
37.51316 |
| 2017-04-23 |
45.68421 |
| 2017-04-24 |
NA |
| 2017-04-25 |
NA |
| 2017-04-27 |
NA |
| 2017-04-28 |
NA |
#不符合台灣標準之測站地點：24小時平均值不得超過125 μg／m3
dfTWSd <- rawdata %>% group_by(SiteName)
result2 = dfTWSd %>% filter(PM10SubIndex>125)
write.table(result2, "result2")
data = read.table("result2")
kable(data)
| SiteId |
SiteName |
MonitorDate |
PSI |
SO2SubIndex |
COSubIndex |
O3SubIndex |
PM10SubIndex |
NO2SubIndex |
#不符合歐盟標準之測站地點：24小時平均值不得超過50 μg／m3
dfEUSd <- rawdata %>% group_by(SiteName)
result3 = dfEUSd %>% filter(PM10SubIndex>50)
EUPM10=result3 %>% group_by(SiteName)
result4=summarise(EUPM10)
write.table(result4, "result4")
data = read.table("result4")
kable(data)
SiteName
安南
板橋
埔里
菜寮
潮州
大里
大寮
大同
大園
淡水
斗六
二林
豐原
鳳山
復興
古亭
觀音
湖口
嘉義
金門
林口
龍潭
崙背
馬公
馬祖
麥寮
美濃
苗栗
南投
楠梓
平鎮
屏東
朴子
前金
前鎮
橋頭
仁武
三義
三重
沙鹿
善化
士林
松山
臺南
臺西
桃園
頭份
土城
萬華
萬里
西屯
線西
小港
新港
新營
新竹
新莊
永和
彰化
中壢
中山
忠明
竹東
竹山
左營