The UNprecedented Simulated Extreme ENsemble (UNSEEN, Thompson et al., 2017) approach is an increasingly popular method that exploits seasonal prediction systems to assess and anticipate climate extremes beyond the observed record. The goal of this package is to provide easy evaluation of the UNSEEN method. This package is developed as part of the ECMWF Summer of Weather Code (ESoWC 2020) using in particular for the ECMWF SEAS5 hindcast data.
You can install the development version from GitHub with:
# install.packages("devtools")
devtools::install_github("timokelder/UNSEEN")
This is a basic example which shows you how to plot the UNSEEN timeseries and apply the independence, stability and fidelity tests using an example dataset:
data("EOBS_UK","SEAS5_UK")
library(UNSEEN)
## basic example code
unseen_timeseries(ensemble = SEAS5_UK, obs = EOBS_UK, ylab = 'UK February precipitation (mm/d)')
independence_test(ensemble = SEAS5_UK)
#> Warning: Removed 1625 rows containing non-finite values (stat_ydensity).
#> Warning: Removed 1625 rows containing non-finite values (stat_boxplot).
stability_test(ensemble = SEAS5_UK, lab = 'UK February precipitation (mm/d)')
#> Warning: Removed 4 row(s) containing missing values (geom_path).
fidelity_test(obs = EOBS_UK$rr, ensemble = SEAS5_UK$tprate)