模拟LGCP kppm时如何指定固定数量的模拟点？

原生接口的现状说明 #

• ppm类对象的模拟支持通过调整rmhcontrol选项直接设置固定模拟点数，但目前没有找到LGCP对应的同类内置配置选项。
• 结合文档说明与源码逻辑，rLGCP的模拟点数生成规则已经明确：
  • 首先计算强度函数lambda与待模拟窗口面积的积分值，对应代码中的mu=integrate(lambda)调用逻辑，该规则与rpoispp的文档说明完全一致
  • 上述积分值会作为rpois的均值参数，通过nn=rpois(mu,nsims)的随机抽样结果，确定每一次模拟的最终点数

固定点数控制的实现方案 #

你提出的两类思路都是可落地的实现路径，并非仅能通过修改底层逻辑实现固定点数：

• 改造底层逻辑的直接控制方案：通过修改函数内部逻辑，将原本由lambda积分后泊松抽样得到的nn值替换为用户自定义常量，确实可以跳过随机抽点步骤，直接实现固定点数的LGCP模拟。这种方案没有冗余计算，生成效率更高，但需要自行封装修改后的模拟函数，无法直接调用原生的rLGCP接口。
• 随机稀疏的变通方案：你考虑的random thinning（随机稀疏）方案不需要改动原生代码，实现成本更低。具体操作流程为：先按照原生rLGCP的流程完成模拟，再对生成的点集做等概率无放回随机抽样，删减到目标点数即可。这种方式得到的点集会完整保留LGCP原有的空间相关性特征，不会破坏模拟的统计性质。

内容的提问来源于stack exchange，提问作者Eric Robinson