Description Usage Arguments Value Author(s) References See Also Examples
Given an ordered sample of either exceedances or upper order statistics which is to be modeled using a GPD with distribution function F, this function provides Segers' estimator of the shape parameter γ, see Segers (2005). Precisely, for k = \{1, …, n-1\}, the estimator can be written as
\hat γ^k_{\rm{Segers}}(H) = ∑_{j=1}^k \Bigl(λ(j/k) - λ((j-1)/k)\Bigr) \log \Bigl(H^{-1}((n-\lfloor cj \rfloor)/n)-H^{-1}((n-j)/n) \Bigr)
for H either the empirical or the distribution function based on the log–concave density estimator and λ the mixing measure given in Segers (2005), Theorem 4.1, (i). Note that for any k, \hat γ^k_{\rm{Segers}} : R^n \to (-∞, ∞). If \hat γ_{\rm{Segers}} \not \in [-1,0), then it is likely that the log-concavity assumption is violated.
1 | generalizedPick(est, c, gam0, ks = NA)
|
est |
Log-concave density estimate based on the sample as output by |
c |
Number in (0,1), determining the spacings that are used. |
gam0 |
Number in R \setminus 0.5, specifying the mixing measure. |
ks |
Indices k at which Falk's estimate should be computed. If set to |
n x 3 matrix with columns: indices k, Segers' estimator using the smoothing method, and the ordinary Segers' estimator based on the order statistics.
Kaspar Rufibach (maintainer), kaspar.rufibach@gmail.com,
http://www.kasparrufibach.ch
Samuel Mueller, samuel.mueller@sydney.edu.au,
www.maths.usyd.edu.au/ut/people?who=S_Mueller
Kaspar Rufibach acknowledges support by the Swiss National Science Foundation SNF, http://www.snf.ch
Mueller, S. and Rufibach K. (2009). Smooth tail index estimation. J. Stat. Comput. Simul., 79, 1155–1167.
Segers, J. (2005). Generalized Pickands estimators for the extreme value index. J. Statist. Plann. Inference, 128, 381–396.
Other approaches to estimate γ based on the fact that the density is log–concave, thus
γ \in [-1,0], are available as the functions pickands
, falk
, falkMVUE
.
1 2 3 4 5 6 7 8 9 10 11 |
Loading required package: logcondens
Warning messages:
1: In rgl.init(initValue, onlyNULL) : RGL: unable to open X11 display
2: 'rgl_init' failed, running with rgl.useNULL = TRUE
3: .onUnload failed in unloadNamespace() for 'rgl', details:
call: fun(...)
error: object 'rgl_quit' not found
k logcon order
[1,] 1 0.00000000 0.00000000
[2,] 2 0.08432864 0.07351167
[3,] 3 0.09070362 -0.41318757
[4,] 4 0.09789028 -0.22954196
[5,] 5 0.01234323 -0.23533598
[6,] 6 -0.03262626 -0.06223357
[7,] 7 0.02163967 -0.08026050
[8,] 8 -0.01932424 -0.10072821
[9,] 9 -0.03159206 -0.11917529
[10,] 10 -0.04979492 -0.17685572
[11,] 11 -0.06656938 0.07462661
[12,] 12 -0.03467781 -0.17065951
[13,] 13 -0.10350782 -0.23046877
[14,] 14 -0.07534733 -0.01223848
[15,] 15 -0.08980431 -0.23660671
[16,] 16 -0.08342789 -0.21717061
[17,] 17 -0.18599160 -0.34169971
[18,] 18 -0.17658036 -0.19632088
[19,] 19 -0.21815990 -0.24165160
[20,] 20 NA NA
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.