obs_uncertainty: Measurements Uncertainty

Description Usage Arguments Value References Examples

Description

U_obs_95: Measurements uncertainty

RMS_U_obs: Root Mean Square Uncertainty

Usage

1
2
3
4
5
6
7
8
U_obs_95(obs, pollutant = NULL, U_RV95r = NULL, alpha = NULL,
  RV = NULL)

RMS_U_obs(obs, pollutant = NULL, U_RV95r = NULL, alpha = NULL,
  RV = NULL)

U_obs_95_year(obs, pollutant = NULL, U_RV95r = NULL, alpha = NULL,
  RV = NULL, Np = NULL, Nnp = NULL)

Arguments

obs

numeric vector of observed values (yearly averaged for U_obs_95_year)

pollutant

one of "NO2", "O3", "PM10", "PM2.5"

U_RV95r

parameter U^{RV}_{95,r}; optional if pollutant is given, otherwise compulsory

alpha

parameter α; optional if pollutant is given, otherwise compulsory

RV

parameter RV; optional if pollutant is given, otherwise compulsory

Np

parameter N_p; optional if pollutant is given, otherwise compulsory

Nnp

parameter N_np; optional if pollutant is given, otherwise compulsory

Value

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

References

Janssen et al., 2017. "Guidance Document on Modelling Quality Objectives and Benchmarking. Version 2.1"

Examples

 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"))

jobonaf/dartle documentation built on May 29, 2019, 4:52 p.m.