利用近似贝叶斯计算估算巨寄生虫种群进化参数

我们使用近似贝叶斯计算（ABC）估算宿主内巨寄生虫种群的随机过程模型参数。宿主的免疫力是一个未观察到的模型变量，只有在宿主牺牲时成虫巨寄生虫被计数。在非常有限的数据条件下，过程速率可以相对精确地推断。建模涉及一个三变量的马尔可夫过程，其中观察到的数据似然计算在计算上是不可处理的。当似然在分析或计算上不可处理时，ABC方法特别有用。我们提出的ABC算法基于序贯蒙特卡洛，具有自适应特性，克服了以前ABC方法的一些缺点。该算法在一个涉及自回归模型模拟数据的测试示例上进行了验证，然后用于推断实验数据的马尔可夫过程模型参数。拟合的模型解释了由于一个短期存在于宿主中的零一免疫变量引起的观察到的额外超二项变异。