[论文解读] Efficient simulation of individual-based population models: the R Package IBMPopSim

Daphné Giorgi, Sarah Kaakaï+1 more|arXiv (Cornell University)|Mar 10, 2023

demographic modeling and climate adaptationDecision Sciences参考文献 26被引用 3

一句话总结

该论文介绍了 IBMPopSim，一个 R 包，通过 Rcpp 实现的 C++ 加速分层算法，实现了对随机个体基础种群模型（IBMs）的高效、精确模拟。该工具支持具有年龄结构的个体、动态事件强度及相互作用的复杂种群动态，为精算学、生态学和人口学建模提供了可扩展的解决方案，具有经过验证的计算效率和通过编译 C++ 代码实现的用户友好型 R 集成。

ABSTRACT

The R Package IBMPopSim aims to simulate the random evolution of heterogeneous populations using stochastic Individual-Based Models (IBMs). The package enables users to simulate population evolution, in which individuals are characterized by their age and some characteristics, and the population is modified by different types of events, including births/arrivals, death/exit events, or changes of characteristics. The frequency at which an event can occur to an individual can depend on their age and characteristics, but also on the characteristics of other individuals (interactions). Such models have a wide range of applications in fields including actuarial science, biology, ecology or epidemiology. IBMPopSim overcomes the limitations of time-consuming IBMs simulations by implementing new efficient algorithms based on thinning methods, which are compiled using the Rcpp package while providing a user-friendly interface.

研究动机与目标

• 解决传统个体基础模型（IBM）在具有动态事件强度的大规模异质种群中计算效率低下的问题。
• 为连续时间中的随机 IBM 提供可扩展、高性能的模拟框架，特别适用于精算学和人口学应用。
• 将先进的随机模拟技术——特别是泊松随机测度和分层方法——集成到用户友好的 R 包中，通过 Rcpp 实现无缝的 C++ 编译。
• 在单一、可扩展的计算环境中支持复杂的种群动态，包括个体相互作用、与年龄相关的事件发生率以及状态依赖的转移。
• 使研究人员能够以极少的代码定义自定义事件和强度，同时通过跳跃过程的路径表示确保精确模拟。

提出的方法

• 使用泊松随机测度对 IBM 进行路径表示，将事件发生建模为连续时间马尔可夫过程中的跳跃。
• 通过用更简单、有界均匀的强度过程支配原始强度，应用分层方法精确模拟事件时间。
• 利用 Rcpp 编译事件生成和强度评估的 C++ 实现，实现高性能计算，同时保持 R 层的易用性。
• 通过三种事件类型建模种群动态：出生/进入、死亡/退出以及状态变化（如特征转移），每类事件的强度函数均依赖于个体年龄、特征及种群层面的相互作用。
• 实现多变量计数过程框架，通过所有潜在事件之间的竞争性分层确定事件时间，确保精确模拟，无需拒绝采样。
• 通过允许强度函数依赖于多个个体的联合状态，支持非相互作用和相互作用种群，且通过支配论证提供了存在性和非爆炸性的形式证明。

实验结果

研究问题

• RQ1 如何在连续时间中高效模拟具有复杂、状态依赖事件强度和相互作用的个体基础模型？
• RQ2 哪些计算策略能够在不依赖离散时间近似或拒绝采样法的前提下，实现对大规模异质种群的精确模拟？
• RQ3 通过 Rcpp 集成 C++ 能在多大程度上提升结构化种群随机模拟的性能，同时保持 R 层的易用性？
• RQ4 如何系统性地将分层方法应用于具有年龄和相互作用依赖强度的一般类个体基础模型？
• RQ5 在存在种群层面相互作用的情况下，确保模拟算法非爆炸性和路径正确性的理论基础是什么？

主要发现

• IBMPopSim 包通过 Rcpp 编译的基于分层的算法实现了对个体基础模型的精确模拟，与朴素实现相比显著减少了运行时间。
• 该方法通过与更简单计数过程的路径支配确保了过程的非爆炸性，理论依据通过 Feller 条件验证提供。
• 该包支持复杂的动态，包括与年龄相关和与相互作用相关的事件强度，能够实现对保险组合和人口系统的现实建模。
• 通过泊松随机测度的路径表示，可在一般状态空间中精确模拟跳跃过程，避免了时间离散化带来的偏差。
• 通过 Rcpp 集成 C++ 代码，实现了高性能计算，同时保持了高层、用户友好的 R 接口，用于定义事件和强度。
• 理论结果证实，模拟过程与原始 IBM 具有相同的概率分布，确保了模拟框架的正确性。