Returns an object of class `"gpdbiv"`

representing the fit
of a bivariate POT (peaks over thresholds) model for joint
excesses over thresholds.

1 2 |

`data1` |
first data vector |

`data2` |
second data vector |

`u1` |
threshold for |

`u2` |
threshold for |

`ne1` |
number of upper extremes to be used for |

`ne2` |
number of upper extremes to be used for |

`global` |
should a global maximisation of the likelihood with respect to marginal and dependence parameters be undertaken. The default alternative is a two-stage local fit where first the marginal parameters are estimated and then the dependence parameter. This is much faster than a global fit. |

`method` |
the optimization method (see |

`...` |
other arguments passed to |

This function implements a model suggested by Richard Smith (see references below). The marginal excess distributions are GPD distributions, as suggested by univariate EVT and implemented in gpd. The dependence specification is known as the logistic or Gumbel dependence structure, but it would be easy to program alternatives.

An object of class `"gpdbiv"`

representing the fit and including
parameter estimates and standard errors.

Multivariate Threshold Methods, Richard L. Smith,
in *Extreme Value Theory and Applications*, ed. J. Galambos,
published by Kluwer, pages 225-248, 1994.

Markov Chain Models for Threshold Exceedances,
R.L. Smith, J.A. Tawn, S.G. Coles,
*Biometrika* **84**, 249-268, 1997.

`gpd`

, `plot.gpdbiv`

,
`interpret.gpdbiv`

1 2 3 4 |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.