[论文解读] A J-function for inhomogeneous point processes

M. N. M. van Lieshout|arXiv (Cornell University)|Aug 26, 2010

Point processes and geometric inequalities被引用 4

一句话总结

本文提出了一种用于强度加权矩平稳点过程的新型非齐次 $J$-函数，将经典的 $J$-函数推广至考虑空间和时间非齐次性的场景。通过生成函数和条件强度表达该统计量，可在强度非平稳条件下检测空间相互作用（如抑制或聚集），并为泊松过程、对数高斯 Cox 过程及稀释硬核过程提供了显式计算公式，通过减法采样估计器和模拟实验进行了验证。

ABSTRACT

We propose new summary statistics for intensity-reweighted moment stationary point processes that generalise the well known J-, empty space, and nearest-neighbour distance distribution functions, represent them in terms of generating functionals and conditional intensities, and relate them to the inhomogeneous reduced second moment function. Extensions to space time and marked point processes are briefly discussed.

研究动机与目标

• 将经典的 $J$-函数推广至强度随空间或时间变化的非齐次点过程。
• 开发一种能捕捉相互作用结构（聚集或抑制）并同时考虑强度非平稳性的汇总统计量。
• 基于生成函数和条件强度，为新的 $J_{\text{inhom}}$ 函数提供理论框架。
• 推导出在关键模型（泊松过程、对数高斯 Cox 过程、稀释硬核过程）下 $J_{\text{inhom}}$ 的显式表达式。
• 提出一种减法采样估计器以实现实际估计，并通过模拟示例验证其性能。

提出的方法

• 利用 $n$-阶点过程相关函数和点过程的生成函数，定义非齐次 $J$-函数 $J_{\text{inhom}}(t)$。
• 通过条件强度和强度加权矩平稳性，表示 $J_{\text{inhom}}$。
• 通过两个生成函数的比值表达 $J_{\text{inhom}}$：$G_B^{!x}(1 - u_t^x)$ 与 $G(1 - u_t^x)$，其中 $u_t^x$ 编码空间距离和强度信息。
• 推导出三种模型类（泊松过程、对数高斯 Cox 过程、稀释硬核过程）下 $J_{\text{inhom}}$ 的显式形式。
• 提出一种 $K_{\text{inhom}}$ 的减法采样估计器，用于从观测数据中估计非齐次 $J$-函数。
• 通过调整相关结构并定义 $J^{ST}_{\text{inhom}}$ 和 $J^{B}_{\text{inhom}}$ 函数，在保持适当平移不变性假设的前提下，将框架扩展至时空点过程和标记点过程。

实验结果

研究问题

• RQ1 如何将经典 $J$-函数推广以适应点过程中的空间和时间非齐次性？
• RQ2 非齐次 $J$-函数在生成函数和条件强度方面的理论表达是什么？
• RQ3 新 $J_{\text{inhom}}$ 函数在不同点过程模型（如泊松过程、对数高斯 Cox 过程、稀释硬核过程）下的行为如何？
• RQ4 能否从观测数据中构建一个一致且实用的 $J_{\text{inhom}}$ 估计器，其在模拟中的表现如何？
• RQ5 如何在保持可解释性和统计一致性的同时，将该框架扩展至时空点过程和标记点过程？

主要发现

• 对于具有非齐次强度的泊松点过程，$J_{\text{inhom}}(t) = 1$，与预期一致，表明无相互作用。
• 对于具有底层高斯场正相关性的对数高斯 Cox 过程，$J_{\text{inhom}}(t)$ 低于泊松参考值，表明在考虑非齐次性后存在聚集现象。
• 对于稀释硬核过程，$J_{\text{inhom}}(t)$ 在初始段呈平坦状态且低于泊松参考值，直至 $r \approx 0.2$，反映出由于硬核距离导致的抑制。
• 对 $K_{\text{inhom}}$ 的减法采样估计器能正确识别出与泊松模型的偏离：在对数高斯 Cox 情况下，$t > 0.13$ 时估计值低于 $\pi t^2$。
• 在稀释硬核过程中，估计器通过在 $t \approx 0.2$ 以内显示低于泊松参考值的结果，证实了抑制现象，与 $J_{\text{inhom}}$-函数的行为一致。
• 将框架扩展至时空和标记点过程是可行的：当截断至 $n=1$ 时，$J^{ST}_{\text{inhom}}$ 退化为 Gabriel 和 Diggle 的 $K^*$ 方法；$J^{B}_{\text{inhom}}$ 允许在给定标记集条件下的相互作用分析。