Learn R Programming

gamair (version 1.0-2)

swer: 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)

Arguments

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.

Examples

Run this code
# NOT RUN {
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)
# }

Run the code above in your browser using DataLab