Description Usage Arguments Value References Examples
U_obs_95
: Measurements uncertainty
RMS_U_obs
: Root Mean Square Uncertainty
1 2 3 4 5 6 7 8 |
obs |
numeric vector of observed values (yearly averaged for |
pollutant |
one of |
U_RV95r |
parameter U^{RV}_{95,r}; optional if |
alpha |
parameter α; optional if |
RV |
parameter RV; optional if |
Np |
parameter N_p; optional if |
Nnp |
parameter N_np; optional if |
U_obs_95
: measurements uncertainty U_95,
as in eq.10, p.20 in Janssen et al., 2017
RMS_U_obs
: root mean square uncertainty RMS_U,
as in eq.11, p.20 in Janssen et al., 2017
U_obs_95_year
: measurements uncertainty for yearly averaged values U_95,
as in eq.12, p.20 in Janssen et al., 2017
Janssen et al., 2017. "Guidance Document on Modelling Quality Objectives and Benchmarking. Version 2.1"
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | # concentrations
cc <- 1:400
# uncertainties for time series
Uo3 <- U_obs_95(cc, pollutant = "O3")
Upm10 <- U_obs_95(cc, pollutant = "PM10")
Upm25 <- U_obs_95(cc, pollutant = "PM2.5")
Uno2 <- U_obs_95(cc, pollutant = "NO2")
plot(cc, Uo3, type="l", log = "xy", xlab="measured concentration",
ylab="uncertainty", ylim=c(1,100))
lines(cc,Upm10,lty=2)
lines(cc,Upm25,lty=3)
lines(cc,Uno2,lty=4)
legend("topleft",lty=1:4,legend=c("O3", "PM10", "PM2.5", "NO2"))
# uncertainties for yearly averages
Uo3 <- U_obs_95_year(cc, pollutant = "O3")
Upm10 <- U_obs_95_year(cc, pollutant = "PM10")
Upm25 <- U_obs_95_year(cc, pollutant = "PM2.5")
Uno2 <- U_obs_95_year(cc, pollutant = "NO2")
plot(cc, Uo3, type="l", log = "xy", xlab="measured concentration",
ylab="uncertainty", ylim=c(1,100))
lines(cc,Upm10,lty=2)
lines(cc,Upm25,lty=3)
lines(cc,Uno2,lty=4)
legend("topleft",lty=1:4,legend=c("O3", "PM10", "PM2.5", "NO2"))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.