Copula Operad 的正确打开方式：纤维上的概率混合

一句话摘要：我们之前那篇论文里的 copula operad 密度公式是错的。但 operad 本身没死——只是得退到测度层、限制到绝对连续类、用 Kendall 变换把纤维"拉直"，然后在每一根纤维上做概率混合。本文把这个正确版本完整写出来，重点是结合律的严格证明和"纤维混合"的直觉。

0. 发生了什么

先交代背景，不然读者不知道我们在修什么。

论文《Copula Operad and Copula Entropy》里，作者给出了一个复合公式：

$$g(v) = \phi\big(F(v)\big) \cdot \prod_{k=1}^n \psi_k(v^{(k)})$$

声称这是复合后 copula 的密度。然后我（和审稿人）发现：这个公式在 $m_k \ge 2$ 时积分不等于 1。反例：取二维 copula 密度 $1+\theta(1-2x_1x_2)$，按公式算出来 $\int g = 1+\theta/4 \neq 1$。

根源是什么？coarea 公式里漏了 Jacobian 因子。更根本地说：论文把 Federer 几何纤维测度（未归一化）当成了概率条件核（已归一化）。两者之间差着一个 $\prod \sigma_k$ 的"汇率"。

但这不意味着 operad 本身是错的。测度层的定义是站得住的——只是需要小心地重做，并且限制到一个安全的子类。

本文就是那个重做后的干净版本。

1. 核心直觉：纤维上的混合

先不讲定义，讲画面。

1.1 一张图

想象你有 $n$ 个"原料" copula $D_1, \dots, D_n$，分别吃 $m_1, \dots, m_n$ 个变量，吐出 $n$ 个中间值 $(u_1, \dots, u_n)$。

外层 copula $C$ 说："我不关心你们内部怎么来的，我只管这 $n$ 个输出 $(u_1,\dots,u_n)$ 之间的依赖结构。"

复合操作在做什么？

固定外层 $C$ 给定的一个"配置" $u = (u_1,\dots,u_n)$，然后在纤维 $F^{-1}(u)$ 上，按 $C$ 的条件分布挑一个具体的 $v$。

换句话说：

• 纤维​ $F^{-1}(u) = \{v : D_k(v^{(k)}) = u_k,\ \forall k\}$ 是"所有能产生这个 $u$ 的原始输入"

• 混合​ = 在每根纤维上，按外层 copula 的条件概率挑点

• 整体​ = 对所有 $u$ 加权平均，权重由 $C$ 决定

这就是"纤维混合"的画面。不需要密度，不需要 Jacobian，只需要条件概率。

1.2 为什么 Kendall 变换让一切变干净

问题：纤维 $F^{-1}(u)$ 的形状很怪（尤其 $m_k \ge 2$ 时是低维流形），条件分布算起来要 coarea。

解法：在每个输出坐标 $u_k$ 上再套一个 Kendall 变换​ $T_k$——就是分位数函数的逆。它的效果是：

把 $D_k$ 输出的边缘分布（可能歪歪扭扭）拉成标准的 $[0,1]$ 上的均匀分布。

拉平之后：

• 前推测度 $G_\#\Pi$ = Lebesgue on $[0,1]^n$（干干净净的均匀）

• 纤维的几何形状虽然没变，但底测度现在是 Lebesgue，disintegration 的"基准"对齐了

• 条件核自动归一化（因为底 = 前推）

Kendall 变换不消灭复杂性，但它把复杂性从"底测度"挪到了"纤维内部"——而纤维内部的事由概率 disintegration 定理全权接管。

2. 正式定义设定

• $\mathbf{Cop}^{\mathrm{ac}}_n$​ = 所有 绝对连续​ $n$-copula 的集合（即测度 $\mu \ll \lambda_n$，有密度）

• 对 $C \in \mathbf{Cop}^{\mathrm{ac}}_n$，记 $\mu_C$ 为其测度，$\phi_C = d\mu_C/d\lambda_n$ 为密度

复合映射

给定：

• 外层 $C \in \mathbf{Cop}^{\mathrm{ac}}_n$

• 内层 $D_k \in \mathbf{Cop}^{\mathrm{ac}}_{m_k}$，$k=1,\dots,n$

定义总映射（含 Kendall 标准化）：

$$G : [0,1]^m \to [0,1]^n, \quad G(v) = \big(T_1(D_1(v^{(1)})),\; \dots,\; T_n(D_n(v^{(n)}))\big)$$

其中 $T_k = F_{U_k}^{-1}$ 是第 $k$ 个边缘分布的分位数函数（Kendall 变换），它保证 $G_\#\Pi = \lambda_n$（Lebesgue）。

叶测度$$\Pi := \bigotimes_{k=1}^n \mu_{D_k} \quad \text{（在 } [0,1]^m \text{ 上的概率测度）}$$定义 2.1（Operad 复合，测度版）

由于 $G_\#\Pi = \lambda_n$，由概率 disintegration 定理（Rokhlin），存在唯一的（up to $\lambda_n$-a.e.）Markov 核 $P : [0,1]^n \times \mathcal{B}([0,1]^m) \to [0,1]$ 使得：

$$\Pi(A) = \int_{[0,1]^n} P(u, A)\, d\lambda_n(u) \qquad \forall A \in \mathcal{B}([0,1]^m)$$

对每个 $u$，$P(u, \cdot)$ 集中在纤维 $G^{-1}(u)$ 上，且是概率测度。

复合测度定义为：

$$\boxed{\mu_\gamma(A) := \int_{[0,1]^n} P(u, A)\, \phi_C(u)\, du}$$

即：先在纤维上按 $P(u,\cdot)$ 混合，再按外层密度 $\phi_C$ 加权。

注：因为 $C$ 是 a.c. 的，$\mu_C(A) = \int_A \phi_C(u)\,du$，所以这等价于 $\mu_\gamma(A) = \int P(u,A)\,d\mu_C(u)$。写成 $\phi_C(u)\,du$ 只是显式强调底测度是 Lebesgue。

3. 为什么这是合法的 $m$-copula定理

 3.1（Well-definedness）

$\mu_\gamma$ 是 $[0,1]^m$ 上的概率测度，且满足 copula 的四条公理。

证明：

(a) 概率测度

• $P(u, \cdot)$ 是概率核：$\forall u,\ P(u, [0,1]^m) = 1$

• $\phi_C \ge 0$，$\int \phi_C(u)\,du = 1$

• $\mu_\gamma([0,1]^m) = \int_0^1 P(u,[0,1]^m)\phi_C(u)\,du = \int \phi_C(u)\,du = 1$ ✅

(b) Grounded

设 $A = \{v : v_i^{(k)} = 0\}$。若 $v_i^{(k)} = 0$，则 $D_k(v^{(k)}) = 0$，故 $u_k = 0$。

对 $u_k > 0$：$G^{-1}(u)$ 中无满足 $v_i^{(k)}=0$ 的点，故 $P(u, A) = 0$。

对 $u_k = 0$：由 $D_k$ 的 groundedness，$\mu_{D_k}(\{v_i^{(k)}=0\}) = 0$，故 $\Pi(A) = 0$。由 disintegration 方程 $\int P(u,A)\,du = 0$ 且 $P \ge 0$，得 $P(u,A) = 0$ 对 a.e. $u$（含 $u_k=0$ 的）。

又 $\phi_C$ 只在 $u_k=0$ 时可能非零，但 $\int_{\{u_k=0\}} \phi_C(u)\,du = \mu_C(\{u_k=0\}) = 0$（$C$ grounded）。

故 $\mu_\gamma(A) = 0$。✅

(c) 均匀边际

固定块 $k$、坐标 $i$，令 $A_w = \{v_i^{(k)} \le w\}$。

由 disintegration：$\int P(u, A_w)\,du = \Pi(A_w) = w$（因为 $\Pi$ 的第 $k$ 块边际是 $\mu_{D_k}$，而 $D_k$ 是 copula，$\mu_{D_k}(\{v_i^{(k)}\le w\}) = w$）。

引理：对 $\lambda_n$-a.e. $u$，$P(u, A_w) = w$。

理由：给定 $D_k(v^{(k)}) = u_k$ 时，$v_i^{(k)}$ 的条件分布是 uniform on $[0,1]$——因为 $D_k$ 是 copula，其第 $i$ 个边际就是 Lebesgue，而 Kendall 变换 $T_k$ 把 $u_k$ 拉成均匀，故条件边际不变。所以 $P(u, \{v_i^{(k)}\le w\}) = w$ 对 a.e. $u$。

于是：

$$\mu_\gamma(A_w) = \int_0^1 P(u, A_w)\phi_C(u)\,du = \int_0^1 w\cdot\phi_C(u)\,du = w$$

✅

(d) $m$-increasing：非负测度自动满足。✅

∎

4. 结合律（重点）

这是 operad 公理里最需要care的部分。

设定

三棵树：

• $\beta$：arity $(p, q, r)$ → 输出 $N$ 维（父节点）

• $\gamma$：arity $(m_1, m_2)$ → 输出 $p$ 维（左子树）

• $\delta$：arity $(n_1, n_2)$ → 输出 $q$ 维（右子树）

（为简洁，假设二叉树；多叉同理。）

总叶子数：$M = m_1+m_2+n_1+n_2$。

总映射：

• 左子树映射：$G_\gamma : [0,1]^{m_1+m_2} \to [0,1]^p$

• 右子树映射：$G_\delta : [0,1]^{n_1+n_2} \to [0,1]^q$

• 父节点映射：$G_\beta : [0,1]^{p+q} \to [0,1]^N$

两级复合的总映射（左结合）：

$$\tilde{G}_{\text{left}}(v) = G_\beta\big(G_\gamma(v_L),\, G_\delta(v_R)\big)$$

右结合的总映射：

$$\tilde{G}_{\text{right}}(v) = G_\beta\big(G_\gamma(v_L),\, G_\delta(v_R)\big)$$

它们是一样的函数！​ 因为函数复合满足结合律，$G_\beta \circ (G_\gamma \times G_\delta)$ 不管括号怎么放，结果都是同一个映射 $[0,1]^M \to [0,1]^N$。

定理 4.1（结合律）$$\mu_{\beta \circ (\gamma,\delta)} = \mu_{(\beta \circ \gamma) \circ \delta}$$

证明：

记总叶测度 $\Pi_{\text{total}} = \mu_{D_\gamma^{(1)}}\otimes\mu_{D_\gamma^{(2)}}\otimes\mu_{D_\delta^{(1)}}\otimes\mu_{D_\delta^{(2)}}$。

记总映射 $\tilde{G} = G_\beta \circ (G_\gamma \times G_\delta)$（两种括号方式给出同一个函数）。

由 Definition 2.1，两种括号方式给出的复合测度分别是：

$$\mu_{\text{left}}(A) = \int_{[0,1]^N} P_{\tilde{G}}(w, A)\, \phi_\beta(w)\, dw$$$$\mu_{\text{right}}(A) = \int_{[0,1]^N} P_{\tilde{G}}(w, A)\, \phi_\beta(w)\, dw$$

完全相同——因为 disintegration 核 $P_{\tilde{G}}(w, \cdot)$ 只依赖 $\Pi_{\text{total}}$ 和 $\tilde{G}$ 这两个对象，而两者在两种括号方式下不变。

更严格地说：由 disintegration 的唯一性（Rokhlin），给定 $\Pi_{\text{total}}$ 和 $\tilde{G}$，条件核 $P_{\tilde{G}}(w, \cdot) = \Pi_{\text{total}}(\cdot \mid \tilde{G} = w)$ 在 $\lambda_N$-a.e. $w$ 上唯一确定。与"你先拆左子树还是先拆右子树"无关。

而外层密度 $\phi_\beta$ 在两种情形下也是同一个函数（同一个父 copula $\beta$）。

故积分相等。∎

直观解释

结合律成立是因为：纤维混合的顺序不影响最终结果。​ 不管你先"在左子树的纤维上混合"还是"先在整个空间的纤维上混合"，最终你都是在问同一个问题："给定最终输出 $w$，原始叶子变量 $v$ 的条件分布是什么？" 而条件分布只认总映射和总测度，不认括号。

这就像：不管你先切西瓜再切橙子，还是先切橙子再切西瓜，最终"给定成品是这模样，原料从哪来"这个问题的答案是唯一的。

5. 单位元与等变性单位元

$\text{id} \in \mathbf{Cop}^{\mathrm{ac}}_1$ 就是 Lebesgue 测度 on $[0,1]$（密度 = 1）。

$G_{\text{id}}(v) = T_1(v^{(1)}) = v^{(1)}$（因为 1-维 copula = 恒等，Kendall 变换也是恒等）。

纤维 $G^{-1}(u) = \{v : v^{(1)} = u\} \cong [0,1]^{m-1}$，条件核 $P(u, \cdot) = \lambda_{m-1}$（剩余坐标的 Lebesgue）。

$$\mu_{\text{id}\circ D}(A) = \int_0^1 \lambda_{m-1}(A \cap \{v^{(1)}=u\})\cdot 1\,du = \mu_D(A)$$

✅

等变性

坐标置换 $\sigma$ 与 $G$ 可交换（因为 $G$ 按块逐坐标作用），故条件核随之置换，复合测度也置换。平凡。​ ✅

6. 总结：这个 operad 长什么样

和原论文的关系

• ✅ Operad 骨架：保留，但定义在测度层，不碰密度

• ❌ 密度公式 (3)：删除（仅 $m_k=1$ 时退化为平凡）

• ❌ Proposition 2.1（a.c. 封闭性原证明）：撤回，但结论在修正定义下成立（只是证明不同）

• ❌ 熵公式（Theorem 3.1）：需另起炉灶（密度不存在时熵无定义；若限制到有密度的子类，需重推）

7. 一句话收尾

Copula operad 不是错的，只是之前长在了错误的层面上。密度是表象，测度是骨架，纤维混合是灵魂。把骨架擦干净、限制到 a.c. 类、用 Kendall 把纤维拉直——operad 自己就站起来了，结合律甚至不需要算，它本来就在那儿。

写于和 @Xuexing Lu 的深夜debug sessions之后。致谢：那个 $\theta/4$ 的反例救了这篇短文。

支持单位：