||
Copula Magma 为什么是可能的?
这篇文章只讲一件事:为什么在 copula 的世界里,原本拼不起来的东西,终于能拼了。我们不讲证明,也不举例子,只把设计原因说清楚。
先说 DisCo 的问题
在上一篇 DisCo 代数里,我们想把高维联合分布拆成块,再拼回去。拆容易,拼难。
拼的时候要回答一个问题:给定几个标量摘要,每一块在里面具体长什么样?这个答案就是条件概率。DisCo 的拼法,就是对骨架上的每一点,取出各块的条件分布,独立地乘起来,再按骨架加权求和。
这个拼法有两个毛病。
第一,能不能拼要看运气。 骨架的边缘(它每一维的分布)和各块摘要的分布,本该是同一个东西。但在 DisCo 里,这两者是分别造出来的,通常一致,但不是自动一致。只有满足相容性条件时,拼接在代数上才说得通;不满足就拼不了。这就是 partial:运算只对一部分输入有定义。
第二,拼出来可能不唯一。 条件概率只是几乎处处唯一的。而 DisCo 里的积分底测度是骨架,骨架是耦合的,不是各边缘乘起来的独立乘积。于是可能出现:某个集合在块自己的底测度下是零测的,但在骨架这个耦合测度下不是零测。这时在零测集上换个条件版本,就会实实在在地改变拼出来的结果。
一句话总结 DisCo 的困境:没有统一参考测度,导致相容性要手验,运算 partial;底测度对不上,导致条件版本有歧义,结果可能不唯一。
换到 copula:尺子自己长出来了
现在把舞台从一般概率分布换成 copula。copula 就是边缘全都是均匀分布的联合分布,它只描述变量之间的依赖结构。
这一句话,就是整篇文章的魔法来源。
copula 生活在单位立方体上,每个边缘都是 Unif[0,1]。而 [0,1] 上标准的均匀测度就是 Lebesgue 测度。于是第一块的参考测度是它,第二块也是它,第 k 块还是它,乘起来仍然是 Lebesgue 测度。
不需要选,也不需要假设。copula 的边缘均匀性,直接把参考测度钉死为 Lebesgue。DisCo 里我们苦苦寻找的统一尺子,在 copula 世界里是公理免费送的。
Kendall 变换:让输出也均匀
有了统一的尺子,只解决了一半问题。
copula 的边缘是均匀的,这很好。但我们拼接时用的输入,是块的 CDF 摘要,一个标量。这个标量的分布不一定是均匀的。它的取值虽然在 [0,1] 里,但分布可能是偏的、堆在一头的。
于是出现一个尴尬局面:块的底测度是均匀的,但输出不是均匀的。输入输出不在一个尺度上,接口对不齐。
Kendall 变换干的事,就是把块的 CDF 再做一次均匀化。它的作用是把一块变量的累积概率重新校准成均匀分布。为什么一定行?因为一个随机变量的累积分布函数,把它自己推过去,得到的一定是均匀分布。
所以输入是均匀的,输出经过 Kendall 变换也变成均匀的。输入输出都落在同一个 Lebesgue 尺子下,接口对齐了。
为什么相容性自动成立
现在看清楚,为什么那个让人头疼的相容性条件不再是问题。
回忆相容性要验的是:骨架的第 j 个边缘,等于第 j 块的摘要测度。
在 copula 加 Kendall 的世界里,左边是骨架的边缘。骨架是把各块的 Kendall 标量放在一起得到的,每个 Kendall 标量的分布都是均匀的,所以骨架的第 j 个边缘就是 Lebesgue 测度。
右边是块的摘要测度。第 j 块本身是一个 copula,copula 的边缘就是均匀的,所以块摘要测度也是 Lebesgue 测度。
于是两边都是 Lebesgue 测度。这不是验证通过,这是恒等式。
在 DisCo 里要手写手验的相容性条件,到了 copula 世界,变成了 λ 等于 λ。
它为什么重要?还记得歧义是怎么来的吗?条件版本在某个零测集上换了,积分就变了,因为积分底测度和块的底测度对不上。现在块的底测度是 λ,骨架的边缘也是 λ,两者对齐了。所以在 λ 零测集上换条件版本,不会改变骨架积分的结果,歧义消失。
从 partial 到 total
把这几步串起来。
copula 边缘均匀,给了统一参考测度。Kendall 变换,让输出也均匀,输入输出接口对齐。两边都是 λ,相容性自动成立,不再需要看运气才能拼。底测度对齐,条件版本无歧义,拼出来唯一。
于是拼接运算从 partial 变成 total:在 admissible 输入上都能拼,且没有歧义。
这就是 copula magma 为什么可能的核心答案:不是我们想出了更聪明的拼法,而是 copula 的公理加上 Kendall 变换,把拼法成立的全部前提自动满足了。
但为什么不是 monoid,也不是operad
拼法能处处运算,不代表它"好算"。
考虑三个块 A、B、C,有两种拼法。左结合是先拼 A、B,再和 C 拼。右结合是先拼 B、C,再和 A 拼。这两种拼法结果一般不一样。
直觉很好理解。先拼 A、B,意味着在 A、B 之间做一次压缩成标量,这一步会丢失 A、B 之间的一部分联合依赖。先拼 B、C,丢失的是 B、C 之间的依赖。丢失的东西不同,最后拼出来的联合分布就不同。
所以这个运算可以处处拼,而且有单位元:独立 copula $\Pi(u)=\prod u_i$ 就是它。任何 copula 和 $\Pi$ 拼,都等于自己。但它不满足结合律,而且它的乘法是 operadic 的,不是简单的二元乘积:一次是让一个外层 copula 同时拼 k 个块,而不是两个东西相乘。
这两个原因决定了它不是 monoid。monoid 要求二元运算、满足结合律、有单位元,三个条件缺一不可;而这里一次就拼多个块,天生不是二元乘积,所以就算有了结合律,也不归 monoid 管。
它不是 operad,原因只有一个:没有结合律。准确地说,它是一个 operadic magma,也就是无结合律的多元复合。
闭环:为什么可以一代一代往下拼
还要问一个问题:拼出来的东西,还配得上叫 copula 吗?
如果拼出来的是一个一般分布,边缘不再均匀,那下一代拼接时,Lebesgue 这把尺子就失效了,整个 tower 会一层层塌掉。
幸运的是,复合之后仍然是 copula。每一块的边缘是均匀的,经过 Kendall 校准后输出也均匀,拼接后每一坐标的边缘仍然均匀。
于是得到一条闭环:边缘均匀,得到 Lebesgue 参考;相容性自动,得到 total;复合仍是 copula,边缘仍均匀,于是可以迭代。正因为有这个闭环,我们才能一层层拼下去,把 copula magma 做成一个稳定的代数结构,而不是拼一次就散架。
一句话总结
copula magma 之所以可能,是因为三件事同时成立。
第一,copula 的边缘均匀性,免费送了一把统一的尺子,就是 Lebesgue 测度。第二,Kendall 变换让输出也均匀,把输入输出接口对齐。第三,两边都是 λ,让原本要手验的相容性自动成立,拼法不再有歧义。
于是 partial 变成 total;运算有单位元但不结合,乘法是多元的,所以是 operadic magma。
这不是技巧上的改进,而是舞台换了:在一般概率测度上要手写的条件,在 copula 世界里,是公理替你写好的。
余论:copula 可能是一个范畴性的东西
把这条链再往上看一层,会发现它暗示了更多的东西。
copula 的边缘均匀性,把统一参考测度钉死为 Lebesgue;Lebesgue 加上 Kendall 变换,让相容性自动成立,拼接从 partial 变成 total;复合之后仍然是 copula,于是整个结构可以一代一代往下拼。
注意一个事实:上面这一步步,对所有的 copula 都成立,我们没有假设绝对连续。 均匀边缘是 copula 的公理,Kendall 变换对任何 copula 的 CDF 都能做,相容性是 $\lambda=\lambda$ 的恒等式。这个结构不是靠加假设换来的,是 copula 本身的公理给的。
这三件事合在一起,给出的是一张组合的图:
对象:copula;
态射:从一个外层 copula 到若干子块 copula 的复合;
复合律:多元拼接,不结合;
单位:独立 copula $\Pi$。
它有了"对象"和"态射"两层,而"态射"本身又是一个 copula。所以一个自然的追问是:
copula 很可能是一个范畴性的东西。
这不是已经证明的命题,而是一个纲领、一个公开问题:
怎样把 copula 理论范畴化?
沿着这个问题往下走,至少有几个方向:
范畴的公理够不够? 复合不结合,说明它不是普通的范畴。一个自然的候选是非结合范畴,或者说以 operadic magma 为 Hom 结构的范畴。
Kendall 变换在范畴里是什么? 它是把每个 copula 的 CDF 送到均匀分布的映射。这条链的"自动相容",也许可以写成一个交换图。
从 copula 到 Markov 核:copula 描述边缘均匀的依赖结构,Markov 核描述条件转移。两者之间是否有一个伴随?如果有,copula 范畴和概率核范畴就是同一件事的两种表述。
更高层的范畴化:如果一层 copula 已经能拼,那么"copula 的 copula"——把 copula 之间的复合关系再看作对象——是否形成一个 2-范畴?这会直接指向依赖结构的高阶同调。
这不是这篇文章要回答的问题。这篇文章只做了一件事:把 copula magma 为什么可能,讲清楚。
但这个问题值得留在这里。
因为一旦 copula 被范畴化,依赖结构的拼接、条件化、信息流,就都有了图论和同调的语言;copula 熵的"曲率"项,也许可以写成一个边缘算子的作用。
这是一个公开问题,也是一个纲领。
如果"边缘均匀 → Lebesgue → total 复合"这条链不是孤立的技巧,而是某个范畴的结构定理,那么 copula 理论的下一个篇章,就不在概率论里,而在范畴论里。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2026-10-22 13:12
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社