形式陈述
一个二维Copula理路Copula 与 Sklar 定理Sklar's theorem · Copula · 联结函数用均匀边缘的联合CDF分离边缘与依赖,证明一般Sklar表示、连续边缘下唯一性,以及离散边缘的不同延拓。需要一整张联合CDF。阿基米德构造先选一个一维函数,再通过“变换、相加、反变换”得到这张CDF。
设 连续、严格递减,且 。记
定义伪逆 :在 上, 是满足 的唯一数;若 ,令 对所有 成立;另约定 。候选函数为
二维生成元判据:式(1)是Copula,当且仅当 是凸函数。此时称它为由 生成的二维阿基米德Copula。 称严格生成元; 时,伪逆中截为零的部分不能省略。
这个判据只针对二维。仅凭递减和凸性,不能把式(1)任意改写成 后宣称它在所有维度都是Copula。
生成元也不是唯一的参数化:对任意常数 ,用 代替 ,其伪逆为 ,式(1)不变。因此生成元的整体尺度不应当被当作另一个可识别依赖参数。
直觉
把靠近零的累计位置变成大的非负数,把位置一变成零。把两个变换值相加,再通过递减伪逆送回零到一之间,得到两个坐标共同低于阈值的概率。
正确边缘很容易得到:第二个坐标为一时,,所以 。难点在于任意矩形的概率都要非负。生成元凸性控制的正是这个二维增量,而不只是每条坐标线的单调。
伪逆中的零段有概率意义。某些低阈值组合可能根本不能同时发生;在这些地方,联合CDF应为零。若强行使用普通逆函数求值,就可能越过逆函数的定义域。
例子与边界
独立和反单调都属于这个接口
取 ,则 、,式(1)给 ,即独立Copula。
取 ,则 、,得到
这正是二维Fréchet–Hoeffding下界理路Fréchet–Hoeffding 界Fréchet–Hoeffding bounds · Fréchet界固定各坐标边缘后,对联合阈值概率给出逐点紧界,区分二维可同时取等与高维下界函数不再是CDF。,由 实现。它的质量集中在一条直线上,没有普通二维密度。模型的合法性不要求生成元严格凸,也不要求联合密度存在。
用同一个生成元尝试三维,得到 ;其在 上的质量为 ,已经非法。这给出了不能把二维凸性条件直接推广到高维的具体反例。
Clayton正参数族
对 ,取
其导数为 ,二阶导数为 ,所以判据给
在零边界用连续极限定义为零。本页限定 ;不把同一个未截断幂式未经检查地用于负参数。
当 时,
四格质量随之为 ,分别对应低低、低高、高低、高高;求和为一,两个低值边缘为 。
Gumbel族和合法参数边界
对 ,令 。在 上,写 ,则
于是
给独立Copula。若 ,二阶导数在 足够接近一时为负,因此不能通过本判据;递减性本身没有保证合法。
对固定正阈值,令 。不等式 表明,当 时式(3)趋于 。零边界也同为零,因此极限是共单调Copula;有限参数和该极限模型的尾事件仍需分别检查。
去掉凸性会出现负矩形质量
取 ,它连续严格递减,却是凹函数。伪逆为 对 成立,超过一时为零。式(1)仍给正确边缘,但在 上,
这是一份有限的拒绝证书;它不依赖画图判断曲率。
推论与应用
凸性如何保证全部矩形质量
先记录递减反函数的一个性质: 凸,当且仅当其递减逆 在逆函数的有效范围上凸。证明例如从 出发,凸性给
因为逆函数递减,得到
逆向交换两个函数重复同一论证即可。若 ,将递减凸 在 以后延为零仍然凸:跨过 的割线斜率从非正值增加到零,不会向下跳。故可以在整个非负半轴使用其凸性。
对任意固定 ,凸函数理路凸函数Convex function函数在任意凸组合处不超过相同权重下函数值的凸组合。的增量 随 不减。这来自凸函数割线斜率的次序,也可先在分段线性逼近上按斜率求和,再由连续性取极限。
取矩形 ,其中 、。设
在正坐标内部,,其矩形质量为
零坐标的情况由连续极限得到。结合正确边缘与零边界,Copula盒增量判据证明式(1)合法。
必要性同样可以直接看到。若式(1)是Copula,取 ,在两个坐标都使用区间 ,非负矩形质量给
因此 在 上满足中点凸性:对其中任意 ,代入 即可。连续性把二进权重的中点不等式扩到所有凸组合,故 凸;递减逆的上述性质再给 凸。有限端点通过极限补齐,证明没有假设二阶导数存在。
光滑严格子类的密度
进一步假设 、、。在开正方形内有
先对 再对 求导,得到
分母为负,前面的负号使密度非负。开正方形内所有紧矩形的概率由这个混合导数积分恢复;边界质量为零,因为边缘均匀,所以它确实给出整个联合分布的密度。这个论证不适用于前面的非严格反单调边界,不能把式(5)强套到那条奇异直线上。
条件逆采样:先核CDF,再反解
对Clayton的 情形,直接求导可得
固定 ,它从零连续严格增加到一,且
这是一份 给定 的条件CDF理路条件分布Conditional distribution · Regular conditional distribution给定观测值后随机量的概率律,以及它与原联合分布相容的核表示。版本,因为对任意 ,
这条积分恒等式给出条件核与联合分布的相容性,并非把零概率事件 填进初等比值。
算法取独立均匀输入 ,解 。由于分母为正,取正平方根并整理得到
分母大于 ,而输出满足 。给定 ,逆变换保证输出的条件CDF为 ;再积分得到联合CDF 。例如 时 ,代回恰有 。
每个样本只需常数次算术及一次平方根。若换成仅能求值的连续严格递增条件CDF,可以在 二分;要使根的包围区间长度不超过 ,至多需 次比较,实际函数求值和精度成本另计。可靠区间采样还须把平方根或比较误差纳入包围,普通浮点样本不是严格数值证书。
两种参数族在角落处的条件事件可以进一步用尾依赖系数理路尾依赖系数Tail dependence coefficient · 上尾与下尾依赖用同分位阈值下的条件事件极限衡量共同极端,计算上下尾例子,并展示零尾依赖、零相关和极限不存在的不同边界。区分。拟合哪一个族是额外的统计问题;生成元合法只保证存在一份概率模型,不表示它已由数据验证。
参考资料