||
为什么 Copula 的 Kendall 复合还是 Copula
这篇文章说的"复合",专指 Kendall 复合——先用 Kendall 变换把每一块 copula 的输出拉成均匀,再在纤维上做概率混合。
之所以要专门给它一个名字,是因为普通 Sklar 复合 $C(D_1,\dots,D_n)$ 并不可信:它本质上只是概率分布函数的直接代入,没有任何概率保证。这种代入可能出现负密度,也可能积分超过 $1$,还会破坏边缘均匀性,这些都有反例。普通复合只是一个普通函数,不是概率对象。
而 Kendall 复合拼出来的东西,仍然是 copula。这件事不依赖绝对连续性,是 copula 的公理加上 disintegration 定理直接保证的。
一、普通 Sklar 复合为什么不够
普通 Sklar 复合写成
$$C(D_1(v^{(1)}),\dots,D_n(v^{(n)})),$$
外层 $C$ 吃 $n$ 个输入,每个输入是某个 $D_k$ 的输出。它的做法是把几个 CDF 直接套在一起,没有做 Kendall 变换,也没有纤维混合。
既然只是函数代入,它就可能以三种方式翻车。
1. 可能出现负密度
把 $D_k$ 的 CDF 直接喂给外层 $C$,再对 $v$ 求混合偏导,得到的"密度"在 $m_k\ge 2$ 时,某些点上可以是负的。
密度是非负函数,这是概率密度的定义。一旦出现负值,它就已经不是概率密度了。这个反例说明:分布函数的复合,和概率意义上的复合,根本不是一回事。
2. 可能出现积分超过 $1$(或不足 $1$)
如果继续沿用那个错误密度公式
$$g(v)=\phi\big(D(v)\big)\prod_{k=1}^n \psi_k(v^{(k)}),$$
还会出第二个问题:coarea 里漏了 Jacobian,把未归一化的几何纤维测度当成概率条件核。反例 $1+\theta(1-2x_1x_2)$ 算出来
$$\int g = 1+\frac{\theta}{4}\neq 1,$$
概率不守恒。积分大于 $1$、小于 $1$,都是可能的,所以它连概率分布都不是。
3. 破坏边缘均匀性
就算躲过了前两个问题,普通复合也未必是 copula。
当 $D_k$ 是 $m_k\ge 2$ 维 copula 时,$D_k(v^{(k)})$ 是一个联合 CDF 的值。它虽然落在 $[0,1]$,分布却不是均匀的,会偏向 $1$ 附近。外层 $C$ 被喂了非均匀的输入,复合体的对应边缘就不再均匀 → 不是 copula。
就算内部只是一维,只要塞进去的不是均匀变量,而是一个任意的一维 CDF $F(v)$,比如
$$F(v)=v^2,\qquad v\in[0,1],$$
普通复合给出的是 $F(v)$ 自己。它当然是个合法的一维分布函数(单调、从 $0$ 到 $1$),但它的边缘分布是 $F$,不是 Uniform$[0,1]$。
所以它连 copula 的第一条公理——边缘均匀——都过不了。一维也保不住。
小结:三种翻车方式并列,但结论一致:普通复合只是一个普通函数,不是概率对象。
没做 Kendall 变换,输入不均匀,边缘就不均匀 → 不是 copula;
用错密度公式,概率不守恒 → 连概率分布都不是;
直接对分布函数求偏导 → 可能出现负密度。
只有在 Kendall 变换把它拉回均匀之后,它才重新成为 copula。
二、Kendall 复合的定义
记总输入 $v=(v^{(1)},\dots,v^{(n)})$,每块 $v^{(k)}\in[0,1]^{m_k}$,总维度 $m=\sum_k m_k$。
第一步,Kendall 变换:对每块输出套分位数函数 $T_k$,
$$u_k = T_k(D_k(v^{(k)}))\in[0,1].$$
任何随机变量的 CDF 推过去都是均匀的,所以 $u_k\sim\mathrm{Unif}[0,1]$——不管 $D_k$ 是几维。合起来就是总映射
$$G(v)=\big(T_1(D_1(v^{(1)})),\dots,T_n(D_n(v^{(n)}))\big).$$
关键是它让前推底测度变成干净的 Lebesgue:
$$G_\#\Pi = \lambda_n \quad\text{on }[0,1]^n.$$
第二步,纤维混合:固定外层配置 $u$,纤维
$$G^{-1}(u)=\{v:G(v)=u\}$$
是"所有能产出这个 $u$ 的原料"。由 disintegration(Rokhlin),存在唯一(a.e.)Markov 核 $P(u,\cdot)$,集中在纤维上,使得叶测度
$$\Pi:=\bigotimes_{k=1}^n \mu_{D_k}$$
分解为
$$\Pi(A)=\int_{[0,1]^n} P(u,A)\,d\lambda_n(u),\qquad \forall A\in\mathcal{B}([0,1]^m).$$
复合测度定义为
$$\boxed{\mu_\gamma(A):=\int_{[0,1]^n} P(u,A)\,\phi_C(u)\,du}$$
其中 $\phi_C=d\mu_C/d\lambda_n$ 是外层密度(若不假设 a.c.,直接写 $d\mu_C(u)$ 也行)。
注意:这里没有密度复合,只有条件概率复合。
三、证明它仍是 $m$-copula
一个 $m$-copula 要满足四条。
1. 是概率测度(免费)
$P(u,\cdot)$ 对每个 $u$ 是概率核,$\phi_C\ge 0$,$\int\phi_C=1$,所以
$$\mu_\gamma([0,1]^m)=\int P(u,[0,1]^m)\,\phi_C(u)\,du =\int 1\cdot\phi_C(u)\,du=1.$$
✅
2. Grounded(两头堵死)
取 $A=\{v:v_i^{(k)}=0\}$。
若 $v_i^{(k)}=0$,因 $D_k$ grounded,$D_k(v^{(k)})=0$,又 $T_k(0)=0$,故 $u_k=0$。
对 $u_k>0$:纤维 $G^{-1}(u)$ 要求第 $k$ 块产出正数,不可能含"坐标等于 $0$"的点,故 $P(u,A)=0$。
对 $u_k=0$:$D_k$ 在那条低维边上测度为 $0$,故 $\Pi(A)=0$;由 disintegration 得 $P(u,A)=0$ 对 a.e. $u$ 成立。又外层 $C$ grounded,$\mu_C(\{u_k=0\})=0$。
两层都归零,所以
$$\mu_\gamma(A)=0.$$
对所有块、所有坐标成立 → grounded。✅
3. 边缘均匀(每一刀都是 $w$)
固定块 $k$、块内坐标 $i$,令 $A_w=\{v_i^{(k)}\le w\}$。
叶测度上,第 $k$ 块的边缘是 $\mu_{D_k}$,而 $D_k$ 是 copula,故
$$\Pi(A_w)=\mu_{D_k}(\{v_i^{(k)}\le w\})=w.$$
disintegration 给出
$$\int_{[0,1]^n} P(u,A_w)\,d\lambda_n(u)=w.$$
引理:对 $\lambda_n$-a.e. $u$,
$$P(u,A_w)=w.$$
理由:$D_k$ 的第 $i$ 边缘本来就是 Lebesgue;Kendall 变换只把 $u_k$ 重新均匀化,不改变条件于 $u_k$ 时块内坐标的均匀边缘。所以给定 $u_k$,块内第 $i$ 坐标的条件分布仍是 $\mathrm{Unif}[0,1]$,累积到 $w$ 就是 $w$。
于是
$$\mu_\gamma(A_w)=\int w\,\phi_C(u)\,du = w\int\phi_C(u)\,du = w.$$
每个坐标边缘都是均匀 → 满足 copula 边缘公理。✅
4. $m$-increasing(自动)
$\mu_\gamma$ 是非负概率测度,对超矩形的 $m$ 维体积自动非负。✅
四、为什么全程不需要绝对连续
用到的只有:
copula 的 grounded 与边缘均匀(公理);
$\Pi=\bigotimes \mu_{D_k}$ 是乘积概率测度;
disintegration 定理;
$G_\#\Pi=\lambda_n$(Kendall 变换给的)。
没有一处要求密度存在。绝对连续性是后面写显式密度、算熵、算 KL 时才用得到的东西;它让"可算",但不参与"拼出来还是不是 copula"。
Kendall 变换的真正作用就一句:把底测度钉成 Lebesgue,让纤维混合的基准对齐,于是相容性自动满足。
五、所以为什么叫 Kendall 复合
普通 Sklar 复合:只是函数代入,没有概率保证 → 可能出负密度、可能积分不为 $1$、可能边缘不均匀。
Kendall 复合 = Kendall 变换(拉回均匀)+ 纤维混合(测度层拼接)。
它拼出来的东西:
是概率测度;
grounded(两头堵死);
边缘均匀(每一刀厚度都是 $w$);
$m$-increasing(自动)。
而且对所有 copula 成立,不挑绝对连续,可以一代一代迭代下去。
闭环是这样的:
$$\boxed{ \begin{aligned} &\text{边缘均匀}\;\to\;\text{Lebesgue 参考}\\ \to\;&\text{Kendall 把块输出拉回均匀}\;\to\;\text{底测度对齐}\\ \to\;&\text{disintegration 条件核自动归一化}\;\to\;\text{Kendall 复合仍是 copula}\\ \to\;&\text{边缘仍均匀}\;\to\;\text{可继续复合}. \end{aligned}}$$一句话总结
普通复合只是函数代入,不是概率对象;它可能出现负密度,也可能积分超过 $1$,还可能破坏边缘均匀。Kendall 复合之所以正确,是因为 Kendall 变换强制每块输出回到均匀,让外层永远吃到合法输入;而拼出来仍是 copula,是因为叶测度的均匀边缘让每一刀厚度都是 $w$、钉死坐标时测度归零——两件事都不依赖绝对连续性,是公理给的。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2026-10-22 13:12
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社