目录

Agreement and Statistical Efficiency in Bayesian Perception Models


作者Yash Deshpande, Elchanan Mossel, Youngtak Sohn
机构MIT (Department of Mathematics / IDSS)
版本arXiv:2205.11561v3 [math.ST], 2023-08-09
链接https://arxiv.org/abs/2205.11561

一群人在图上互相发言,每个人每轮只从自己的后验里抽一个样本发出去(甚至只发一个 0/1 的硬币结果),而不是像经典经济学模型那样汇报后验均值。本文证明:只要图连通、而且每个人知道对方是在抽样,那么「听到的样本频率 ≈ 对方的整个后验」,于是等价于直接交换后验——所有人最终一致,且一致到的那个后验恰好等于把所有私有信号汇总后的最优贝叶斯后验。反过来,如果听者误以为对方发来的是世界给的新独立观测,整个系统数学上退化成两个互相扔球的 Pólya urn:大家仍然达成一致,但收敛到一个随机的点质量(方差→0 = 过度自信),且几乎必然不等于正确答案。

✓一致性:任意连通图
✓统计有效:= 全信息最优后验
1 bit每轮通信量足够
✗误规约:一致但过度自信

研究动机

术语速查

这篇论文的抽象符号都有非常具体的对应物,先把它们落地,后面才看得懂。

符号 / 术语 在本文中具体是什么 LLM 类比
\(\theta\) 未知世界状态。Bernoulli 模型里就是一枚硬币的偏差 \(\theta\sim\mathrm{Unif}([0,1])\);高斯模型里是 \(\theta\sim\mathsf N(0,1)\) 的一个标量。 被问的那个问题的真实答案
\(S_i\) agent \(i\) 的初始私有信号,给定 \(\theta\) 条件独立。Bernoulli:抛一次 \(\mathrm{Ber}(\theta)\);高斯:\(\mathsf N(\langle a_i,\theta\rangle,\Sigma_i)\),\(a_i\) 是「信号强度」。 这个 LLM 自己的预训练数据
\(G=(V,E)\) 社交网络图,\(|V|=n\) 有限。唯一要求是连通(不连通就对每个连通分量分别套结论)。最小的非平凡例子就是一条边(2 个 agent)。 谁在和谁对话
\(\nu_i(t)\) agent \(i\) 在第 \(t\) 轮的后验分布 \(\mathrm{Law}(\theta\mid H_i(t))\)。注意:agent 的后验计算是精确的,不做近似推断。 模型当前对答案的完整不确定性
\(Y_i(t)\) 发言:从 \(\nu_i(t)\) 抽的一个样本(给定历史后与他人条件独立)。sampling model 中直接把这个实数发给邻居。 这一轮采样出来的那条回答
\(X_i(t)\) 二元表态:Bernoulli 模型里再多抛一次硬币 \(X_i(t)\sim\mathrm{Ber}(Y_i(t))\),只把这 1 个 bit 发给邻居。 只输出「是/否」的极端压缩
\(H_i(t), K_i(t)\) agent \(i\) 到 \(t\) 时刻看到的全部信息(\(\sigma\)-代数):自己的 \(S_i\) 加上所有邻居历史上发来的消息。\(H_i(\infty)=\cup_t H_i(t)\)。 预训练数据 + 完整对话历史
最优贝叶斯后验 \(\mathrm{Law}(\theta\mid S_1,\dots,S_n)\):假想有个上帝把所有人的私有信号汇总后算出的后验。这是理论上限。 把所有模型的训练数据合起来重训一遍
probability matching 按后验概率比例随机行动,而非取后验均值/最大值。心理学实验反复观察到的人类行为,也是本文 sampling model 的建模出发点。 temperature>0 的采样解码
misspecification(误规约) agent 1 认为收到的 \(\tilde X_2(s)\) 是 i.i.d. \(\mathrm{Ber}(\theta)\)(世界给的新观测),而真实是 \(\mathrm{Ber}(\tilde Y_2(s))\)(对方后验的采样)。模型对消息怎么产生的这件事搞错了。 把别的模型的输出当成新爬到的网页

核心方法

1. sampling model:只传一个采样值的社会学习

设定:未知状态 \(\theta\in\mathbb R^d\) 服从公共先验 \(\nu\);\(t=0\) 时每个 agent \(i\) 拿到私有信号 \(S_i\sim P_i(\cdot\mid\theta)\),给定 \(\theta\) 条件独立(这一族 \(\{P_i(\cdot\mid\theta)\}\) 叫 signal structure)。此后每一轮 \(t\ge1\):

  1. 精确算后验:\(\nu_i(t):=\mathrm{Law}(\theta\mid H_i(t))\),其中 \(H_i(t)=\sigma\big(S_i,(Y_j(s))_{j\in\partial i,\,s<t}\big)\)(自己的信号 + 邻居过去发的全部消息)。
  2. 抽一个样本:\(Y_i(t)\sim\nu_i(t)\),给定 \(H_i(t)\) 后与其他一切条件独立。
  3. 发给邻居:把 \(Y_i(t)\) 广播给所有 \(j\in\partial i\)。所有 agent 同步进行。

和 belief model 的区别只有一处:发出去的是后验的一个样本,不是后验均值。但这一处区别让整个过程变成一个带噪声的随机系统,收敛性完全不显然——作者强调,哪怕 \(G\) 只是一条边,收敛性质也是高度非平凡的。

flowchart LR
  S1["private signal S1"] --> P1["posterior nu_1(t)"]
  P1 -->|"draw one sample"| Y1["Y1(t)"]
  Y1 -->|"broadcast"| P2
  S2["private signal S2"] --> P2["posterior nu_2(t)"]
  P2 -->|"draw one sample"| Y2["Y2(t)"]
  Y2 -->|"broadcast"| P1
  Y1 -.->|"extra coin flip"| X1["X1(t) in 0,1"]
  Y2 -.->|"extra coin flip"| X2["X2(t) in 0,1"]
          
drag to pan · scroll to zoom
sampling model 的一轮:各自精确算后验 → 抽一个样本 → 广播给邻居。虚线是 Bernoulli 变体:对采样值再抛一次硬币,只发 1 个 bit。

2. 定理 2.1:一致性,靠「模仿原理」

Theorem 2.1:对每个 agent \(i\),随机后验 \(\nu_i(t)\) 几乎必然弱收敛到 \(\nu_i(\infty)=\mathrm{Law}(\theta\mid H_i(\infty))\);而且若 \(i,j\) 相邻则 \(\nu_i(\infty)=\nu_j(\infty)\) a.s.。因此 \(G\) 连通时所有人渐近达成一致。

收敛部分是鞅收敛定理(Lévy upward theorem 作用在分布函数 \(F_{it}(u)=\mathbb P(\theta\le u\mid H_i(t))\) 上,先在可数稠密集 \(\mathbb Q^d\) 上收敛,再推到全部连续点)。一致性部分才是难点,用的是 sampling 版本的 imitation principle(模仿原理):

直观理解:因为 \(Y_j(t)\) 在给定 \(H_j(t)\) 后是条件独立采样的,所以经验频率 \(\frac1T\sum_{t\le T}\mathbf 1(Y_j(t)\le u)\) 会收敛到 \(F_j(u)\)——agent \(i\) 只要听得够久,就能把邻居 \(j\) 的整个极限后验估出来。于是如果 \(i\) 和 \(j\) 的极限信念不同,\(i\) 就可以构造一个凸组合估计量 \(\lambda F_i(u)+\frac{1-\lambda}{T}\sum_t\mathbf 1(Y_j(t)\le u)\),它是 \(H_i(\infty)\)-可测且无偏的,而且(由引理 3.1)严格降低估计 \(\mathbf 1(\theta\le u)\) 的均方误差。但 \(F_i(u)\) 本身就是条件期望,是 \(H_i(\infty)\) 下最小均方误差的估计——矛盾。

证明骨架(引理 3.1 + 三项 Minkowski 分解)

Lemma 3.1:若 \(Z_1,Z_2\) 二阶矩有限,且 \(\min_{\lambda\in[0,1]}\mathbb E[Z(\lambda)^2]=\mathbb E[Z_1^2]=\mathbb E[Z_2^2]\)(其中 \(Z(\lambda)=\lambda Z_1+(1-\lambda)Z_2\)),则 \(Z_1=Z_2\) a.s.。证法:把 \(\mathbb E[Z(\lambda)^2]\) 展成 \(\lambda\) 的二次函数,它在 \(\lambda=0\) 和 \(\lambda=1\) 同时取最小 \(\Leftrightarrow\) 它是常函数 \(\Leftrightarrow M_{11}=M_{12}=M_{22}\Rightarrow\mathbb E[(Z_1-Z_2)^2]=0\)。

取 \(Z_1=F_i(u)-\mathbf 1(\theta\le u)\)、\(Z_2=F_j(u)-\mathbf 1(\theta\le u)\),反设存在 \(\lambda,\varepsilon\) 使凸组合的 MSE 比 \(F_i\) 小 \(\varepsilon\)。对估计量 \(Z_1(T)=\lambda F_i(u)+\frac{1-\lambda}{T}\sum_{t=1}^T\mathbf 1(Y_j(t)\le u)\) 做 Minkowski 分解成三项:

  • \(\Xi_1\le(\mathbb E[(F_i(u)-\mathbf 1(\theta\le u))^2]-\varepsilon)^{1/2}\)(反设)
  • \(\Xi_2\le T^{-1/2}\):因为 \(\{\mathbf 1(Y_j(t)\le u)-F_{jt}(u)\}_t\) 在不同时刻互不相关(采样机制 + tower property),采样噪声以 \(1/\sqrt T\) 被平均掉
  • \(\Xi_3\le T_0(\delta)/T+\delta^{1/2}\):因为 \(F_{jt}(u)\to F_j(u)\)(鞅收敛 + 控制收敛)

先取 \(\delta\ll\varepsilon\) 再取 \(T\) 充分大,便与「条件期望是 MSE 最优」矛盾。

3. 定理 2.4:一致到的那个后验就是最优后验

一致不等于学得对(action model 就是一致但不最优)。要证统计有效性需要对 signal structure 加两个温和条件:

Theorem 2.4:在上述条件下且 \(G\) 连通,极限共同后验 \(\nu(\infty)=\mathrm{Law}(\theta\mid S_1,\dots,S_n)\) a.s.——每个人最终学得和看到了所有人的私有信息一样好。

证明骨架(Chebyshev 和不等式的等号条件 ⇒ 可测性)

令 \(Z_i=\log\frac{\mathbb P(\theta\in A\mid S_i)}{\mathbb P(\theta\in A^c\mid S_i)}\)、\(Z=\sum_i Z_i\),由 Bayes law 对 \(A\) 和 \(A^c\) 各用一次再相除消掉 \(f\),得 \(W:=\mathbb P(\theta\in A\mid S_1..S_n)=\frac{\beta e^Z}{1+\beta e^Z}\)。令 \(X:=\mathbb P(\theta\in A\mid H_i(\infty))\),由 tower property 有 \(X=\mathbb E[W\mid H_i(\infty)]=\mathbb E[W\mid X]\)。又因 \(Z_i\) 是 \(H_i(\infty)\)-可测的,

\[\mathbb E[Z_iW\mid X]=\mathbb E\big[Z_i\,\mathbb E[W\mid H_i(\infty)]\mid X\big]=\mathbb E[Z_i\mid X]\cdot\mathbb E[W\mid X]\]

对 \(i=1..n\) 求和得 \(\mathbb E[Z\cdot g(Z)\mid X]=\mathbb E[Z\mid X]\,\mathbb E[g(Z)\mid X]\),其中 \(g(x)=\frac{\beta e^x}{1+\beta e^x}\) 严格单调。这正好是 Chebyshev 和不等式取等号的条件,因此 \(Z\) 必须是 \(\sigma(X)\)-可测的,从而 \(W\) 也是,于是 \(X=W\) a.s.。证明思路继承 [MST11] 的 economics signal model,但那里是二元状态;这里要处理任意(可能无界支撑的)先验,所以才引入 Def 2.2/2.3。

4. Bernoulli sampling model:每轮只发 1 个 bit 也够

把通信压到极限:\(d=1\),\(\theta\sim\mathrm{Unif}([0,1])\),\(S_i\sim\mathrm{Ber}(\theta)\);每轮 agent 先抽 \(Y_i(t)\sim\nu_i(t)\),再抛一次 \(X_i(t)\sim\mathrm{Ber}(Y_i(t))\),只发出 \(X_i(t)\in\{0,1\}\)。既然只传 1 bit,能指望学到的最多是后验均值(它是这个模型的充分统计量)。

Theorem 2.9:后验均值 \(m_i(t)=\mathbb E[\theta\mid K_i(t)]\) 几乎必然且 \(L^1\) 收敛,所有 agent 的极限相同,且 \(m_i(\infty)=\mathbb E[\theta\mid S_1,\dots,S_n]=\frac{\sum_i S_i+1}{n+2}\)(由 Lemma 3.2:均匀先验 + \(m\) 次 Bernoulli 观测的后验是 \(\mathrm{Beta}(\sum Z_i+1,\,m+1-\sum Z_i)\))。

直观理解:额外那次抛硬币虽然再加了一层噪声,但它是无偏的:\(\mathbb E[X_j(t)\mid K_j(t)]=\mathbb E[Y_j(t)\mid K_j(t)]=m_j(t)\),且不同时刻的噪声互不相关。所以「长期平均邻居发来的 0/1」照样能无损还原邻居的后验均值,定理 2.1 的模仿论证原封不动可用。

5. 误规约版本:同一个模型退化成互相扔球的 Pólya urn

现在只改一件事:\(G\) 是一条边,agent 1 误以为 agent 2 发来的 \(\tilde X_2(s)\) 是 i.i.d. \(\mathrm{Ber}(\theta)\)(世界给的新独立样本),于是用 \(\tilde p(\theta\mid S_1,(\tilde X_2(s))_{s<t})\) 更新;agent 2 对称。

Proposition 2.10:\(\tilde\nu_i(t)\xrightarrow{w}\delta_{M(\infty)}\) a.s.(收敛到点质量,即后验方差 →0),两个 agent 的极限相同,但 \(\mathbb P\big(M(\infty)\neq\mathbb E[\theta\mid S_1,S_2]\big)=1\)——几乎必然学错,而且无限自信。

为什么?由 Lemma 3.2,误规约后验是 Beta 分布,均值 \(M_1(t)=\frac{S_1+\sum_{s<t}\tilde X_2(s)+1}{t+2}\)。所以 agent 1 这一轮发出的消息分布恰是 \(\mathrm{Ber}\big(\frac{S_1+\sum_{s<t}\tilde X_2(s)+1}{t+2}\big)\)——这完全就是一个 Pólya urn:罐子 \(i\) 初始装 \(S_i+1\) 个红球、\(2-S_i\) 个蓝球;每轮两个罐子各均匀摸一个球,把摸到的球送进对方罐子,如此重复。作者说文献里([Pem07] 综述、Bernard Friedman's urn [Fre65])找不到这个「互相喂球的双罐系统」,所以这个概率论结果本身可能独立有趣。

flowchart LR
  U1["urn 1: S1+1 red, 2-S1 blue"] -->|"draw uniformly"| B1["ball X1(t)"]
  B1 -->|"put into other urn"| U2["urn 2: S2+1 red, 2-S2 blue"]
  U2 -->|"draw uniformly"| B2["ball X2(t)"]
  B2 -->|"put into other urn"| U1
  U1 --> C["Beta posterior, variance shrinks to zero"]
  U2 --> C
  C --> R["random point mass M_inf, not Bayes optimal"]
          
drag to pan · scroll to zoom
误规约的 Bernoulli sampling model 等价于两个互相送球的 Pólya urn。每条被误当成「新证据」的消息都在强化已有的偏向,于是后验方差被虚假地压到 0,落点 \(M(\infty)\) 是个随机变量而非正确答案。
证明骨架(鞅 + \(L^2\) 递推 + urn 的无原子性)
  • 平均值是鞅:\(M(t)=\frac{M_1(t)+M_2(t)}2\)。由 \(\mathbb E[X_i(t)\mid\mathcal F(t)]=M_i(t)\) 得 \(\mathbb E[M_1(t+1)\mid\mathcal F(t)]=\frac{t+2}{t+3}M_1(t)+\frac1{t+3}M_2(t)\)(对称地有 \(M_2\)),两式相加得 \(\mathbb E[M(t+1)\mid\mathcal F(t)]=M(t)\)。有界鞅 ⇒ Doob 收敛到某随机变量 \(M(\infty)\in[0,1]\)。
  • 两人的分歧消失:\(D(t)=M_1(t)-M_2(t)\) 满足 \(a_{t+1}=\mathbb E[D(t+1)^2]\le\frac{t(t+2)}{(t+3)^2}a_t+\frac1{(t+3)^2}\)。迭代得 \((t+2)a_t\lesssim\log t\),即 \(a_t\lesssim\frac{\log t}{t}\to0\)(\(L^2\) 收敛)。再用超鞅 \(A(t)=D(t)^2-\sum_{i\le t}(i+2)^{-2}\) 升级到几乎必然收敛,得 \(D(t)\to0\) a.s.。思路取自 Bernard Friedman's urn 的分析 [Fre65]。
  • 方差塌缩:由 \(\tilde\nu_i(t)=\mathrm{Beta}\big((t+2)M_i(t),(t+2)(1-M_i(t))\big)\),方差 \(=\frac{M_i(t)(1-M_i(t))}{t+3}\to0\)。所以极限是点质量,不是分布——这就是「过度自信」的数学形式。
  • 几乎必然不等于最优:urn 模型的标准结论(time-dependent Pólya urn,[Pem90, Thm 3];也见 [AMR16])给出:给定 \(S_1,S_2\),\(M(\infty)\) 的条件分布没有原子。而 \(\mathbb E[\theta\mid S_1,S_2]\) 是 \(S_1,S_2\) 的确定性函数,所以两者相等的概率为 0。

6. 最关键的一句话

只要听者知道说话者是在从自己的后验里采样,重复采样的无偏性就让「听频率」等价于「直接交换后验」,于是一致性与最优学习都成立;一旦听者把这些样本误当成世界给的新独立证据,同一个过程就变成自我强化的 Pólya urn,收敛到一个随机的、确定的、错的答案。

主要实验结果

这是一篇纯理论论文(math.ST),没有 benchmark。「实验」有两类:§4.1 的解析反例对比,和 Appendix B.2 对高斯模型的数值模拟。本节另外补充一组本地复现的误规约 urn 模拟。

1. sampling model 严格优于有限行动的 action model(§4.1)

最小例子:\(G\) 是一条边,\(\theta\in\{0,1\}\) 均匀先验;信号 \(S_i\) 的密度在 \(\theta=1\) 时为 \(2x\)、在 \(\theta=0\) 时为 \(2-2x\)(\(x\in[0,1]\))。两个 agent 每轮只能发 0 或 1。

模型 每轮发什么 agent 1 的极限后验 \(\mathbb P(\theta=1\mid\cdot)\) 是否最优
action model 当前更可能的状态(取 argmax) \(\dfrac{6x_1}{2+4x_1}\) 否
sampling model(Bernoulli) 以当前后验为概率抛硬币 \(\dfrac{x_1x_2}{x_1x_2+(1-x_1)(1-x_2)}\) 是(= 真后验)

道理很简单:当 \(x_1,x_2>0.5\) 时,action model 里两人会永远发 1,于是 agent 1 关于 \(x_2\) 只学到「\(x_2>1/2\)」这一件事,后验卡在 \(\frac{6x_1}{2+4x_1}\) 不动;而 sampling model 里硬币会以后验为概率翻面,频率把 \(x_2\) 的确切信息慢慢漏出来,最终拼出真后验。作者指出只要行动空间有限,这个现象可以推广到其他网络和信号结构。

固定 \(x_2=0.75\),横轴是 agent 1 的私有信号 \(x_1\)。action model(红)系统性偏离真后验(蓝 = sampling model 的极限),而且偏差方向不固定——信号弱时过度自信,信号强时反而低估。两条线只在特殊点相交。

2. 高斯模型的数值模拟(Appendix B.2)

高斯情形下后验有闭式解(Lemma B.1):\(\nu_i(t)=\mathsf N(X_i(t),\sigma_i(t))\),其中 \(X_i(t)=\frac{\mu_i(t)^T\Sigma_{ii}(t)^{-1}H_i(t)}{\mu_i(t)^T\Sigma_{ii}(t)^{-1}\mu_i(t)+1}\),\(\sigma_i(t)=\frac1{\mu_i(t)^T\Sigma_{ii}(t)^{-1}\mu_i(t)+1}\),并可按 (30) 逐步递推,所以可以直接数值验证。

Figure 1 left
Figure 1(左):\(G\) 为一条边、信号强度 \(a_1=1,a_2=2\) 的一次实现,\(t\le500\)。真值 \(\theta\approx-0.502\),\(S_1\approx-0.371\)、\(S_2\approx-1.08\);虚线是最优贝叶斯后验均值 \(\frac{a_1S_1+a_2S_2}{a_1^2+a_2^2+1}\approx-0.423\)。纵轴尺度与信号方差相当。
Figure 1 right
Figure 1(右):同一次实现放大后看波动。两个 agent 的后验均值从 \(X_1(0)\approx-0.371\)、\(X_2(0)\approx-0.542\) 出发,\(t=500\) 时分别到 \(-0.406\) 与 \(-0.407\)——既互相一致,也贴住了最优值。采样噪声让轨迹一直抖动,但抖动幅度随 \(t\) 衰减。
Figure 2
Figure 2:\(n=7\)、\(a_i=i\) 时,完全图(clique,蓝)与环(cycle,橙)中第一个 agent 的后验方差,\(t\le20\)。虚线是最优后验方差 \((\sum_{i=1}^7 i^2+1)^{-1}\approx0.0071\)。方差单调下降(因为 \(X_i(t)\) 是鞅,方差 \(=\sigma_i(t)(1-\sigma_i(t))\)),且连接更密的网络收敛更快——符合直觉,但收敛速率的理论分析仍是开放问题。

3. 本地复现:误规约 urn 的「一致但错」

论文没给 Proposition 2.10 的数值图。这里按 §2.3 的递推直接模拟:取 \(S_1=1,S_2=0\),此时最优后验均值 \(\mathbb E[\theta\mid S_1,S_2]=\frac{1+1}{2+2}=0.5\)。

4 次独立运行中 agent 1 的误规约后验均值 \(M_1(t)\),\(t\le300\)。虚线是正确答案 0.5。每条轨迹都很快稳定下来——但稳定在哪里完全由早期几次抛硬币的运气决定(0.20 / 0.29 / 0.51 / 0.60),这正是 \(M(\infty)\) 为随机变量的含义。
20000 次运行(\(t=1500\))后极限值 \(M(\infty)\) 的分布直方图。它是一条散开的连续分布而不是集中在 0.5 的尖峰:落在 \([0.45,0.55]\) 的只有约 20%,而每一次运行的后验方差都已收缩到 \(O(1/t)\)。「一致 + 无原子 + 零方差」合起来就是过度自信的精确含义。
import random def misspecified_edge(S1, S2, T, rng): """§2.3: each agent treats the other's bit as a fresh i.i.d. Ber(theta).""" s1, s2 = S1, S2 # own signal + all received bits for t in range(1, T + 1): M1 = (s1 + 1) / (t + 2) # Beta posterior mean, Lemma 3.2 M2 = (s2 + 1) / (t + 2) X1 = 1 if rng.random() < M1 else 0 X2 = 1 if rng.random() < M2 else 0 s1 += X2 # "ball" sent from urn 2 into urn 1 s2 += X1 return M1, M2 # agree with each other, differ from 0.5 rng = random.Random(0) print(misspecified_edge(1, 0, 400, rng)) # -> (0.2139, 0.2040), optimal = 0.5

注意代码里唯一的错误就是分母 \(t+2\):agent 把 \(t-1\) 条「别人后验的采样」当成了 \(t-1\) 条独立观测。正确模型下这些消息的信息量会被折扣,分母不该这样涨。

4. 四个模型横向对比

模型 每轮传递 渐近一致? 学到最优后验? 出处
belief model(信念) 精确的后验均值 是 是 [Aum76, GP82, DS99]
action model(行动) 有限行动集中的 argmax 是 一般不是 [GK03, RSV09, MST14]
sampling model(本文) 后验的一个采样 \(Y_i(t)\) 是(Thm 2.1) 是(Thm 2.4) 本文
Bernoulli sampling(本文) 1 个 bit \(X_i(t)\sim\mathrm{Ber}(Y_i(t))\) 是(Thm 2.9) 是(后验均值层面) 本文
误规约 Bernoulli(本文) 同上,但听者误解生成机制 是(但收敛到点质量) 几乎必然不是 本文 Prop 2.10

局限与展望