甘胜进;游文杰;涂开仁
【摘 要】Copula function can relate marginal distribution with joint distribution .It is an effective tool to describe the dependency between variables .Gaussian distribution plays a important role in practice .This paper studied the extreme tail dependence coefficient and extreme lower tail Copula of Gussian distribution .Gaussian copula's extreme value Copula is derived .%Copula函数把边缘分布函数与联合分布函数联系起来,是研究变量间相依性的一种有效工具。而高斯分布在实际应用中占有重要的地位,本文主要研究高斯分布的极尾相依系数、极高斯Copula函数,推导出极值高斯Copula函数。
【期刊名称】《佳木斯大学学报(自然科学版)》
【年(卷),期】2015(000)001
【总页数】3页(P149-151)
【关键词】高斯Copula;下尾相依系数;极下尾Copula;极值Copula
【作 者】甘胜进;游文杰;涂开仁
【作者单位】福建师范大学福清分校电子与信息工程学院 福建 福清 350300;福建师范大学福清分校电子与信息工程学院 福建 福清 350300;福建师范大学福清分校经济与管理学院,福建 福清350300
【正文语种】中 文
【中图分类】O213
0 引 言
p 维连续型随机变量X=(X1,X2,…,Xp)T 的联合分布函数为F(x1,x2,…,xp),边缘分布函数分别为F1(x1),F2(x2),…,Fp(xp),则有:
其中诸Fi(Xi)~U[0,1],Fi(xi)=ui,(u1,u2,…,up)∈[0,1]p.由此可见Copula 函数把联合分布函数与其边缘分布函数连接起来,是[0,1]上均匀分布随机变量的联合分布函数,能够有效研究变量间相依性,其在经济金融、风险管理等方面有着这广泛的应用.高斯Copula 函数是多元正态分布的Copula 形式,由于随机变量的单调递增变换不会改变其Copula 函数形式,故不妨设X ~Np(0,R),其中R(>0)为相关系数矩阵,φ(x),φ(x)分别表示标准正态分布的分布函数和密度函数,则高斯Copula 函数为:
其Copula 密度为:
X 的联合密度函数:
1 利用高斯Copula 函数研究二元正态变量之间的相关性
为计算方便,不妨对二元高斯Copula 进行探讨,类似的方法可以推广到多元,即
1.1 尾相依系数
称p(X2 ≤x|X1 ≤x)为度量在随机变量X1 实现条件下X2 实现的概率大小,当x →-∞称此时极限为极下尾相依系数(TDC),记λL 为,为了便于利用Copula 函数计算,假设X1,X2 具有相同的边缘分布为F(x),则:
同理可以定义上极尾相依系数其中为(F1(X1),F2(X2))的生存Copula 函数,也是(1-F1(X1),1-F2(X2))的分布函数,与C(u,v)之间关系是v-1+C(1-u,1-v).当(X1,X2)是对称分布时,λL=λu,故对于高斯Copula 只讨论极下尾相依系数:
性质1:λL=0.证明:
说明两正态分布在尾部是渐进独立的.
性质2: 当-1 ≤ρ1 ≤ρ2 ≤1,Cρ1(u,v)≤Cρ2(u,v),并且C-1(u,v)=max(u+v-1,0),C1(u,v)=min(u,v),进而当ρ ≥0,uv ≤Cρ(u,v)≤min(u,v).
此性质的证明需要用到以下引理:
引 理[2]:(X1,X2,…,Xp)T, (X1′,X2′,…,Xp′)T 分别服从Np(0,Σ),Np(0,
Σ′)多元正态分布,若Σ ≥Σ′(≥表示对应元素大小关系),则对于任意的实数t1,t2,…,tp,有
下图是使用R 软件绘出u=v 时,C(u,u)对ρ=0.5,0,-0.5 图像
1.2 尾Copula 函数.
另外一种描述尾部相依特征的是尾相依Copula 函数,下尾Copula 函数定义如下:
其中
图1 Cρ(u,u)分别对ρ=0.5,0,-0.5 图像
类似地可以定义上尾Copula 函数,在此不再赘述.当u →0 时,称)为极下尾Copula函数,而这往往是诸多学者研究的重点[3 ~4].关于极尾Copula 计算,最初由文献[5]给出相关计算定理,文献[6]在此基础上稍加改动,更加一般化,定理如下:
定理1[6]:C(u,v)为[0,1]2 连续函数,并且对于任意的u >0,v >0,有C(u,v)>0,若下列极限存在:
则
其中
事实上,依据条件可得
故
关于高斯分布的极下尾Copula 函数计算参照文献[5]方法,可得,故y)=xy,恰好渐进独立.与下极尾相关系数得到的结论相同.
2 高斯分布的极值copula 函数.
二维高斯CopulaC(u,v)的极值Copula 函数为,为计算极值Copula 函数,作以下定义:称b(w1,w2,…,wp)=为X=(X1,X2,…,XP)T尾相依函数,显然当p =1 时,b(w1)=w1,并且b(w1,w2,…,wp)为阶数为1 的齐次函数,即b(tw1,tw2,…,twp)=tb(w1,w2,…,wp),两边关于t 求导,然后令t=1 代入得到:
由尾相依函数定义可知:
令
wI 表示由wi,i ∈I 构成的集合.文献[8]证明了C(u1,u2,…,up)的极值Copula 函数:
性质3: 二元正态分布的极值Copula 函数CEV(u1,u2)=u1u2
证明:
其中上式第三个等式由于连续性以及应用罗必塔法则:
第五个等式是由)得到,故a(w1,w2)=w1+w2
极值Copula 函数CEV(u1,u2)=u1u2,表明极值变量间也是相互独立的.
参考文献:
[1] Gabriel Frahm.On the Extremal Dependence Coefficient of Multivariate[J].Statistics&Probability Letters,2006,76:1470–1481.
[2] Yin Chan and Haijun Li.Tail Dependence for Multivariate T -Copulas and Its Monotonicity[J].Insurance:Mathematics and Economics,2008,42(2):763–770.
[3] Demarta S,McNeil A J.The t Copula and Related Copulas[J].International statistical review,2005,73(1):111-129.
[4] Bortot P.Tail Dependence in Bivariate Skew-Normal and Skew-t Distributions[J].Available Online:www2.stat.unibo.it/bortot/ricerca/paper-sn-2.pdf,2010.
[5] Alessandro Juri and Mario V.Wüthrich.Tail Dependence from a Distributional Point of View[J].Extremes,2003,6(3):213-246.
[6] Elena Di Bernardino and V_eronique Maume-Deschamps.Estimating Bivariate Tail:a Copula Based Approach[J].Journal of Multivariate Analysis,2013,119:81–100.
[7] Roger B.Nelsen.An Introduction to Copulas[M].Second Edition.New York:Springer Science Business Media,2006.
[8] Aristidis K.Nikoloulopoulos and Harry Joe.Extreme Value Properties of Multivariate t Copulas[J].Extremes,2009,12(2):129-148.
因篇幅问题不能全部显示,请点此查看更多更全内容