留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

一类重尾极值指数估计

上一篇

下一篇

张绿云, 陈守全. 一类重尾极值指数估计[J]. 西南大学学报(自然科学版), 2021, 43(3): 89-94. doi: 10.13718/j.cnki.xdzk.2021.03.013
引用本文: 张绿云, 陈守全. 一类重尾极值指数估计[J]. 西南大学学报(自然科学版), 2021, 43(3): 89-94. doi: 10.13718/j.cnki.xdzk.2021.03.013
ZHANG Lyu-yun, CHEN Shou-quan. A Class of Extreme Value Index Estimators for Heavy-Tailed Distributions[J]. Journal of Southwest University Natural Science Edition, 2021, 43(3): 89-94. doi: 10.13718/j.cnki.xdzk.2021.03.013
Citation: ZHANG Lyu-yun, CHEN Shou-quan. A Class of Extreme Value Index Estimators for Heavy-Tailed Distributions[J]. Journal of Southwest University Natural Science Edition, 2021, 43(3): 89-94. doi: 10.13718/j.cnki.xdzk.2021.03.013

一类重尾极值指数估计

  • 基金项目: 国家自然科学基金项目(11571283)
详细信息
    作者简介:

    张绿云,硕士研究生,主要从事极值统计的研究 .

    通讯作者: 陈守全,博士,副教授
  • 中图分类号: O211.4

A Class of Extreme Value Index Estimators for Heavy-Tailed Distributions

  • 摘要: 提出了如下一类重尾极值指数估计 $ \tilde{\gamma}_{n}^{\left(\beta_{1} \cdot \beta_{2}\right)}(k)=\left\{\left(\beta_{2}-\beta_{1}\right)\left[\frac{{\mathop x\limits^ \wedge }^{\left(\beta_{1}\right)}(k)}{{\mathop x\limits^ \wedge }^{\left(\beta_{2}\right)}(k)}-1\right]^{-1}+1-\beta_{1}\right\}^{-1} $ 其中 $ {\mathop x\limits^ \wedge }^{(\beta)}(k):=\frac{1}{k} \sum\limits_{i=0}^{k-1}\left(\frac{X_{n-i, n}}{X_{n-k, n}}\right)^{1-\beta} $ 证明了其弱相合性和渐近正态性,并在均方误差最小的情况下,给出了 \lt inline-formula \gt $\tilde{\gamma}_{n}^{\left(\beta_{1} \cdot \beta_{2}\right)}(k) $ \lt /inline-formula \gt 中调节参数 \lt i \gt k \lt /i \gt 的最优选择.
  • 加载中
  • [1] FISHER R A, TIPPETT L H C. Limiting Forms of the Frequency Distribution of the Largest or Smallest Member of a Sample [J]. Mathematical Proceedings of the Cambridge Philosophical Society, 1928, 24(2): 180-190. doi: 10.1017/S0305004100015681
    [2] JENKINSON A F. The Frequency Distribution of the Annual Maximum (or Minimum) Values of Meteorological Elements [J]. Quarterly Journal of the Royal Meteorological Society, 1955, 81(348): 158-171. doi: 10.1002/qj.49708134804
    [3] doi: http://biomet.oxfordjournals.org/external-ref?access_num=10.1214/aos/1176343247&link_type=DOI HILL B M. A Simple General Approach to Inference about the Tail of a Distribution [J]. The Annals of Statistics, 1975, 3(5): 1163-1174.
    [4] doi: http://biomet.oxfordjournals.org/external-ref?access_num=10.1214/aos/1176343003&link_type=DOI PICKANDS J. Statistical Inference Using Extreme Order Statistics [J]. The Annals of Statistics, 1975, 3(1): 119-131.
    [5] DEKKERS A L M, EINMAHL J H J, DE HAAN L. A Moment Estimator for the Index of an Extreme-Value Distribution [J]. The Annals of Statistics, 1989, 17(4): 1833-1855.
    [6] BEIRLANT J, VYNCKIER P, TEUGELS J L. Excess Functions and Estimation of the Extreme-Value Index [J]. Bernoulli, 1996, 2(4): 293-318. doi: 10.2307/3318416
    [7] FRAGA ALVES M I, GOMES M I, DE HAAN L, et al. Mixed Moment Estimator and Location Invariant Alternatives [J]. Extremes, 2009, 12(2): 149-185. doi: 10.1007/s10687-008-0073-3
    [8] 王嫣然, 彭作祥. 基于分块的重尾指数估计量的推广[J]. 西南师范大学学报(自然科学版), 2019, 44(1): 25-28. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-XNZK201901006.htm
    [9] BERAN J, SCHELL D, STEHLÍK M. The Harmonic Moment Tail Index Estimator: Asymptotic Distribution and Robustness [J]. Annals of the Institute of Statistical Mathematics, 2014, 66(1): 193-220. doi: 10.1007/s10463-013-0412-2
    [10] DE HAAN L, FERREIRA A. Extreme Value Theory [M]. New York: Springer, 2006.
  • 加载中
计量
  • 文章访问数:  675
  • HTML全文浏览数:  675
  • PDF下载数:  111
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-12-18
  • 刊出日期:  2021-03-20

一类重尾极值指数估计

    通讯作者: 陈守全,博士,副教授
    作者简介: 张绿云,硕士研究生,主要从事极值统计的研究
  • 西南大学 数学与统计学院,重庆 400715
基金项目:  国家自然科学基金项目(11571283)

摘要: 提出了如下一类重尾极值指数估计 $ \tilde{\gamma}_{n}^{\left(\beta_{1} \cdot \beta_{2}\right)}(k)=\left\{\left(\beta_{2}-\beta_{1}\right)\left[\frac{{\mathop x\limits^ \wedge }^{\left(\beta_{1}\right)}(k)}{{\mathop x\limits^ \wedge }^{\left(\beta_{2}\right)}(k)}-1\right]^{-1}+1-\beta_{1}\right\}^{-1} $ 其中 $ {\mathop x\limits^ \wedge }^{(\beta)}(k):=\frac{1}{k} \sum\limits_{i=0}^{k-1}\left(\frac{X_{n-i, n}}{X_{n-k, n}}\right)^{1-\beta} $ 证明了其弱相合性和渐近正态性,并在均方误差最小的情况下,给出了 \lt inline-formula \gt $\tilde{\gamma}_{n}^{\left(\beta_{1} \cdot \beta_{2}\right)}(k) $ \lt /inline-formula \gt 中调节参数 \lt i \gt k \lt /i \gt 的最优选择.

English Abstract

  • 设{Xnn≥1}是一列独立同分布的随机变量序列,其共同的分布函数为F(x),X1,n≤…≤XnnX1,…,Xn的顺序统计量. 若对非退化分布函数G(x),存在规范化常数an>0,$b_{n} \in \mathbb{R} $,对一切G(x)的连续点有

    成立,则称F属于G的吸引场,记为FD(G). 文献[1-2]证明G(x)的形式必为

    极值理论的应用非常广泛,比如在自然灾害、巨额保险赔付以及其他领域如水文学或工程中的罕见事件都时常用到. 形状参数γ被称作极值指数(EVI),它衡量右尾函数$\bar{F}(x):=1-F(x) $的尾部性质. 当分布函数未知时,对极值指数γ的估计是极值理论的重要组成部分. 因此,对于极值指数γ的估计受到学者的广泛关注[3-8].

    在本文中,我们针对极值指数γ>0的分布函数进行研究,这种分布函数也称作重尾分布函数. 如果1-F是指数为$ -\frac{1}{\gamma}$的正则变换(记为$1-F \in R V_{-\frac{1}{\gamma}} $),即

    则称分布函数F是重尾分布函数.

    对于极值指数的研究,前人已经得到较好的结论. 当γ>0时,文献[3]提出了著名的的Hill估计量. 而作为Hill估计的一个推广,文献[9]提出了如下定义的调和矩估计

    其中1≤kn-1且β>0是调谐参数.

    本文受文献[9]中估计量的启发,提出定义如下的一类关于重尾极值指数的新估计

    其中β1β2是调谐参数β的两个取值,且

    当序列k=k(n)→∞满足$ \frac{k}{n} \rightarrow 0$,此时容易得到$ \tilde{\gamma}_{n}^{\left(\beta_{1}, \beta_{2}\right)}(k)$的弱相合性. 此外,存在A(t)使得当t→∞

    成立时,我们考虑$ \tilde{\gamma}_{n}^{\left(\beta_{1}, \beta_{2}\right)}(k)$的渐近性质,其中ρ<0是二阶参数且|A(t)|∈RVρ(参见文献[10],推论2.3.5).

  • 在下文中,我们假设FD(G),γ>0,这等价于$U:=\left(\frac{1}{1-F(x)}\right)^{ \leftarrow } $是指数为γ的正则变换(记为URVγ),其中F(x):=inf{yF(y)≥x}是F的广义逆函数.

    FD(G)当且仅当

    对所有x>0和某些正可测函数a(t)成立. 注意,如果γ=0,则(4)式的右边为lnx.

    为简化定理的表示,本文规定以下记号:

    下面给出本文的主要结果.

    定理1假设条件(4)成立,序列{k(n)}满足k=k(n)→∞,$\frac{k(n)}{n} \rightarrow 0(n \rightarrow \infty) $. 则

    此外,若二阶条件(3)成立,且存在$ \lambda \in \mathbb{R}$使得

    则对$\beta_{1}>1-\frac{1}{2 \gamma}, \beta_{2}>1-\frac{1}{2 \gamma} $,有

    其中

    Vβ1β2(γ)如(5)定义.

    下面考虑序列k的最优选择,寻找使$\widetilde{\gamma}_{n}^{\left(\beta_{1}, \beta_{2}\right)}(k) $的均方误差$\operatorname{MSE}_{\infty}\left(\tilde{\gamma}_{n}^{\left(\beta_{1}, \beta_{2}\right)}(k)\right) $可以达到最小的k.

    定理2  1) 假设条件(3)对ρ<0成立,中间序列{k0(n)}表示使$\operatorname{MSE}_{\infty}\left(\tilde{\gamma}_{n}^{\left(\beta_{1}, \beta_{2}\right)}(k)\right)=A^{2}\left(\frac{n}{k}\right) \mu^{2}+\frac{\sigma^{2}}{k} $达到最小的序列.

    此外,

    其中ssRV2ρ-1的广义逆函数,且

    2) 假设A(t)=C tρ,则1)中的k0

    其中└x┘表示不超过x的最大整数. 此外,

  • Y1Y2,…,Yn是独立同分布的标准Pareto随机变量序列,其共同的分布函数为$ 1-\frac{1}{y}$y≥1,Y1,nY2,n≤…≤Ynn表示Y1Y2,…,Yn的顺序统计量. 易知$\left\{X_{i}\right\}_{i=1}^{n} \stackrel{\mathrm{d}}{=}\left\{U\left(Y_{i}\right)\right\}_{i=1}^{n} $,且$ \left\{\frac{Y_{n-i, n}}{Y_{n-k, n}}\right\}_{i=0}^{k-1} \stackrel{\mathrm{d}}{=}\left\{Y_{k-i, k}\right\}_{i=0}^{k-1}$.

    下面对本文的主要结果进行证明.

    定理1的证明  由二阶条件可得

    因此

    a<1时,$\mathrm{E}\left(Y^{a}\right)=\frac{1}{1-a} $,当$a<\frac{1}{2} $时,$\operatorname{var}\left(Y^{a}\right)=\frac{a^{2}}{(1-a)^{2}(1-2 a)}=: \sigma_{a}^{2} $σa2且当a<1,a+b<1时

    那么对于$ \beta>1-\frac{1}{2 \gamma}$,当n→∞时

    其中N(0,1)表示标准正态分布. 由Pn(β)的定义可得

    那么

    从而

    现定义

    fk(t)表示Qn(β1β2)(k)的特征函数. 由Pn(β)的表达式,有

    其中$\beta_{1}>1-\frac{1}{2 \gamma}, \beta_{2}>1-\frac{1}{2 \gamma} $Vβ1β2(γ)如前定义.

    因此

    若存在$\lambda \in \mathbb{R} $满足$\sqrt{k} A\left(\frac{n}{k}\right) \rightarrow \lambda $,则当n→∞时

    其中

    由于

    故现在利用映射g(θ)来构造一个关于γ的估计量. 应用Delta方法,可得

    其中μσ2如前定义.

    定理2的证明  类似文献[10]中的定理可得结果.

参考文献 (10)

目录

/

返回文章
返回