[论文解读] A Monad for Probabilistic Point Processes

Swaraj Dash, Sam Staton|arXiv (Cornell University)|Jan 26, 2021

Point processes and geometric inequalities被引用 4

一句话总结

本文通过组合Giry函子$G$（用于概率测度）与袋函子$B$（用于多重集），引入了用于概率点过程的单子$GB = G \circ B$，并利用单子之间的分配律证明其构成单子。关键贡献在于证明了从点过程到其期望计数的强度测度映射是一个单子态射，从而实现了强度测度的组合式计算，并使Wald引理由成为其特例。

ABSTRACT

A point process on a space is a random bag of elements of that space. In this paper we explore programming with point processes in a monadic style. To this end we identify point processes on a space X with probability measures of bags of elements in X. We describe this view of point processes using the composition of the Giry and bag monads on the category of measurable spaces and functions and prove that this composition also forms a monad using a distributive law for monads. Finally, we define a morphism from a point process to its intensity measure, and show that this is a monad morphism. A special case of this monad morphism gives us Wald's Lemma, an identity used to calculate the expected value of the sum of a random number of random variables. Using our monad we define a range of point processes and point process operations and compositionally compute their corresponding intensity measures using the monad morphism.

研究动机与目标

• 为使用单子建模概率点过程提供范畴论的、组合式的框架。
• 将点过程形式化为可测空间$X$上多重集（袋）的概率测度，即$G(B(X))$中的元素。
• 确立点过程的强度测度是一个单子态射，支持组合式推理。
• 通过单子态射的性质，证明该框架可作为Wald引理由的特例。
• 在概率编程与数据库等实际统计建模中，弥合范畴论与实践应用之间的鸿沟。

提出的方法

• 将点过程形式化为$G(B(X))$中的元素，即在可测空间$X$上对袋函子应用Giry函子。
• 利用Beck的单子间分配律理论，证明$GB$构成一个单子。
• 将强度测度定义为复合态射$GB \hookrightarrow MM \xrightarrow{\mu} M$，其中$M$是所有测度的单子。
• 通过验证相关图表的交换性，证明强度测度映射$\mathbb{E}: GB \to M$是一个单子态射。
• 利用单子态射的性质，以函数式编程风格（如do-notation）实现强度测度的组合式计算。
• 利用自然性与函子性，推导出复合分布的Wald引理等恒等式。

实验结果

研究问题

• RQ1 能否在范畴论中系统地将点过程建模为单子，以支持组合式构造？
• RQ2 是否存在从点过程到其强度测度的自然单子态射，保持单子结构？
• RQ3 Wald引理由关于随机和期望的结论，能否作为单子态射性质的推论？
• RQ4 如何利用单子结构组合式地计算复合点过程的强度测度？
• RQ5 该范畴论框架能否应用于概率逻辑编程与概率数据库，以支持结构化建模？

主要发现

• 复合函子$GB = G \circ B$在可测空间范畴上构成一个单子，为概率点过程提供了范畴论基础。
• 强度测度映射$\mathbb{E}: GB \to M$是一个单子态射，确保强度计算与单子组合可交换。
• Wald引理由作为单子态射性质的特例被恢复，表明$\mathbb{E}(\underline{N} \gg t. \underline{X}) = \mathbb{E}(\underline{N}) \gg t. \mathbb{E}(\underline{X})$。
• 对于单位正方形上强度为$\Lambda$的泊松点过程，强度测度为$\lambda W \mapsto \Lambda \cdot |W|$，可通过单子态射组合式计算。
• 该框架支持使用函数式编程语法（如do-notation）对复杂点过程进行组合式模拟与分析。
• 该结果提供了一种基于范畴论的、原则性的、组合式方法，用于计算强度测度，无需手动积分。