Florian Puchhammer公司;阿马尔·本·阿卜杜拉;L’Ecuyer，皮埃尔 随机生物和化学反应网络的tau-leaping模拟的array-RQMC方差减少。 （英语） Zbl 1470.92125号 牛市。数学。生物。 83，第8号，第91号论文，第31页（2021年）. 小结：我们探讨了使用Array-RQMC（一种用于模拟马尔可夫链的随机拟蒙特卡罗方法）来减少使用跳变模拟随机生物或化学反应网络时的方差。任务是通过样本路径上的样本平均值估计给定未来时间（T）的分子拷贝数函数的期望值，目标是减少该样本平均值估计器的方差。我们发现，当该方法得到适当应用时，可以获得数千倍的方差减少。这些因素远大于其他作者以前对相同示例尝试RQMC方法时观察到的因素。Array-RQMC模拟马尔可夫链的一系列实现，并需要一个排序函数来根据这些链的状态在每一步之后重新排序。排序函数的选择是该方法效率的关键因素，尽管在我们的实验中，无论排序方法如何，Array RQMC都从未比普通蒙特卡洛差过。每一步每种类型的预期反应次数也会对效率增益产生影响。 引用于2文件 理学硕士： 92立方厘米 系统生物学、网络 92C40型 生物化学、分子生物学 92E20型 化学中的经典流动、反应等 60J22型 马尔可夫链中的计算方法 65二氧化碳 蒙特卡罗方法 关键词：化学反应网络;随机生物系统;方差减少;准蒙特卡罗;阵列-RQMC;牛叉;连续时间马尔可夫链;吉莱斯皮 软件：晶格生成器;MRG32k3a型;RandQMC公司;SSJ公司;索波尔.cc;算法955 PDF格式BibTeX公司 XML格式引用 \textit{F.Puchhammer}等人，公牛。数学。生物学83，第8期，论文91，31页（2021；Zbl 1470.92125） 全文： 内政部 arXiv公司 链接 参考文献： [1] Anderson，WJ，《连续时间马尔可夫链：面向应用的方法》（1991），纽约：Springer，纽约·Zbl 0731.60067号 ·doi:10.1007/978-1-4612-3038-0 [2] Anderson，DF，将跳跃后检查纳入tau-leaping，《化学物理杂志》，128，5，054-103（2008）·doi:10.1063/1.2819665 [3] 安德森，D。；Higham，D.，《连续时间马尔可夫链的多级蒙特卡罗及其在生化动力学中的应用》，多尺度模型模拟，10，1，146-179（2012）·Zbl 1262.60072号 ·数字对象标识代码：10.1137/10840546 [4] 安德森，DF；Kurtz，TG；Koepl，H.等人。；Densmore，D。；塞蒂，G。；di Bernardo，M.，化学反应网络的连续时间马尔可夫链模型，生物分子电路的设计与分析，3-42（2011），纽约：Springer，纽约·doi:10.1007/978-1-4419-6766-4_1 [5] 比恩特杰斯，CHL；Baker，RE，应用于随机生物系统中τ跳跃的准蒙特卡罗方法，Bull Math Biol，812931-2959（2019）·Zbl 1422.92066号 ·doi:10.1007/s11538-018-0442-2 [6] Ben Abdelah A，L’Ecuyer P，Puchhammer F（2019）随机波动模型下期权定价的阵列RQMC。摘自：2019年冬季模拟会议记录。IEEE出版社，第440-451页。https://www.informs-sim.org/wsc19papers/429.pdf [7] Breiman，L.，《随机森林》，《马赫学习》，第45、1、5-32页（2001年）·Zbl 1007.68152号 ·doi:10.1023/A:1010933404324 [8] 曹毅。；Gillespie，DT；Petzold，LR，《慢尺度随机模拟算法》，《化学物理杂志》，122,1014-116（2005）·数字对象标识代码：10.1063/1.1824902 [9] de Boor，C.，《样条曲线实用指南》（2001），纽约：Springer，纽约·Zbl 0987.65015号 [10] Demers V，L'Ecuyer P，Tuffin B（2005）随机准蒙特卡罗与分裂相结合的稀土模拟。摘自：2005年欧洲模拟和建模会议记录。欧洲，比利时根特，第25-32页 [11] Dick，J。；Pillichshammer，F.，《数字网络和序列：差异理论和准蒙特卡罗积分》（2010），剑桥：剑桥大学出版社，剑桥·Zbl 1282.65012号 ·doi:10.1017/CBO9780511761188 [12] Dick，J。；伊利诺伊州斯隆；王，X。；Woźniakowski，H.，具有一般权重的加权Korobov空间中的好格规则，数值数学，103，63-97（2006）·Zbl 1097.65004号 ·doi:10.1007/s00211-005-0674-6 [13] Dion M，L'Ecuyer P（2010）随机准蒙特卡罗模拟的美式期权定价。摘自：2010年冬季模拟会议记录，第2705-2720页 [14] El Haddad，R。；莱科特，C。；L'Ecuyer，P。；凯勒，A。；Heinrich，S。；Niederreiter，H.，多维状态空间上离散时间马尔可夫链的准蒙特卡罗模拟，蒙特卡罗和准蒙特卡洛方法2006，413-429（2008），柏林：Springer，柏林·Zbl 1141.65311号 ·doi:10.1007/978-3-540-74496-2_24 [15] 福克斯，BL；Glynn，PW，模拟有限时域马尔可夫过程的离散时间转换，SIAM J Appl Math，501457-1473（1990）·Zbl 0825.60043号 ·doi:10.1137/0150085 [16] Gerber，M。；肖邦，N.，《序贯准蒙特卡洛》，J R Stat Soc Ser B，77，第3部分，509-579（2015）·Zbl 1414.62109号 ·doi:10.1111/rssb.12104 [17] Giles，MB，算法955:逆泊松累积分布的近似，ACM Trans Math Softw，42，1-22（2016）·Zbl 1347.65019号 ·doi:10.145/2699466 [18] Gillespie，DT，耦合化学反应的精确随机模拟，《物理化学杂志》，81，25，2340-2361（1977）·doi:10.1021/j100540a008 [19] Gillespie，DT，化学Langevin方程，化学物理杂志，113，1，297-306（2000）·doi:10.1063/1.481811 [20] Gillespie，DT，化学反应系统的近似加速随机模拟，化学物理杂志，115，4，1716-1733（2001）·数字对象标识代码：10.1063/1.1378322 [21] Hellander，A.，基于均匀化和准蒙特卡罗的化学主方程瞬态解的有效计算，化学物理杂志，128，154，109（2008） [22] FJ希克内尔；Hellekalek，P。；Larcher，G.，《莱迪思规则：它们的衡量标准如何？》？，随机和准随机点集，109-166（1998），纽约：Springer，纽约·Zbl 0920.65010号 ·doi:10.1007/978-1-4612-1702-23 [23] FJ希克内尔；方，KT；FJ希克内尔；尼德雷特，H.，获得格求积规则的（{O（N^{-2+\epsilon}）}收敛性，蒙特卡罗和准蒙特卡罗方法2000，274-289（2002），柏林：斯普林格，柏林·Zbl 1002.65009号 ·数字对象标识代码：10.1007/978-3-642-56046-0_18 [24] FJ希克内尔；Hong，HS；L'Ecuyer，P。；Lemieux，C.，拟蒙特卡罗求积的可扩展格序列，SIAM科学计算杂志，22，3，1117-1138（2001）·Zbl 0974.65004号 ·doi:10.1137/S1064827599356638 [25] Higham，DJ，化学反应建模与模拟，SIAM Rev，50，2，347-368（2008）·Zbl 1144.80011号 ·数字对象标识代码：10.1137/060666457 [26] Joe，S。；Kuo，FY，用更好的二维投影构建Sobol序列，SIAM科学计算杂志，30，5，2635-2654（2008）·Zbl 1171.65364号 ·doi:10.1137/070709359 [27] 金，JK；Josić，K。；Bennett，MR，《随机和确定性准静态近似之间的关系》，BMC系统生物学，9，87，1-13（2015）·doi:10.1186/s12918-015-0218-3 [28] 科洛登，PE；Platen，E.，随机微分方程的数值解（1992），柏林：Springer，柏林·Zbl 0752.60043号 ·doi:10.1007/978-3-662-12616-5 [29] Koh，W。；Blackwell，KT，《改进的基于梯度扩散的空间直接法以保持完全扩散波动》，《化学物理杂志》，137，15，154-111（2012）·doi:10.1063/1.4758459 [30] 莱科特，C。；Coulibaly，I.，《使用线性Boltzmann方程网络的准蒙特卡罗方法》，SIAM J Numer Ana，35，1，51-70（1998）·Zbl 0913.65137号 ·doi:10.137/S0036142995290051 [31] L'Ecuyer，P.，组合多递归随机数生成器的良好参数和实现，Oper Res，47，1，159-164（1999）·Zbl 1042.65505号 ·doi:10.1287/opre.47.1159 [32] L'Ecuyer，P.，金融中的准蒙特卡罗方法及其应用，Fin Stoch，13，3，307-349（2009）·Zbl 1199.65004号 ·doi:10.1007/s00780-009-0095-y [33] L'Ecuyer，P。；温柔点，JE；Haerdle，W。；Mori，Y.，《随机数生成》，《计算统计手册》，35-71（2012），柏林：施普林格出版社，柏林·Zbl 1243.62001 ·doi:10.1007/978-3-642-21551-33 [34] L'Ecuyer P（2016）SSJ：Java中的随机模拟。http://simul.iro.umontreal.ca/ssj/ [35] L'Ecuyer，P。；Glynn，普华永道；Owen，AB，《随机拟蒙特卡罗：从业者简介》，蒙特卡罗和拟蒙特卡洛方法：MCQMC 2016，29-52（2018），柏林：施普林格，柏林·Zbl 1407.65013号 ·doi:10.1007/978-3-319-91436-7_2 [36] L'Ecuyer P，Build E（2005）《Java与SSJ的仿真》。摘自：2005年冬季模拟会议记录。新泽西州皮斯卡塔韦IEEE出版社，第611-620页 [37] L'Ecuyer，P。；Lemieux，C.，通过格规则减少方差，《管理科学》，46，9，1214-1235（2000）·Zbl 1232.65008号 ·doi:10.1287/mnsc.46.9.1214.12231 [38] L'Ecuyer，P。；Lemieux，C。；德罗，M。；L'Ecuyer，P。；Szidarovszky，F.，《随机拟蒙特卡罗方法的最新进展，建模不确定性：随机理论、方法和应用的检验》，419-474（2002），波士顿：Kluwer学术出版社，波士顿·Zbl 0979.00016号 ·数字对象标识代码：10.1007/0-306-48102-220 [39] L'Ecuyer，P。；芒格，D。；沃罗日·尼亚科夫斯基，H。；Plaskota，L.，《关于随机移动晶格规则的优值》，蒙特卡罗和准蒙特卡罗方法2010，133-159（2012），柏林：施普林格，柏林·Zbl 1271.65045号 ·doi:10.1007/978-3-642-27440-46 [40] L'Ecuyer P，Munger D（2016）Algorithm 958:lattice builder:构建秩1格规则的通用软件工具。ACM Trans Mathe Softw 42（2）：第15条 [41] L’Ecuyer，P。；莱科特，C。；Tuffin，B。；尼德雷特，H。；Talay，D.，具有有序状态空间的马尔可夫链的随机准蒙特卡罗模拟，蒙特卡罗和准蒙特卡罗方法2004，331-342（2006），柏林：施普林格，柏林·Zbl 1097.65016号 ·doi:10.1007/3-540-31186-619 [42] L'Ecuyer P、Demers V、Tuffin B（2007），《拉里火山、分裂和准蒙特卡洛》。ACM Trans Model Comput Simul 17（2）：第9条，45页·Zbl 1281.62085号 [43] L'Ecuyer，P。；莱科特，C。；Tuffin，B.，马尔可夫链的随机准蒙特卡罗模拟方法，Oper Res，56，4，958-975（2008）·Zbl 1167.90673号 ·doi:10.1287/opre.1080.0556 [44] L'Ecuyer，P。；莱科特，C。；阿尔切夫·古德特，A。；L'Ecuyer，P。；Owen，AB，《马尔可夫链的数组上RQMC：映射备选方案和收敛速度》，蒙特卡罗和准蒙特卡罗方法2008，485-500（2009），柏林：施普林格，柏林·Zbl 1184.65018号 ·doi:10.1007/978-3642-04107-5_31 [45] L'Ecuyer，P。；芒格，D。；莱科特，C。；Tuffin，B.，Array-RQMC的排序方法和收敛速度：一些经验比较，数学计算模拟，143191-201（2018）·Zbl 1482.65005号 ·doi:10.1016/j.matcom.2016.07.010 [46] L'Ecuyer P、Marion P、Godin M、Fuchhammer F（2020）QMC和RQMC点集自定义构建工具。在：蒙特卡罗和准蒙特卡罗方法：MCQMC 2020，提交了手稿。可在网址：http://www.iro.umontreal.ca/lecuyer/myftp/papers/mcqmc20latnet.pdf [47] Lemieux，C.，蒙特卡罗和准蒙特卡罗抽样（2009），柏林：施普林格，柏林·Zbl 1269.65001号 [48] Lemieux C，Cieslak M，Luttmer K（2004）RandQMC用户指南：C中随机准蒙特卡罗方法的软件包。软件用户指南。可在http://www.math.uwaterloo.ca网站/clemieux/randqmc.html [49] Matousěk，J.，关于锚固箱的（L_2）差异，J Complex，14527-556（1998）·Zbl 0942.65021号 [50] Niederreiter H（1992）随机数生成和拟蒙特卡罗方法。摘自：SIAM CBMS-NSF应用数学区域会议系列，第63卷。SIAM公司·兹比尔0761.65002 [51] Owen，AB，扰乱均匀分布求积的蒙特卡罗方差，SIAM J Numer Ana，34，5，1884-1910（1997）·Zbl 0890.65023号 ·doi:10.137/S0036142994277468 [52] Owen，AB，光滑函数积分的加扰净方差，Ann Stat，25，4，1541-1562（1997）·Zbl 0886.65018号 ·doi:10.1214/aos/1031594731 [53] Owen，AB，用于超高维模拟的拉丁超立方体采样，ACM Trans Model Comput Simul，8，1，71-102（1998）·Zbl 0917.65022号 ·数字对象标识代码：10.1145/272991.273010 [54] 帕吉特，JMA；Ilie，S.，反应扩散系统随机模拟的自适应τ叉方法，AIP Adv，6，3，035-217（2016）·doi:10.1063/1.4944952 [55] Pollock DSG（1993）使用三次样条曲线进行平滑。伦敦大学玛丽女王学院和伦敦韦斯特菲尔德学院技术报告 [56] Rao，CV；阿金，AP，《随机化学动力学和准静态假设：Gillespie算法的应用》，《化学物理杂志》，118，11，4999-5010（2003）·数字对象标识代码：10.1063/1.1545446 [57] 西内斯库，V。；L'Ecuyer，P.，随机移位格规则的方差界和存在性结果，计算应用数学杂志，2363296-3307（2012）·Zbl 1264.65005号 ·doi:10.1016/j.cam.2012.02.034 [58] 伊利诺伊州斯隆；Joe，S.，《多重积分的晶格方法》（1994），牛津：克拉伦登出版社，牛津·Zbl 0855.65013号 [59] Sobol’，IM，立方体中点的分布和积分的近似计算，苏联计算数学物理，7，4，86-112（1967）·Zbl 0185.41103号 ·doi:10.1016/0041-5553（67）90144-9 [60] 斯特雷尔，R。；Ilie，S.，慢动力学和快动力学反应扩散系统的混合随机模拟，化学物理杂志，143，23，108-234（2015）·数字对象标识代码：10.1063/1.4937491 [61] 托马斯·P。；斯特劳贝，AV；Grima，R.，《慢尺度线性噪声近似：时间尺度分离条件下生物化学网络的精确、简化随机描述》，BMC系统生物学（2012）·doi:10.1186/1752-0509-6-39 [62] 瓦希特，C。；凯勒，A。；凯勒，A。；Heinrich，S。；Niederreiter，H.，《马尔可夫链的高效同时模拟》，蒙特卡罗和准蒙特卡罗方法，2006年，669-684（2008年），柏林：施普林格出版社，柏林·Zbl 1141.65325号 ·doi:10.1007/978-3-540-74496-2_41 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。