The Rmetrics "fExtremes" package is a collection of functions to analyze and model extreme events in Finance and Insurance.
The fExtremes
package provides functions for analyzing
and modeling extreme events in financial time Series. The
topics include: (i) data pre-processing, (ii) explorative
data analysis, (iii) peak over threshold modeling, (iv) block
maxima modeling, (v) estimation of VaR and CVaR, and (vi) the
computation of the extreme index.
Data Sets:
Data sets used in the examples of the timeSeries packages.
Data Preprocessing:
These are tools for data preprocessing, including functions to separate data beyond a threshold value, to compute blockwise data like block maxima, and to decluster point process data.
blockMaxima extracts block maxima from a vector or a time series
findThreshold finds upper threshold for a given number of extremes
pointProcess extracts peaks over Threshold from a vector or a time series
deCluster de-clusters clustered point process data
This section contains a collection of functions for explorative data analysis of extreme values in financial time series. The tools include plot functions for empirical distributions, quantile plots, graphs exploring the properties of exceedances over a threshold, plots for mean/sum ratio and for the development of records. The functions are:
emdPlot plots of empirical distribution function
qqparetoPlot exponential/Pareto quantile plot
mePlot plot of mean excesses over a threshold
mrlPlot another variant, mean residual life plot
mxfPlot another variant, with confidence intervals
msratioPlot plot of the ratio of maximum and sum
recordsPlot Record development compared with iid data
ssrecordsPlot another variant, investigates subsamples
sllnPlot verifies Kolmogorov's strong law of large numbers
lilPlot verifies Hartman-Wintner's law of the iterated logarithm
xacfPlot plots ACF of exceedances over a threshold
Parameter Fitting of Mean Excesses:
normMeanExcessFit fits mean excesses with a normal density
ghMeanExcessFit fits mean excesses with a GH density
hypMeanExcessFit fits mean excesses with a HYP density
nigMeanExcessFit fits mean excesses with a NIG density
ghtMeanExcessFit fits mean excesses with a GHT density
GPD Distribution:
A collection of functions to compute the generalized Pareto distribution. The functions compute density, distribution function, quantile function and generate random deviates for the GPD. In addition functions to compute the true moments and to display the distribution and random variates changing parameters interactively are available.
dgpd returns the density of the GPD distribution
pgpd returns the probability function of the GPD
qgpd returns quantile function of the GPD distribution
rgpd generates random variates from the GPD distribution
gpdSlider displays density or rvs from a GPD
GPD Moments:
gpdMoments computes true mean and variance of GDP
GPD Parameter Estimation:
This section contains functions to fit and to simulate processes that are generated from the generalized Pareto distribution. Two approaches for parameter estimation are provided: Maximum likelihood estimation and the probability weighted moment method.
gpdSim generates data from the GPD distribution
gpdFit fits data to the GPD istribution
GPD print, plot and summary methods:
print print method for a fitted GPD object
plot plot method for a fitted GPD object
summary summary method for a fitted GPD object
GDP Tail Risk:
The following functions compute tail risk under the GPD approach.
gpdQPlot estimation of high quantiles
gpdQuantPlot variation of high quantiles with threshold
gpdRiskMeasures prescribed quantiles and expected shortfalls
gpdSfallPlot expected shortfall with confidence intervals
gpdShapePlot variation of GPD shape with threshold
gpdTailPlot plot of the GPD tail
GEV Distribution:
This section contains functions to fit and to simulate processes that are generated from the generalized extreme value distribution including the Frechet, Gumbel, and Weibull distributions.
dgev returns density of the GEV distribution
pgev returns probability function of the GEV
qgev returns quantile function of the GEV distribution
rgev generates random variates from the GEV distribution
gevSlider displays density or rvs from a GEV
GEV Moments:
gevMoments computes true mean and variance
GEV Parameter Estimation:
A collection to simulate and to estimate the parameters of processes generated from GEV distribution.
gevSim generates data from the GEV distribution
gumbelSim generates data from the Gumbel distribution
gevFit fits data to the GEV distribution
gumbelFit fits data to the Gumbel distribution
print print method for a fitted GEV object
plot plot method for a fitted GEV object
summary summary method for a fitted GEV object
GEV MDA Estimation:
Here we provide Maximum Domain of Attraction estimators and visualize the results by a Hill plot and a common shape parameter plot from the Pickands, Einmal-Decker-deHaan, and Hill estimators.
hillPlot shape parameter and Hill estimate of the tail index
shaparmPlot variation of shape parameter with tail depth
GEV Risk Estimation:
gevrlevelPlot k-block return level with confidence intervals
Two functions to compute Value-at-Risk and conditional Value-at-Risk.
VaR computes Value-at-Risk
CVaR computes conditional Value-at-Risk
A collection of functions to simulate time series with a known extremal index, and to estimate the extremal index by four different kind of methods, the blocks method, the reciprocal mean cluster size method, the runs method, and the method of Ferro and Segers.
thetaSim simulates a time Series with known theta
blockTheta computes theta from Block Method
clusterTheta computes theta from Reciprocal Cluster Method
runTheta computes theta from Run Method
ferrosegersTheta computes theta according to Ferro and Segers
exindexPlot calculates and plots Theta(1,2,3)
exindexesPlot calculates Theta(1,2) and plots Theta(1)
The fExtremes
Rmetrics package is written for educational
support in teaching "Computational Finance and Financial Engineering"
and licensed under the GPL.
Package: fExtremes
Type: Package
License: GPL Version 2 or later
Copyright: (c) 1999-2014 Rmetrics Association
URL: https://www.rmetrics.org