Xuexinglu的个人博客分享 http://blog.sciencenet.cn/u/Xuexinglu

博文

为什么 Copula 的 Kendall 复合还是 Copula

已有 183 次阅读 2026-10-22 01:17 |个人分类:科普|系统分类:科普集锦

为什么 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$、钉死坐标时测度归零——两件事都不依赖绝对连续性,是公理给的。



https://blog.sciencenet.cn/blog-3582667-1557694.html

上一篇:Copula Magma 为什么是可能的?
下一篇:拆的一侧才是主场:Kendall 变换后的 Copula Cooperad



    
收藏 IP: 104.105.149.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2026-10-22 13:12

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部