fit.rho: Estimator of the second order tail index parameter

In lbelzile/mev: Modelling of Extreme Values

Estimator of the second order tail index parameter

Description

Estimator of the second order tail index parameter

Usage

fit.rho(xdat, k, method = c("fagh", "dk", "ghp"), ...)

Arguments

xdat	vector of positive observations
k	number of highest order statistics to use for estimation
method	string for the estimator
...	additional arguments passed to individual routinescurrently ignored.

Examples

# Example with rho = -0.2
n <- 1000
xdat <- mev::rgp(n = n, shape = 0.2)
kmin <- floor(n^0.995)
kmax <- ceiling(n^0.999)
rho_est <- fit.rho(
   xdat = xdat,
   k = n - kmin:kmax)
rho_med <- mean(rho_est$rho)

lbelzile/mev documentation built on June 14, 2025, 6:40 p.m.