Skip to content

模型Model

二维阿基米德 Copula

Bivariate Archimedean copula · 阿基米德联结函数

从递减生成元和伪逆构造二维均匀边缘联合分布,证明凸性充要条件,区分严格光滑密度、奇异边界和条件逆采样。

形式陈述 ​

一个二维Copula需要一整张联合CDF。阿基米德构造先选一个一维函数,再通过“变换、相加、反变换”得到这张CDF。

设 φ:(0,1]→[0,∞) 连续、严格递减,且 φ(1)=0。记

A=limu↓0φ(u)∈(0,∞],φ(0)=A.

定义伪逆 h:[0,∞]→[0,1]:在 0≤t<A 上,h(t) 是满足 φ(h(t))=t 的唯一数;若 A<∞,令 h(t)=0 对所有 t≥A 成立;另约定 h(∞)=0。候选函数为

(1)Cφ(u,v)=h(φ(u)+φ(v)),0≤u,v≤1.

二维生成元判据:式(1)是Copula,当且仅当 φ 是凸函数。此时称它为由 φ 生成的二维阿基米德Copula。A=∞ 称严格生成元;A<∞ 时,伪逆中截为零的部分不能省略。

这个判据只针对二维。仅凭递减和凸性,不能把式(1)任意改写成 h(∑i=1dφ(ui)) 后宣称它在所有维度都是Copula。

生成元也不是唯一的参数化:对任意常数 a>0,用 aφ 代替 φ,其伪逆为 h(t/a),式(1)不变。因此生成元的整体尺度不应当被当作另一个可识别依赖参数。

直觉

φ 把靠近零的累计位置变成大的非负数,把位置一变成零。把两个变换值相加,再通过递减伪逆送回零到一之间,得到两个坐标共同低于阈值的概率。

正确边缘很容易得到:第二个坐标为一时,φ(1)=0,所以 C(u,1)=h(φ(u))=u。难点在于任意矩形的概率都要非负。生成元凸性控制的正是这个二维增量,而不只是每条坐标线的单调。

伪逆中的零段有概率意义。某些低阈值组合可能根本不能同时发生;在这些地方,联合CDF应为零。若强行使用普通逆函数求值,就可能越过逆函数的定义域。

例子与边界

独立和反单调都属于这个接口 ​

取 φ(u)=−log⁡u,则 A=∞、h(t)=e−t,式(1)给 C(u,v)=uv,即独立Copula。

取 φ(u)=1−u,则 A=1、h(t)=max(1−t,0),得到

C(u,v)=max(u+v−1,0).

这正是二维Fréchet–Hoeffding下界,由 (T,1−T) 实现。它的质量集中在一条直线上,没有普通二维密度。模型的合法性不要求生成元严格凸,也不要求联合密度存在。

用同一个生成元尝试三维,得到 max(u+v+w−2,0);其在 (1/2,1]3 上的质量为 −1/2,已经非法。这给出了不能把二维凸性条件直接推广到高维的具体反例。

Clayton正参数族 ​

对 θ>0,取

φθ(u)=u−θ−1θ,hθ(t)=(1+θt)−1/θ.

其导数为 φθ′(u)=−u−θ−1<0,二阶导数为 (θ+1)u−θ−2>0,所以判据给

(2)Cθ(u,v)=(u−θ+v−θ−1)−1/θ.

在零边界用连续极限定义为零。本页限定 θ>0;不把同一个未截断幂式未经检查地用于负参数。

当 θ=1 时,

C1(u,v)=uvu+v−uv,C1(1/2,1/4)=15.

四格质量随之为 1/5,3/10,1/20,9/20,分别对应低低、低高、高低、高高;求和为一,两个低值边缘为 1/2,1/4。

Gumbel族和合法参数边界 ​

对 θ≥1,令 φθ(u)=(−log⁡u)θ。在 0<u<1 上,写 L=−log⁡u>0,则

φθ′(u)=−θLθ−1u,φθ″(u)=θu2(Lθ−1+(θ−1)Lθ−2)≥0.

于是

(3)Cθ(u,v)=exp⁡(−((−log⁡u)θ+(−log⁡v)θ)1/θ).

θ=1 给独立Copula。若 0<θ<1,二阶导数在 u 足够接近一时为负,因此不能通过本判据;递减性本身没有保证合法。

对固定正阈值,令 a=−log⁡u,b=−log⁡v。不等式 max(a,b)≤(aθ+bθ)1/θ≤21/θmax(a,b) 表明,当 θ→∞ 时式(3)趋于 exp⁡(−max(a,b))=min(u,v)。零边界也同为零,因此极限是共单调Copula;有限参数和该极限模型的尾事件仍需分别检查。

去掉凸性会出现负矩形质量 ​

取 φ(u)=1−u,它连续严格递减,却是凹函数。伪逆为 h(t)=1−t2 对 0≤t≤1 成立,超过一时为零。式(1)仍给正确边缘,但在 (15/16,1]2 上,

ΔC=C(1,1)−2C(15/16,1)+C(15/16,15/16)=1−2⋅1516+34=−18.

这是一份有限的拒绝证书;它不依赖画图判断曲率。

推论与应用

凸性如何保证全部矩形质量 ​

先记录递减反函数的一个性质:φ 凸,当且仅当其递减逆 h 在逆函数的有效范围上凸。证明例如从 x=h(s),y=h(t) 出发,凸性给

φ((1−α)x+αy)≤(1−α)s+αt.

因为逆函数递减,得到

h((1−α)s+αt)≤(1−α)h(s)+αh(t).

逆向交换两个函数重复同一论证即可。若 A<∞,将递减凸 h 在 A 以后延为零仍然凸:跨过 A 的割线斜率从非正值增加到零,不会向下跳。故可以在整个非负半轴使用其凸性。

对任意固定 a≥0,凸函数的增量 h(x+a)−h(x) 随 x 不减。这来自凸函数割线斜率的次序,也可先在分段线性逼近上按斜率求和,再由连续性取极限。

取矩形 (u0,u1]×(v0,v1],其中 u0≤u1、v0≤v1。设

s=φ(u1)+φ(v1),a=φ(u0)−φ(u1),b=φ(v0)−φ(v1).

在正坐标内部,a,b≥0,其矩形质量为

(4)h(s)−h(s+a)−h(s+b)+h(s+a+b)=(h(s+b+a)−h(s+b))−(h(s+a)−h(s))≥0.

零坐标的情况由连续极限得到。结合正确边缘与零边界,Copula盒增量判据证明式(1)合法。

必要性同样可以直接看到。若式(1)是Copula,取 0≤s<t<A,在两个坐标都使用区间 (h(t),h(s)],非负矩形质量给

h(2s)−2h(s+t)+h(2t)≥0.

因此 h 在 [0,A) 上满足中点凸性:对其中任意 a<b,代入 s=a/2,t=b/2 即可。连续性把二进权重的中点不等式扩到所有凸组合,故 h 凸;递减逆的上述性质再给 φ 凸。有限端点通过极限补齐,证明没有假设二阶导数存在。

光滑严格子类的密度 ​

进一步假设 A=∞、φ∈C2((0,1))、φ′<0。在开正方形内有

φ(C(u,v))=φ(u)+φ(v).

先对 u 再对 v 求导,得到

(5)c(u,v)=∂u∂vC(u,v)=−φ″(C(u,v))φ′(u)φ′(v)[φ′(C(u,v))]3≥0.

分母为负,前面的负号使密度非负。开正方形内所有紧矩形的概率由这个混合导数积分恢复;边界质量为零,因为边缘均匀,所以它确实给出整个联合分布的密度。这个论证不适用于前面的非严格反单调边界,不能把式(5)强套到那条奇异直线上。

条件逆采样:先核CDF,再反解 ​

对Clayton的 θ=1 情形,直接求导可得

K(u,v)=∂uC1(u,v)=v2(u+v−uv)2,0<u<1.

固定 u,它从零连续严格增加到一,且

∂vK(u,v)=2uv(u+v−uv)3>0(0<v<1).

这是一份 V 给定 U=u 的条件CDF版本,因为对任意 0≤a<b≤1,

∫abK(u,v)du=C1(b,v)−C1(a,v).

这条积分恒等式给出条件核与联合分布的相容性,并非把零概率事件 {U=u} 填进初等比值。

算法取独立均匀输入 U,Z∈(0,1),解 K(U,V)=Z。由于分母为正,取正平方根并整理得到

(6)V=UZ1−(1−U)Z.

分母大于 U>0,而输出满足 0<V<1。给定 U=u,逆变换保证输出的条件CDF为 K(u,⋅);再积分得到联合CDF C1。例如 u=1/2,z=1/4 时 v=1/3,代回恰有 K(1/2,1/3)=1/4。

每个样本只需常数次算术及一次平方根。若换成仅能求值的连续严格递增条件CDF,可以在 [0,1] 二分;要使根的包围区间长度不超过 ε>0,至多需 max{0,⌈log2⁡(1/ε)⌉} 次比较,实际函数求值和精度成本另计。可靠区间采样还须把平方根或比较误差纳入包围,普通浮点样本不是严格数值证书。

两种参数族在角落处的条件事件可以进一步用尾依赖系数区分。拟合哪一个族是额外的统计问题;生成元合法只保证存在一份概率模型,不表示它已由数据验证。

参考资料
  • Paul Embrechts、Filip Lindskog、Alexander McNeil,Modelling Dependence with Copulas and Applications to Risk Management,2001-09-10作者稿,§§6.1–6.2,印页31–33:伪逆、二维生成元判据、参数例子及密度式(6.4)。本文给出凸增量与中点必要性的完整证明,条件采样则直接反解本页指定的Clayton条件CDF。
关系图谱14 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系

使用的工具

被这些条目使用