| swer {gamair} | R Documentation |
Swiss 12 hour extreme rainfall
Description
Records the most extreme 12 hourly total rainfall each year for 65 Swiss weather stations. The data period is 1981-2015, although not all stations start in 1981.
Usage
data(swer)
Format
The swer data frame has the following columns
- year
The year of observation.
- exra
The highest rainfall observed in any 12 hour period in that year, in mm.
- nao
Annual North Atlantic Oscillation index, based on the difference of normalized sea level pressure (SLP) between Lisbon, Portugal and Stykkisholmur/Reykjavik, Iceland. Positive values are generally associated with wetter and milder weather over Western Europe.
- location
The measuring station location name.
- code
Three letter code identifying the station.
- elevation
metres above sea level.
- climate.region
One of 12 distinct climate regions.
- N
Degrees north.
- E
Degrees east.
Details
The actual extreme rainfall measurements are digitized from plots in the MeteoSwiss reports for each station. The error associated with digitization can be estimated from the error in the digitized year values, since the true values are then known exactly. This translates into a mean square error in rainfall of about 0.1% of the station maximum, and a maximum error of about 0.3% of station maximum.
Source
Mostly from the MeteoSwiss website:
NAO data from:
Hurrell, James & National Center for Atmospheric Research Staff (Eds). Last modified 16 Aug 2016. "The Climate Data Guide: Hurrell North Atlantic Oscillation (NAO) Index (station-based)."
Examples
require(gamair);require(mgcv)
data(swer)
## GEV model, over-simplified for speed...
system.time(b <- gam(list(exra~s(elevation)+ climate.region,
~s(elevation),~1),family=gevlss,data=swer))
plot(b,pages=1,scale=0,scheme=1)