Tangent exponential model approximation for the GEV distribution


The function gev.tem provides a tangent exponential model (TEM) approximation for higher order likelihood inference for a scalar parameter for the generalized extreme value distribution. Options include location scale and shape parameters as well as value-at-risk (or return levels). The function attempts to find good values for psi that will cover the range of options, but the fail may fit and return an error.


gev.tem(param = c("loc", "scale", "shape", "VaR"), dat, psi = NULL,
  p = NULL, n.psi = 50, plot = TRUE)



parameter over which to profile


sample vector for the GEV distribution


scalar or ordered vector of values for the interest parameter. If NULL (default), a grid of values centered at the MLE is selected


probability associated with the (1-p)th quantile for return levels if param="VaR".


number of values of psi at which the likelihood is computed, if psi is not supplied (NULL). Odd values are more prone to give rise to numerical instabilities near the MLE


logical indiating whether plot.fr should be called upon exit


an object of class fr (see tem) with elements

  • normal: maximum likelihood estimate and standard error of the interest parameter psi

  • par.hat: maximum likelihood estimates

  • par.hat.se: standard errors of maximum likelihood estimates

  • th.rest: estimated maximum profile likelihood at (psi,\hat{λ})

  • r: values of likelihood root corresponding to ψ

  • psi: vector of interest parameter

  • q: vector of likelihood modifications

  • rstar: modified likelihood root vector

  • param: parameter


Leo Belzile, from code by A. Davison from the hoa package

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.