（管理科学与工程专业论文）分位数回归理论及其应用.pdf

中文摘要 分位数回归是给定回归变量五估计响应变量y 条件分位数的一个基本方法。 它不仅可以度量回归变量在分布中心的影响，而且还可以度量在分布上尾和下尾 的影响，因此较之经典的最小二乘回归具有独特的优势。本文主要对分位数回归 的理论、c o p u l a 分位数回归、极端分位数以及分位数回归在各个领域的应用进 行了深入研究。论文的主要工作如下： 1 论文介绍了极值的基本理论，为以后的各章提供了理论基础。并选取 l o g i s t i c 分布，应用二元超阈值模型和二元点过程模型度量沪深股市收益率的尾 部相关性。结果表明：沪深股市收益率在尾部具有很强的相关性，并且这两种模 型都不失为一种很好的建模方法。 2 论文构建线性条件分位数回归模型，分析澳大利亚西部f 玎e m a n t l e 港地区 在1 8 9 7 1 9 8 9 年间年最高海平面高度与时间及年平均南方涛动指数之间的线性变 化趋势，并与经典的最小二乘回归拟合进行比较。结果表明：在不同分位数下年 最高海平面高度与时间及南方涛动指数之间所呈现的线性趋势是不同的，分位数 回归比经典的最小二乘回归能够提供更多的信息，因此对于我们进行预测和防范 具有十分重要的意义。 3 论文研究了c o p u l a 分位数回归，推导出几种常见c o p u l a 的分位数曲线， 并应用模拟研究的方法说明分位数回归估计方法的精确性。在此基础之上，选取 c l a y t o nc o p u l a ，应用c 0 p u l a 非线性分位数回归模型度量沪深股市收益率在不同 分位数下的风险相关性，并与由极值理论方法得到的结果进行比较。结果表明： 在不同分位数下沪深股市具有不同的相关关系，比普通的回归方法能更全面的描 述不同区域的风险相关关系，而极值理论方法侧重于极端情况下尾部指标的估 计。 4 论文通过研究极端分位数的估计方法及渐近性质，把极端分位数所具有 的行为特征应用到v ，i r 的研究中，建立上海股市收益率的条件分位数模型，描 述其在极端分位数下的变化趋势。并选取适当的尾部模型，在此基础之上应用外 推法预测非常极端分位数下的条件服，并与直接由分位数回归模型预测的结果 进行比较。结果表明：两种方法得到的结果变化趋势都是一致的，由外推法预测 的结果相对小一些。 关键词：分位数回归极值理论相关结构函数尾部相关性极端分位数风险 度量 a b s t r a c t ( 沁锄t i l er e g r e s s i o n i sab a s i ct o o lf 0 r e s t i m a t i i 毽c o n d i t i o n a lq u a n t i l e so fa r e 8 p o n s ev 撕a b l eyg i v e nav e c t o ro fr e g r e s s o r sz i tc a nb eu s e dt om e a s u r en l ee 虢c t o fr e 伊e s s o r sn o to i l l yi nm ec e n t e ro f ad i s t r i 鼬i o 玛b u ta l s oi nt h eu p p e ra n dl o w e r t a i l & s oi th a sm u c hm o r ea d v 锄t a g e s 也a n 也ec l a s s i c a ll e a s ts q l l a r er e g r e s s i o n t h e t l l 巧o fq 删1 er e g r e s s i o 玛c o ”l aq u a n t i l er e g r e s s i o 玛c ) 【仃e i n a lq u a ：蚯l e s 锄d 印p l i c a t i o 璐o fq u a m i l er e g r e s s i o ni nn l l n yf i e l d sa r ed i s c u s s e di i l 廿l i sp a p 既t h e n l a i na c m e v e m e n t so f l i sw o f ka r el i s t e da sf o l l o w s ： 1 n l eb a s i ce x 缸e m ev a l u em e o 巧i si n l d u c e d ，w 1 1 i c hi s 也eb a s i so f0 廿1 e r c h 印衄- s w bc h 0 0 s el o g i s t i cd i s 臼曲u 如na n du 8 eb i v 撕a t ee x c e s st l l r e s h o l dm o d e l 趾db i v a r i a t ep o mp r o c e s sm o d e lt 0m e a s u r et l l e 僦d 印e r l d e n c eo fs h a n 班1 a i 柚d s h e n 出e i ls t o c km a r k e t n l er e 8 u n ss h o wm a t 也e 咖r a t e so fs h a n 曲a ia n d s h e m 出e ns t o c km a r k e t sk l v es 仃o n gt a i ld 印e n d e n c ea i l d 也et w om o d e l sa r ee x c e l l e n t f o ra p p l i c a t i o n 2 硼1 e1 i i l e a r 仃e n do ft h ea n 玎u a lm a x i m u ms e al e v c la tf r e m a n t l ep 嘣w e s t e r n a u s 僦i a ，r e l a t e dw i mt i i n e 觚ds o u 廿1 e mo s c i a t i o n 砌e xd 证n g18 9 7 19 8 9i s a n a l y z e db yl i n e a rc o n m t i o n a lq 啪t i l er e g r e s s i o n 皿1 0 d e l a n d 蚀【er e s u l ti sc o m p a r e d w i t l l m to ft h ec l a s s i c a ll e a s ts q u a r er e g r e s s i o n t h er e l 湖t ss h o wt l l a t ，硼【d c rd i 侬黔m q u a n t i l e s ，也el i n e a r 慨do fm e 锄m l a lm a ) 【i n m ms e al e v e lr e l a t e dw 砒lt i i n ea n d s o u 廿1 锄o s c i l l 撕o ni 】m e xi sd i 妇陪f e m ，锄dq u a m i l er e g r e s s i o nc a np r o 、，i d em u c h m o r ei n f o n l 豫t i o n l a nt 1 1 ec l a s s i c a ll e a s t 双l u a r e r e g r e s s i o n s o “i so f 目e a t s i g n i f i c a r l c ef o rp r e d i c t i o na i l dp r e v 咖i o n 3 t h e 也e 0 巧o fc o p u l aq 啪t i l er e g r e s s i o ni ss t i j d i e d 锄dt 1 1 eq u a 矾1 ec u n ，e so f s e v e m lc o m m o nc o p u l a sa r eo b t a h l e d t h ea c c u r a c yo fq u a i i t i l er e 母e s s i o ne s t i 玎1 a t i o n i ss h o w nb ys i i i l _ u l a t i o nr c s e a r c h w ec h o o s ec l a y t o nc o p u l a 锄du s ec o p u l an o n l i i l e a r c o n d i t i o i l a lq u a m i l er e g r e s s i o nm o d e lt 0m e a s u r et h et a i la r e ad s kd 印e n d e l l c ei i l s h 锄曲a i 锄ds h e n z h e i ls t o c km a r k e t s 气n dt 1 1 e nt l l er e s u ho fn l ：i sa p p r o a c hi s c o m p a r e dw i mt h et a i ld 印e 1 1 d e l l c em e a 鲫r eb ye x 愉n ev a l u em e t h o d t h er e s u n s s h o wm a ts h a i l 曲a ia n ds h e i l 疝e ns t o c km a r k e t sh a v ed i 鼬r i s kd 印e n d e n c eu n d e r d i a e r c n tq u 枷l e sa n de x 仃e m ev a l u et h e o 巧m e t h o do i 蚵f o c u s e so nn l ees t ：i m a t i o no f t a i ld 印e n d a l c e 4 b y 咖d y i i 培t 1 1 ee s t i m 撕0 nm e 1 0 d 锄d 蛔，n 】p 跏cb e h a 访o r 8o fe x 仃c 柚a 1 q u 觚t i l e s ，w ea p p l yi t sb e h a v i 邮t 0 血er e s e a r c ho fv 撼t h ec o n d i t i o m l 畔嘣l e r c g r e s s i o nm o d e lo fr e t u l mr a t e so fs h a n g h a i s t o c km r k e ti se 鼬1 i s h e d ，w l l i c h d e s c r i b e st h e 臼e n do fr a t e su n d e re x 仃e r n a lq 啪t i l e s c o n d i t i o n a l 、a ri nv a ye x 胁e q 啪t i l e si sp r e d i c a t e db yu s i n ge x 仃a p o l a t i o nm e m o d su n d e r 山ep r o p e rt a i lm o d e l c o m p a r i s o nw i mt h ep r e d i c t i o no f 也eo r d i n a d rq u 枷l er e 黟e s s i o nm o d e li sa l s 0 g i v e n t h er e 团n t ss h o wt h a t 也et e n d e n c i e so ft l l e 铆op r e d i c d o n sa r es i i i l i l a ra n dt h e v a l u ee s t i l l l 砷e db y 也ee x 订a p o l a t i o nm e t l l o d si sr e l a 吐v e l ys m a u k e y w o i m s ：q 砌咀t i l er e g r e s s i 吗e x 仃锄ev a l u e 也e o c 叩u l a ，伽d e p e n d e n c e ， e 炳伽1 a lq 啪t i l e s ，v - a r 独创性声明 本人声明所呈交的学位论文是本人在导师指导下进行的研究工作和取得的 研究成果，除了文中特别加以标注和致谢之处外，论文中不包含其他人已经发表 或撰写过的研究成果，也不包含为获得丕鲞盘堂或其他教育机构的学位或证 书而使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中 作了明确的说明并表示了谢意。 躲羔够期：加嗡年m 学位论文版权使用授权书 本学位论文作者完全了解苤盗盘堂有关保留、使用学位论文的规定。 特授权墨鲞基堂可以将学位论文的全部或部分内容编入有关数据库进行检 索，并采用影印、缩印或扫描等复制手段保存、汇编以供查阅和借阅。同意学校 向幽家有关部门或机构送交论文的复印件和磁盘。 ( 保密的学位论文在解密后适用本授权说明) 日 导师签名：叟勘产 、 f 签字同期：枷口分年6 月y 日 天津大学博士学位论文 1 1 论文研究背景 第一章绪论 1 1 1 分位数回归理论的演进 1 8 7 0 年，英国的高尔顿在研究人类身高的遗传规律时发现：父母是高个子的， 其子女的身高有低于父母身高的趋势；相反，父母是矮个子的，其子女的身高却往 往有高于父母身高的趋势。从全局来看，高、矮个子人的子女都有“回归”于一 般人身高的趋势，这就是统计学上“回归 的最初涵义。1 8 8 6 年，高尔顿正式提 出了“回归的概念。经过他的学生皮尔逊多年的进一步发展，这个出于生物统 计学领域的概念，便被推广为一般统计方法的重要概念。 回归分析具有悠久的历史，其中以经典的最小二乘回归应用最为广泛。原因 是最小二乘回归的解释与人们的直观想象一致，并且该方法易于计算，特别当假 设误差是正态分布时，它具有无偏性、有效性等优良性质。但是在实际问题中，应 用最小二乘回归需要满足较严格的条件，如等方差性，随机误差间两两不相关等 条件；当需要对回归系数的显著性进行推断时，通常还要假定随机误差的正态性 等；当数据分布是重尾或有异常值时，结果的稳健性较差。特别是对于大量的数 据，应用最小二乘回归法只能得到一条回归曲线，而一条曲线所能提供的信息毕 竟是有限的。因此，人们在使用经典的回归分析同时，也在不断地探索更新更好 的方法。 1 9 7 8 年，k o 即k e r 和b a s s e 仕【1 】首次提出了分位数回归的概念，他们是在1 8 1 8 年l a p l a c e 【2 】提出的中位数回归( 最小绝对偏差估计) 理论的基础上，把中位数回 归推广到了一般的分位数回归。分位数回归相对于最小二乘回归应用的条件更为 宽松，挖掘的信息更丰富，它不仅可以度量回归变量在分布中心的影响，而且还 可以度量回归变量在分布上尾和下尾的影响，即捕捉整个条件分布的特征。特别 当误差为非正态分布时，分位数回归估计量比最小二乘估计量更有效，因此它弥 补了最小二乘回归的不足。 分位数回归理论在近二、三十年的发展过程中，k o e i l l ( e r ，b a s s 甜，p o w e l l ， c h e m o z h u k o v 等都作出了非常重大的贡献。1 9 7 8 年，k d e n k e r 和b 勰s 鲋【3 】提出分 位数回归概念，把响应变量看作是其它变量的线性函数，推导出分位数回归系数 的渐近分布，从而发展了线性分位数回归理论；1 9 8 2 年，他们又研究了分位数 第一章绪论 回归的线性假设检验以及异方差的稳健性检验，为分位数回归的应用提供了保证 刚【5 】；1 9 8 6 年，b a s s e n 等【6 】研究了回归分位数的强相合性等性质；随着线性分位 数回归的发展，同年p o w e l l 【7 】基于删失模型提出了非线性分位数回归；鼬m 和 w h i t e 【8 】研究了非线性分位数回归估计量的一致性等性质；l9 8 7 年，k o e n k e r 等【9 】 提出了线性模型的l 估计法，并且，同年k o e l l l 【e r 等【l0 】提出了关于分位数回归的 有效算法；b u c b 【l l 】【1 2 】【1 3 】在1 9 9 5 年和1 9 9 8 年分别讨论了分位数回归模型渐 近协方差矩阵的估计方法，以及分位数回归最新的一些发展，并应用它分析了美 国女性薪水结构的变化情况；2 0 0 0 年，k o e n l 【e r 和z h i i i ex i a 0 【1 4 懈决分位数回归 过程中存在的特定推断问题；同年飚m 和m u l l e r 【1 5 】关于两步分位数回归的渐近特 性进行了研究；2 0 0 1 年，t 勰c h e 【1 6 j 研究了分位数回归的无偏性；2 0 0 2 年，k 0 e n k 饼 等【1 7 】又讨论了线性异方差模型的l 估计法。 在此基础之上，c h e m o z h u k o v 和h a nh o n g 【1 8 】提出了研究删失分位数回归的 三步评估方法；吴建南和b r e t s c l l i l e i d e r 等【1 9 】用蒙特卡罗( m o m ec a l l o ) 方法产生 1 0 0 个随机数据集合来比较显著加权分析方法与分位数回归的优劣；k 0 t t a s 和 l o ) 和 瓯) ，使得 h m p “! 鱼吐工) ：日( x ) 工r ( 2 2 ) “” 成立，其中日0 ) 是非退化的分布函数，那么日必属于下列三种类型之一： i 型分布：日l ( 力= 唧 一，) ， 啊 o ( 2 3 ) 盯仃 其中，善r ，盯 o ，称作位置参数( 1 0 c a t i o n p a 阳me t i 神，仃称作尺度参数( s c a l e p a r a 加e t e r ) ，孝称作形状参数( s h a p ep 踟e t 神。日称作广义极值分布( g e i l 洲i z e d e x 臼伽n ev a l u ed i s 劬而。吣，简记为g e v 分布。 形状参数的取值决定了极值分布的类型： ( 1 ) 当善= o 时，日，仃，善) 表示极值i 型分布，位置参数和尺度参数不变； ( 2 ) 当孝 o 时，取口= 1 善，日心声，仃，善) 表示极值n 型分布，位置参数为丘一阳， 尺度参数为阳； ( 3 ) 当毒 o ) 和溉) ，使得 f ”( 口。x + 屯) = f ( 工) ( 2 - 4 ) 则称分布函数f ( x ) 是最大值稳定的( m a x s t a b l e ) 。 由此可证明三种极值分布都是最大值稳定分布，并且由定理2 1 和定义2 1 可以进一步得到结论：一个分布函数，( 工) 是最大值稳定分布，当且仅当f ( x ) 是三 种极值分布类型之一。 2 1 2 极值分布分位数与重现水平 定义2 2 称分布函数f ( x ) 的反函数 f 。1 ( p ) ；i 1 1 f ( 工r ：f ( x ) p ) ， 0 p l( 2 5 ) 第二章极值统计理论 为它的分位数函数，而毛= f 一( p ) 称为分布函数f 的p 分位数。 由定义2 1 可知，广义极值分布的p ( o o 和蛾 ，使得 y 一 l i m p r c 尘l 兰x ) = l i m p ( 吒x + 玩) = 日( 功 4 ” 口n ”一 成立，则称随机变量x ( 或x 的分布函数f ) 属于极值分布日( 力的最大值吸引场 ( m a ) 【i n l 啪d o m a i l l o f a t 仃a c t i o n ) 。记作x 四h ( 日) 或f 坦功( 奶。 在研究极值分布的最大值吸引场问题时，v 0 nm i s e s 函数起了重要作用。 定义2 4 设f 是分布函数，上端点，o o ，若存在某个z ，使得f 可以表 示成如下形式 一 1 户( 功= c e x p - 1 4 三j 疵) z o 为尺度参数，f e 尺为形状参数。 类似于广义极值分布，形状参数有时用口来表示，令口= l 善，得到广义p a r e t 0 分布的另一种表达形式： f - 兰 g l 似，砷= 1 一p 9 ， 工， 【o ， 工 o ； i o ，x 0 x p g ，g 2 ，g 3 分别称为p a r e t oi 型、型和型分布。当，础，口= 1 时，称为标 准g p d ，相应的分布简记为g 0 ；孝) 或g l ( 功，g ，口) 和g 3 似口) 。图2 - 1 给出了当 口= o 5 ，口= 1 和口= 2 时，标准p 锄加i 型、i i 型和i i i 型分布的密度函数曲线。 图2 - l 标准g p d 密度函数图 下面我们讨论广义p a r e t o 分布的一个重要性质- p o t ( p e a ko v e rn l i e s h o l d ) 稳定性，即某一分布，对于给定阈值的超出量经过规范化后的渐近分布仍然是它 本身，下面给出严格的定义及性质。 定义2 6 对于给定的分布函数，( 力，如果存在常数气，包，使得对任何实数 x ，都有 e ( 吼x + 吃) = 以力， 其中层( z ) = p x 一”x 陋 “) 是( 2 8 ) 式所定义的超出量分布函数，则称分布函数 ，( 功具有p o t 稳定性( p o t 鼬i l i 够) ，或称分布函数f ( 工) 是p o t 稳定分布。 性质2 1 广义p a r e t 0 分布是p o t 稳定分布，常数气，吃可取为： 气：1 + 孝兰兰吃= 一( 1 + 孝兰坐) 天津大学博士学位论文 定义2 7 设随机变量彳的分布函数f ( z ) ，为尸( x ) 支撑的上端点，石超过阂 值“的超出量分布为e ( x ) ，如果存在广义p a r e t o 分布g ( 工) ，使得 e ( 功= g ( 劝， “ 则称石( 或分布函数以工) ) 属于广义p a r e t o 分布g ( 工) 的p o t 吸引场。一 广义极值分布和广义p a r e t o 分布之间存在着非常密切的关系。 性质2 2 广义极值分布属于广义p a r e t o 分布的p o t 吸引场。 性质2 3 广义p a r e t o 分布属于广义极值分布的最大值吸引场。 2 2 多元极值理论 对于某些比较简单的问题，应用一元极值理论进行分析就足够了，然而实际 中许多问题依赖于多个因素的影响，这就迫切需要多元极值理论的发展。多元极 值理论目前已成为极值理论中较新且发展迅速的领域，并且形成了区组最大值模 型、超阈值模型以及点过程模型等，多元极值分布广泛地应用于气象、水文、环 境、工程、金融等领域。为了简便，我们只讨论二元极值模型，在研究二元极值 分布时，若变量之间是相互独立的，则可以直接把二元极值问题转化为一元极值 问题；若变量之间存在着某种相关性，则我们引入相关结构函数c 0 p u l a 来描述 它们之间的相关关系。 2 2 1 相关结构函数及其性质 c o p u l a 描述了多元随机变量的一维边缘分布与联合分布之间的函数关系， 其严格的数学定义如下： 定义2 8 如果定义在集合 0 ，1 】 0 ，1 上的一个二元函数c 满足以下条件 ( 1 ) c 的定义域为 o ，1 】 0 ，1 】； ( 2 ) c 满足边界条件：c ( “，0 ) = c ( 0 ，) = o ，c ( “，1 ) = “，c ( 1 ，1 ，) = v ； ( 3 ) c 是二增( 2 面c r e a s i i l g ) 函数，即对任意“l ，”2 ，l ，吃 o ，1 】，且“l “2 ， y l 1 ，2 ，有： c ( “2 ，吃) 一c ( “2 ，) 一c ( “l ，v 2 ) + c ( “l ，v 1 ) o 则称函数c 为c o p u l a 。 定理2 2 ( s k l a r 定型7 伽) 设随机向量( x ，d 的联合分布函数为f ( x ，y ) ，边缘分 布函数分别为e ( x ) 和最( 少) ，则对于任意的( 工，y ) r 2 ( r 表示实数集合) ，一定 存在一个c o p u l ac ，使得 第二章极值统计理论 f ( 五少) = c ( 鼻( x ) ，e ( y ) ) ( 2 - 1 1 ) 若e ( x ) ，e ( y ) 都是连续分布函数，则c 是唯一的；否则，c 在啪呲上 是唯一的。反之，如果c 是一个相关结构函数，e 和最是一元分布函数，则由 ( 2 1 1 ) 式定义的函数，也y ) 是一个边缘分布为互和e 的二元联合分布函数。 定理2 2 是c o p u l a 的存在性定理，同时也是已知联合分布，后，求c o p u l a c 的理论依据。令u = 互( x ) ，y = 最( 】，) ，它们都是服从区间 0 ，1 】上的均匀分布， 则( 【厂，y ) 的联合分布函数为 c ( “，v ) = p r ( 【，“，y 1 ，) = f ( e 叫 ) ，e 1 ( v ) )( 2 - 1 2 ) 所以c o p u l ac 可以看作是边缘分布为区间 0 ，1 】上均匀分布随机向量( u ，矿) 的联 合分布函数。 因为c o p u l a 实质上是一个二元分布函数，如果对应的概率密度函数c 存在， 则c 可以由下式得到： c m ，y ) ：丝丝堕 a z 咖 利用c 0 p u l a 与联合分布函数只边缘分布函数e ，最的关系式( 2 - 1 1 ) ，得( x ，d 的联合密度函数为 ，( x ，j ，) = c ( 互( 功，e ) ) 石( 曲以( y ) ， ( 2 1 3 ) 其中z ( x ) ，五( y ) 分别表示随机变量x ，】，的概率密度函数。 c 0 p u l ac 具有如下基本性质： 性质2 4 【7 1 】对于任意( “，功 0 ，1 】2 ，都有形( “，d c ( “，d 鲋( “，d ，其中 m ，功= m i l l 缸，嵋，矽 ，d = m a x 劬+ 1 ，一1 ) 都是c 0 p u l a ，称m ，v ) 和形( “，1 ，) 分别为f r 6 c h e t - 。h o e 尉i i 坞上、下界。 特别地当u ，y 是区间【0 ，1 】上服从均匀分布的随机变量时，若它们的联合分布 函数是m ，则p r ( u = v ) = 1 ，称u ，矿完全正相关( 主对角线) ；若它们的联合分 布函数是形，则p 州+ v = 1 ) = l ，称u ，y 完全负相关( 次对角线) 。完全正相关 和完全负相关分别对应于一个变量是另一个变量的单调递增与单调递减函数。 性质2 5 设连续随机向量( x ，】，) 的c o p u l a 为o - y ，则石与】，相互独立当且 仅当o y ( “，v ) = 删，并称该c o p u l a 为乘积c 0 p u l a ，记作n 。 性质2 6 设连续随机向量( x ，聊的c o p u l a 为q 1 ，则口( x ) 与( y ) 都是严 格单调递增函数，则 c z ，r ( h ，叻= q ( j x ，( y ) ( m ，v ) 天津大学博士学位论文 即在严格单增变化下，随机变量间的c o p u l a 保持不变。如果口( x ) ，厦y ) 都是严 格单调递减函数，则 c 口( 工) ，( y ) ( “，v ) = “+ v l + c x ，r ( 1 一“1 一v ) ， 如果口( x ) 严格单调递增，( 】，) 严格单调递减，则 q ( x ) 。m ( “，v ) = 掰一c 石，y ( “，1 一v ) ， 如果口( 彳) 严格单调递减，觑y ) 严格单调递增，则 巴( j ) ，艄( 甜，功= v g ，y ( 1 一“，y ) ， 性质2 8 令c 是一个c o 叫a ，则对任意的s l ，s 2 ，f 【0 ，1 ，5 l j 2 ，有 c ( j l ，f ) c ( s 2 ，f ) ，c ( f ，s 1 ) c ( f ，s 2 ) 即关于每个变量非降。 对任意“i ，“2 ，1 ，l ，吃 o ，1 ，且“l “2 ，h 吃，有 l c ( “2 ，屹) 一c ( “l ，1 ，1 ) l i “2 一“li + l 吃一h 即满足l i p s c l l i t z 条件。 性质2 9 令c 是一个c o p u l a ，则c 关于掰，的偏导数存在，且 o 导c ( “，v ) 1 ，o 导c ( “，) 1 口“c 性质2 1 0 令q ，c 2 是两个c o p u l a ，口 0 ，1 是权重系数，则( 1 一口) c l + 仍 是一个c o p u l a ，即c o p u l a 的线性组合仍是一个c o p u l a 。 定义2 9 设随机向量( z ，y ) 的联合分布函数为f 瓴y ) ，对应的c o p u l a 为c ， 则称联合生存函数f ( 五y ) = p r ( x 五】， j ，) 对应的c 0 p u l a 为生存c o p u l a ( 翻l n ，i v a lc o p u l a ) ，记作c ，( 毛y ) = c ( 互( x ) ，最( y ) ) ( 2 1 5 ) 生存c o p u l a 也是c o p u l a ，这是因为仑( f ，o ) = e ( o ，f ) = f + o l + c ( 1 一f ，1 ) = o ， c ( f ，1 ) = c ( 1 ，f ) = f + 1 一l + c ( 1 一f ，o ) = f ，满足定义3 1 的边界条件。对任意的 “l ，”2 ，1 ，l ，吃 o ，l 】，上l “l “2 ，v 1 屹，有 c ( “2 ，屹) 一c ( “2 ，m ) 一c ( “l ，屹) + c ( “l ，m ) = c ( 1 一“2 ，1 一吃) 一c ( 1 一“2 ，l v 1 ) 一c ( 1 一，1 一屹) + c ( 1 一“l ，l v 1 ) 0 满足2 增性。同c o p u l a 是连接边缘分布函数与联合分布函数的桥梁一样，生存 c 0 p u l a 是连接边缘生存分布与联合生存分布的桥梁。 生存c o 叫a 与原c o p m a 之间的关系为： 第二章极值统计理论 c ( “，) = “+ ，一1 + c ( 1 一“，l d( 2 一1 6 ) 定义2 1 0 设c 。，c 2 是两个c o p u l a ，如果对于任意的 ，v ) 0 ，1 】2 ，都有 q ( “，d q ( “，1 ，) 则称c l 小于c 2 ，记为c l q ，“ 一称为相关序。 根据c o p u l a 的性质3 1 ，对于任意c o p u l a ，都有矿 c m ，说明f r 6 c h e t - h o e 衄i i 培下界矽是所有相关结构函数中最小的，而上界m 是最大的c o p u l a 。 2 2 2 常见c o p u l a 族 在这一节主要介绍椭圆型c o p u l a 以及阿基米德c o p u l a 。 椭圆型c o p u l a 1 二元正态c 叩u l a c 似v ；加p p 赤e 叶高舡 倍 其中m - 1 ( ) 为标准正态分布的反函数，p 为z ，】，之间的线性相关系数，且 一1 p 1 。由于金融数据大多呈现尖峰厚尾，而二元正态c 叩u l a 具有对称性， 两尾比较薄，因此无法捕捉金融市场之间的相关关系。 2 二元f c o p u l a 。 c “唧= p p 南 l + 筹) 导蛐 ( 2 - 1 8 ) 其中f - 1 ( ) 表示自由度为d 的f 分布的逆函数，p 为疋】，之间的线性相关系数， 且一l p 1 。类似于二元正态c 叩u l a ，二元fc o p u l a 也具有对称性特点，因此 这两个c o p u l a 都无法捕捉到变量间的非对称的相关关系。 在金融风险分析中，由于二元正态c o p u l a 尾部是渐近独立的，如果把相关 变量用正态c o p u l a 描述，则会降低风险。有时出于计算上方便的考虑，用fc 0 p u l a 描述两个变量的相关性，尽管它也具有对称性，但是它所描述的是随机变量的尾 部渐近相关性，在风险度量上，比正态c o p u l a 有所改进。 阿基米德c o p u l a 阿基米德c o p u l a 是c o p u l a 族中最重要的一类。它具有许多特点：( 1 ) 构造 非常简单( 可以由生成元矿来构造) ；( 2 ) 许多c o p u l a 族都属于这类c o p u l a ；( 3 ) 这类c o p u l a 具有许多特点，如联合非对称性，厚尾性等。由于阿基米德c o p u l a 具有建模需要的良好性质，它已在越来越多的模型中被应用。 天津大学博士学位论文 定义2 1 1 若连续严格单调递减函数缈： o ，1 】专 0 ，叫满足伊( 1 ) = 0 ，则伊的拟 逆函数q s e u d o i l w e r s c ) 矿_ 1 1 定义为 扩= 妒：舌竺竺 并且- 1 1 的定义域d d m ( 9 q 1 ) = 【o ，叫，值域尺册( 伊卜1 1 ) = o ，1 】。 所有的二元阿基米德c 0 p u l a 都可以写成下列形式 c ( ，“2 ) = 伊卜u 矿( “i ) + 烈“2 ) ，“l ，“2 o ，1 ( 2 1 9 ) 其中函数缈( ) 称为生成元( g e i l e r a _ t i d r ) ，满足以下三个条件：( 1 ) 在【0 ，1 】上连续，严 减；( 2 ) 矿( 1 ) = o ；( 3 ) 下凸的，即妒1 ( ) o 。如果驴( 0 ) = ，称伊为严格生成元， 有矿卜1 1 = 矿。每个阿基米德c o p u l a 都有一个独有的生成元。当9 ( o ) 时，伊为 非严格生成元。由于非严格生成元在建模中得到的结果不是很明确，我们仅考虑 严格生成元情形，此时有 尹( c ( m ， ，) ) = 矿( “) + 缈( 1 ，) 因此，阿基米德c 0 p u l a 可表示为 c ，力= 尹- 1 ( 矽( f ) + 烈v ) ) 称为严格阿基米德c 0 p u l a 。 下面从相关性角度讨论几种常用的阿基米德c 0 p u l a ： 1 g u m b e lc o p u l a c ，v ；们= e x p 一( ( 一l n “) p + ( 一l n l ，) 9 ) 口)( 2 2 0 ) 其中1 p 佃，当9 = 1 时，阢y 相互独立；当秒一佃时，阢y 完全正相关。 g u i 】n b e lc o p u l a 具有非对称性，它对于变量在分布上尾处的变化十分敏感，因此 能够很好的捕捉金融市场中上尾相关性的变化。 2 c l a y t o nc o p u l a c ( “，v ；口) = ( 材。一+ 厂9 1 ) 一1 7 护( 2 2 1 ) 其中o 目 o 时，u ，矿正相关；当秒 0 时，【，矿负相关；当 9 一。时，c ( “，v ；国= 圳，表明两变量相互独立。f r a n kc 0 p u l a 函数的性质类似 第二章极值统计理论 于正态c o p u l a ，其分布的尾部渐近独立，且具有对称相关性，对于变量间尾部 的相关变化不太敏感。 4 c 1 a y t o n j o ec 叩u l a c ( “，v ；口，回= 1 一( 1 一 ( 1 一( 1 一“) 占) 前+ ( 1 一( 1 一功占) 埘一1 ) 】一1 坩) 1 借( 2 2 3 ) 其中o p “) = 乃 ( 2 - 2 4 ) 存在且非零( 零) ，则称x 和y 上尾渐近相关( 上尾渐近独立) ；若极限 厶 l 姆p “f ( x ) ”j g ( y ) 0 ，y o ，( 2 - 2 9 ) 其中 m = f 一号，半煳毗 ( 2 - 3 0 ) 矿为指数测度，是一l 阶齐次函数【1 唧，满足 y ( 甄缈) = 口- 1 y ( 五y ) ， 墨笋是 o ，1 上均值为圭的分布函数，满足 1 = f 础( 纠= f ( 1 一国) 承( 功 ( 2 3 1 ) 并且称日( x ，少) 为二元极值分布。 二元极值分布与区间 0 ，1 】上满足( 2 - 3 1 ) 式的分布函数k 2 一一对应，由于k 不具有有限参数形式，因此二元极值分布不能用有限参数形式表示。另外，k 可 以可微，也可以不可微。当k 可微时，设密度为七，有 y ) ：f 一 詈，半撇训戗 同样矿o ) 也可以写作 啪尺* ) f t 嚣，等拟州国 令 天津大学博士学位论文 彳 ) = f m a x 国( 1 - f ) ，( 1 一国弦) d k ( 国) ( 2 3 2 ) 其中4 ( f ) 称作p i c k a i l d s 相关函数，满足彳( 0 ) = 彳( 1 ) = l ，且m a x f ，1 一f ) 彳( f ) 1 ， 0 s f 1 ，贝0 有 矿( 戈，力；( 三+ 与么内 x yx 七y 因此二元极值分布也可以写作 j 了( 工，少) = c ) 【p 一二+ 土) 4 ( ! ! 一) ) ( 2 3 3 ) y 工一