[论文解读] An algebra of Stein operators

Robert E. Gaunt, Guillaume Mijoule+1 more|arXiv (Cornell University)|Apr 22, 2016

Random Matrices and Applications参考文献 52被引用 4

一句话总结

本文提出了一种新颖的代数框架，用于构造乘积分布的多项式 Stein 算子，特别关注由独立随机变量乘积构成的随机变量。通过利用系数为多项式的微分算子并运用对偶性原理，作者推导出一大类分布（包括正态分布、贝塔分布、变方差伽马分布和广义伽马分布）的显式 k 阶 Stein 算子，即使目标密度未以闭式表达，也能实现可计算性。一项关键贡献是推导出 k 个独立对称变方差伽马分布随机变量乘积的密度的新公式。

ABSTRACT

We build upon recent advances on the distributional aspect of Stein's method to propose a novel and flexible technique for computing Stein operators for random variables that can be written as products of independent random variables. We show that our results are valid for a wide class of distributions including normal, beta, variance-gamma, generalized gamma and many more. Our operators are $k$th degree differential operators with polynomial coefficients; they are straightforward to obtain even when the target density bears no explicit handle. As an application, we derive a new formula for the density of the product of $k$ independent symmetric variance-gamma distributed random variables.

研究动机与目标

• 开发一种系统且灵活的代数框架，用于推导独立随机变量乘积分布的 Stein 算子。
• 将 Stein 方法的适用范围扩展至密度未显式给出或通过特殊函数定义的分布。
• 为一大类目标分布提供统一的方法，以构造具有多项式系数的 k 阶微分 Stein 算子。
• 推导出 k 个独立对称变方差伽马分布随机变量乘积密度的显式公式。
• 通过在经典和特殊分布（包括通过 Meijer G 函数和 Kummer U 函数定义的分布）中的应用，展示所提方法的可计算性与普适性。

提出的方法

• 作者将多项式 Stein 算子定义为乘法算子 M 与微分算子 D 的单项式有限线性组合，即 A = ∑ a_ij M^i D^j。
• 他们基于 Fokker-Planck 方程和分部积分法运用对偶性论证，从密度的微分方程推导算子，尤其当密度满足线性常微分方程时。
• 对于通过特殊函数（如 Kummer U、Meijer G）定义的分布，该方法利用这些函数已知的微分方程来构造相应的 Stein 算子。
• 通过将算子代数应用于独立分量，将该方法推广至乘积分布，利用联合密度可分解且算子结构尊重独立性的事实。
• 该方法包含一种归一化与缩放策略：对于正态分布和变方差伽马分布，算子在右侧乘以 M；对于 PRR 分布，乘以 M²。
• 通过推导经典分布（如正态分布、伽马分布、贝塔分布、F 分布、逆伽马分布、t 分布、PRR 分布、变方差伽马分布、广义伽马分布）的已知与新 Stein 算子，验证了该框架，结果与现有结论一致。

实验结果

研究问题

• RQ1 当目标密度未显式给出时，如何系统地构造独立随机变量乘积分布的多项式 Stein 算子？
• RQ2 具有多项式系数的 k 阶微分 Stein 算子在一大类分布中构造的代数结构是什么？
• RQ3 能否利用 Fokker-Planck 方程与 Stein 算子之间的对偶性，即使对于通过特殊函数定义的分布，也推导出可计算的算子？
• RQ4 k 个独立对称变方差伽马分布随机变量乘积的 Stein 算子形式为何？
• RQ5 所提方法如何将现有的一阶与二阶算子方法推广至高阶情形？

主要发现

• 本文成功构造了 k 个独立对称变方差伽马分布随机变量乘积的 k 阶多项式 Stein 算子，为该分布的密度提供了新的显式公式。
• 该方法为一大类分布（包括正态分布、贝塔分布、伽马分布、F 分布、逆伽马分布、t 分布、PRR 分布、变方差伽马分布和广义伽马分布）提供了一致且系统化的 Stein 算子推导方法，所有算子均具有显式的多项式系数微分算子形式。
• 对于 μ=0 的变方差伽马分布，Stein 算子为 σ²T₁Tᵣ + 2θMTᵣ/₂，其中 Tᵢ 表示映射 f ↦ x^i f’(x) 的算子。
• 该方法推广了已知结果：例如，经典正态 Stein 算子 Af(x) = f’(x) − xf(x) 可作为该通用框架的特例恢复。
• 该方法对通过特殊函数定义的分布具有鲁棒性：对于 PRR 分布，算子为 sT₁T₂ − M²T₂ₛ，其来源于 Kummer U 函数的微分方程。
• 该框架即使在目标密度未显式给出时，只要其满足线性微分方程，也能推导出 Stein 算子，展示了其在标准参数族之外的广泛应用性。