obs_uncertainty: Measurements Uncertainty

In jobonaf/dartle: Air Quality Model Benchmarking

Description

U_obs_95: Measurements uncertainty

RMS_U_obs: Root Mean Square Uncertainty

Usage

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

# 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.