The package SpatialModelsZAMG was developed for spatial modeling of extreme snow depth (sd) and snow water equivalent (swe) in Austria, but can in principle be used for any kind of extreme events.
There is a function for model selection, and several functions for optimization of the model coefficients, like via smooth modeling, bivariate Huesler-Reiss, bivariate Extremal-Gaussian or bivariate Extremal-t model.
Furthermore, return levels can be calculated and maps can be plotted.
For the bivariate models, there is also the possibility to calculate and plot conditional return levels.
The package provides the following functions:
get_data_from_Robj
: extract the snow data from a saved robj
file
model_selection
: find best linear models according to AIC
optimizer_smooth_model
: optimize the coefficients of the best fitted linear models via smooth modeling
optimizer_biv_hr_model
: optimize the coefficients of the best fitted linear models via bivariate Huesler-Reiss model
optimizer_biv_ext_gauss_model
: optimize the coefficients of the best fitted linear models via bivariate Extremal-Gaussian model
optimizer_biv_ext_t_model
: optimize the coefficients of the best fitted linear models via bivariate Extremal-t model
GEVparameters_from_models
: calculate GEV parameters according to given linear models
GEVparameters_map
: create GEV parameter maps for Austria
returnlevels
: calculate return levels
returnlevel_plot
: create return level plots
returnlevel_map
: create return level maps for Austria
cond_returnlevels
: calculate conditional return levels
cond_returnlevel_plot
: create conditional return level plots
cond_returnlevel_map
: create conditional return level maps for Austria
model_difference_map
: create model difference maps for Austria
The package also includes some data files, which are used in the example codes:
‘lon.at’: longitude data from Austria
‘lat.at’: latitude data from Austria
‘alt.at’: altitude data from Austria
‘border.at’: border data from Austria
‘sample_grid_data’: sample grid data for sd and swe
‘sample_snow_data’: sample snow data from 36 examplary stations distributed over Austria
‘get_data’: exemplary function output from get_data_from_Robj
‘sd_m_select’: exemplary function output from model_selection
for sd
‘swe_m_select’: exemplary function output from model_selection
for swe
‘optim_hr’: exemplary function output from optimizer_biv_hr_model
‘sd_GEVparam’: exemplary function output from GEVparameters_from_models
for sd
‘swe_GEVparam’: exemplary function output from GEVparameters_from_models
for swe
Simon Gstoehl, simon.gstoehl@vol.at
Contact:
University of Innsbruck: Tobias Hell, tobias.hell@uibk.ac.at
ZAMG Innsbruck: Harald Schellander, harald.schellander@zamg.ac.at
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.