×

使用随机模型的近似连续极限描述进行可靠有效的参数估计。 (英语) Zbl 1497.92027号

摘要:基于随机个体的数学模型对于建模生物现象很有吸引力,因为它们自然地捕捉到了生物数据中经常明显的随机性和可变性。这种模型还允许我们跟踪感兴趣人群中个体的运动。不幸的是,捕捉这种微观细节意味着模拟和参数推断可能会变得计算昂贵。克服这一计算限制的一种方法是粗粒度随机模型,以提供一个近似的连续模型,该模型可以用少得多的计算工作量求解。然而,粗粒度连续体模型可能有偏差或不准确,特别是对于某些参数状态。在这项工作中,我们通过演示如何模拟两种常用的实验,将随机数学模型和连续数学模型结合在二维细胞生物学实验的基于格子的模型的背景下:细胞增殖试验屏障测定我们的方法包括建立一个简单的统计模型,以反映昂贵的随机模型和相关的计算成本低廉的粗粒度连续模型之间的差异。我们基于从用户选择分布中抽取的设计点上的有限数量的昂贵随机模型评估,形成了这个统计模型,对应于一个计算机实验设计问题。通过简单的设计点选择方案,我们表明,将差异的统计模型与计算成本低廉的连续体模型结合使用,可以进行预测和推理,同时纠正由于连续体近似导致的偏差和不准确。我们通过模拟增殖试验,其中连续极限模型是众所周知的logistic常微分方程,以及屏障测定其中连续极限模型与著名的Fisher-KPP偏微分方程密切相关。我们构造了一个近似似然函数来进行参数推断,无论有无差异校正项。使用最大似然估计,我们提供了未知参数的点估计,并使用剖面似然来描述这些估计中的不确定性,并形成近似置信区间。对于所考虑的推理问题范围,仅使用连续极限模型会导致有偏差的参数估计和覆盖率较差的置信区间。相反,结合模型差异统计模型中产生的校正项,我们可以以最小的计算开销准确地恢复参数。主要的折衷是,相关的置信区间通常更宽,反映了近似过程引入的额外不确定性。复制这项工作结果所需的所有算法都是用开源Julia语言编写的,可以在GitHub上获得。

MSC公司:

92立方厘米 发育生物学,模式形成
92立方37 细胞生物学
60G99型 随机过程
PDF格式BibTeX公司 XML格式引用
全文: 内政部 链接

参考文献:

[1] 贝克·R·E。;Simpson,M.J.,《纠正出生-死亡运动过程的平均场近似值》,Phys。E版,82,第041905条pp.(2010)
[2] 宾尼,R.N。;A.詹姆斯。;普朗克,M.J.,《包含邻域相关方向偏差的集体细胞运动的空间力矩动力学》,J.R.Soc.Interface,120150228(2015)
[3] Bolker,B。;Pacala,S.W.,《利用矩方程来理解生态系统中随机驱动的空间格局形成》,Theor。大众。生物学,52,179-197(1997)·兹伯利0890.92020
[4] Bosco,D.B。;肯沃西,R。;Zorio,D.A.R。;Sang,Q.-X.A.,人类间充质干细胞通过采用非增殖性成纤维细胞状态对紫杉醇产生耐药性,PLOS ONE,10,文章e0128511 pp.(2015)
[5] 布朗宁,A.P。;McCue,S.W。;Simpson,M.J.,探索细胞增殖试验最佳持续时间的贝叶斯计算方法,公牛。数学。生物学,79,1888-1906(2017)·兹比尔1373.92004
[6] 布鲁纳,M。;Chapman,S.J.,《硬球扩散中的排除体积效应》,《物理学》。E版,85,第011103条pp.(2012)
[7] 考克斯·D·R。;Wermuth,N.,《多元相关性:模型、分析和解释》(1996),Chapman和Hall/CRC·Zbl 0880.62124号
[8] Decaestecker,C。;O.Debeir。;Van Ham,P。;Kiss,R.,抗迁移药物可以在体外筛选吗?细胞迁移定量分析的2D和3D分析综述,《医学评论研究》,27,149-176(2006)
[9] 德鲁尔斯,C。;奥伯特,M。;巴杜阿尔,M。;Grammaticos,B.,《肿瘤细胞迁移建模:从微观模型到宏观模型》,Phys。E版,79,第031917条pp.(2009)
[10] Eaton,M.L.,《多元统计:向量空间方法》,《数学研究所》。统计讲座笔记专题。序列号。,53 (2007) ·Zbl 1160.62326号
[11] Fisher,R.A.,《优势基因的发展浪潮》,《优生学年鉴》,第7355-369页(1937年)
[12] 弗雷斯特,A。;Sobester,A。;Keane,A.,《通过代理建模进行工程设计:实用指南》(2008),John Wiley&Sons
[13] Frasca,M。;Sharkey,K.J.,《相互作用个体群体中流行病传播的离散时间矩闭合》,J.Theor。生物学,399,13-21(2016)·Zbl 1347.92092号
[14] Gillespie,D.T.,耦合化学反应的精确随机模拟,化学杂志。物理。,81, 2340-2361 (1977)
[15] Gramacy,R.B.,《替代品:应用科学的高斯过程建模、设计和优化》(2020年),查普曼和霍尔/CRC
[16] 哈里达斯,P。;彭宁顿,C.J。;麦戈文,J.A。;McElwain,D.L.S。;Simpson,M.J.,《共同培养屏障分析中细胞迁移和细胞增殖速率的量化》,揭示了皮肤和黑色素瘤细胞在黑色素瘤扩散和侵袭过程中的相互作用,J.Theor。生物学,423,13-25(2017)
[17] 海因斯,K.E。;Middendorf,T.R。;Aldrich,R.W.,《非线性生物物理模型中参数可识别性的测定:贝叶斯方法》,《普通物理杂志》。,143, 401-416 (2014)
[18] 姜强。;Fu,X。;Yan,S。;李,R。;杜,W。;曹,Z。;钱,F。;Grima,R.,基因表达的非马尔可夫模型的神经网络辅助近似和参数推断,Nat.Commun。,12, 2618 (2021)
[19] Johnson,S.G.,2021年。Julia的NLopt模块。2022年2月检索到NLopt。
[20] 凯皮奥,J。;Kolehmainen,V.,逆向问题中建模误差和不确定性的近似边缘化,《贝叶斯理论与应用》一书,644-672(2013)
[21] 凯皮奥,J。;Somersalo,E.,《统计与计算反问题》,第160卷(2006),Springer·Zbl 0927.35134号
[22] 凯皮奥,J。;Somersalo,E.,《统计反问题:离散化、模型简化和反犯罪》,J.Compute。应用。数学。,198, 493-504 (2007) ·Zbl 1101.65008号
[23] M.C.肯尼迪。;O'Hagan,A.,《当快速近似可用时预测复杂计算机代码的输出》,Biometrika,87,1-13(2000)·Zbl 0974.62024号
[24] M.C.肯尼迪。;O'Hagan,A.,《计算机模型的贝叶斯校准》,J.R.Stat.Soc.Ser。B: 统计方法。,63, 425-464 (2001) ·Zbl 1007.62021号
[25] Koehler,J.R。;Owen,A.B.,《计算机实验》,《手册统计》,第13期,第261-308页(1996年)·兹比尔0919.62089
[26] 科尔莫戈罗夫,A.N。;彼得罗夫斯基,P.G。;Piskunov,N.S.,《物质数量增加时扩散方程的研究及其在生物问题中的应用》,莫斯科大学数学系。公牛。,1, 1-26 (1937)
[27] 法律,R。;D.J.Murrell。;Dieckmann,U.,《时空中的人口增长:空间逻辑方程》,生态学,84(2003),22-22.www.jstor.org/stable/3108013
[28] Lei,C.L。;Ghosh,S。;惠特克,D.G。;Aboelkassem,Y。;Beattie,K.A。;坎特韦尔,C.D。;Delhaas,T。;休斯顿,C。;Montes Novaes,G。;潘菲洛夫,A.V。;Pathmanathan,P。;Riabiz,M。;韦伯·多斯·桑托斯,R。;Walmsley,J。;Worden,K。;米拉姆斯,G.R。;Wilkinson,R.D.,在校准机械电生理模型时考虑差异,Philos。事务处理。R.Soc.A,37820190349(2020)·Zbl 1470.92057号
[29] Liang,C.C.,Park,A.,Guan,J.L.体外划痕试验:一种分析体外细胞迁移的方便且廉价的方法。国家协议2,329-333。DOI:10.1038/nprot.2007.30。
[30] 利普什塔特,A。;Loinger,A。;新泽西州巴拉班。;Biham,O.,《无合作绑定的遗传拨动开关》,Phys。修订版Lett。,96,第188101条pp.(2006)
[31] F.R.麦克法兰。;牧师,M.A.J。;Lorenzi,T.,探索空间相互作用和抗原识别在实体肿瘤免疫反应中的作用的随机个体模型,J.Theor。生物学,480,43-55(2019)·2014年9月14日
[32] 麦克林,P。;埃杰顿,M.E。;汤普森,A.M。;Cristini,V.,《导管原位癌(DCIS)的基于患者校准试剂的建模:从微观测量到临床进展的宏观预测》,J.Theor。生物学,301122-140(2012)·Zbl 1397.92346号
[33] O.J.麦克拉伦。;Nicholson,R。;英国比亚卡森。;奥沙利文,J.P。;O'Sullivan,M.L.,《将后验信息近似误差纳入分层框架,以促进地热反演问题和不确定性量化的现成MCMC采样》,《水资源》。第56号决议,第24396条pp.(2020年)
[34] Macnamara,C.K。;E.I.米切尔。;Chaplain,M.A.J.,合成基因调控网络的空间随机建模,J.Theor。生物学,468,27-44(2019)·Zbl 1411.92121号
[35] Maini,P.K。;McElwain,D.L.S。;Leavesley,D.I.,《解释人腹膜间皮细胞创伤愈合细胞迁移分析的行波模型》,《组织工程》,第10期,第475-482页(2004年)
[36] Maini,P.K。;McElwain,D.L.S。;Leavesley,D.,《伤口愈合试验中的行波》,应用。数学。莱特。,17, 575-580 (2004) ·Zbl 1055.92025号
[37] McKay,M.D.,《计算机代码输出分析中选择输入变量值的三种方法的比较》,《技术计量学》,21,2,239-245(1979)·Zbl 0415.62011号
[38] 莫特·R·L。;罗斯,R.J.H。;Hainey,K.J。;哈里森,O.J。;Keighren,文学硕士。;兰迪尼,G。;贝克·R·E。;Painter,K.J。;I.J.Jackson。;Yates,C.A.,《用基本数学模型调和不同哺乳动物的色素沉着模式》,Nat.Commun。,7, 1-13 (2016)
[39] Murray,J.D.,《数学生物学I:导论》(2002),施普林格出版社:纽约施普林格出版社·Zbl 1006.92001号
[40] Myers,R.H。;蒙哥马利特区。;Anderson-Cook,C.M.,《响应面方法:使用设计实验优化工艺和产品》(2016),John Wiley&Sons·Zbl 1332.62004号
[41] Owen,A.B.,《拉丁超立方体采样的中心极限定理》,J.R.Stat.Soc.:Ser。B(方法学)。,54, 2, 541-551 (1992) ·兹比尔0776.62041
[42] Pawitan,Y.,《所有可能性:使用可能性的统计建模和推断》(2013),牛津大学出版社:牛津大学出版社·Zbl 1256.62006年
[43] Powell,M.J.D.,《无导数边界约束优化的BOBYQA算法》(2009),应用数学和理论物理系:英国剑桥应用数学与理论物理系,(技术报告)
[44] Raue,A。;克鲁茨,C。;Maiwald,T。;巴赫曼,J。;席林,M。;Klingmüller,美国。;Timmer,J.,《利用剖面可能性对部分观测到的动力学模型进行结构和实际可识别性分析》,生物信息学,25(2009),123-12
[45] Royston,P.,估计和置信区间的剖面似然,Stata J.,7,376-387(2007)
[46] 桑特纳,T.J。;威廉姆斯,B.J。;诺茨,W.I。;Williams,B.J.,《计算机实验的设计与分析》(2003),Springer:Springer纽约·Zbl 1041.62068号
[47] 辛普森,M.J。;兰德曼,K.A。;休斯,B.D。;Fernando,A.E.,流动人口中尺度模式的模型,《物理学A》,389,1412-1424(2010)
[48] 辛普森,M.J。;兰德曼,K.A。;Hughes,B.D.,《基于延时数据的细胞浸润与增殖机制》,Physica A,389,3779-3790(2010)
[49] 辛普森,M.J。;粘合剂,B.J。;哈里达斯,P。;伍德,B.K。;Treloar,K.K。;McElwain,D.L.S。;Baker,R.E.,聚类单层发育的实验和建模研究,Bull。数学。生物学,75,871-889(2013)·Zbl 1273.92020年
[50] 辛普森,M.J。;夏普,J.A。;Baker,R.E.,《区分平均场、力矩动力学和出生死亡运动过程的随机描述》,《物理学A》,395,236-246(2014)·Zbl 1395.60092号
[51] 辛普森,M.J。;布朗宁,A.P。;沃恩,D.J。;O.J.麦克拉伦。;Baker,R.E.,sigmoid人口增长模型的参数可识别性和模型选择,J.Theor。生物学,535,第110998条pp.(2022)·Zbl 1483.92113号
[52] 特雷罗,K.K。;辛普森,M.J。;McElwain,D.L.S。;Baker,R.E.,细胞扩散率和细胞增殖率的体外评估是否对分析几何形状敏感?,J.西奥。生物,35671-84(2014)·Zbl 1412.92027号
[53] Tremel,A。;蔡,A。;Tirtaatmadja,N。;休斯,D.B。;史蒂文斯,G.W。;兰德曼,K.A。;O'Connor,A.J.,单层形成和伤口愈合过程中的细胞迁移和增殖,《化学》。工程科学。,64, 247-253 (2009)
[54] Vellela,M。;Qian,H.,随机化学反应的准平稳分析:Keizer佯谬,Bull。数学。《生物学》,68,1727-1746(2007)·Zbl 1298.92129号
[55] Villaverde,A.F。;Pathirana,D。;Fröhlich,F。;Hasenauer,J。;Banga,J.R.,《动态模型校准协议》,Briefings Bioinf。,23, 1-19 (2022)
[56] 维特代罗,S.T。;McCue,S.W。;古纳辛格,G。;新泽西州哈斯。;Simpson,M.J.,《细胞迁移的实时细胞周期动力学数学模型》,生物物理。J.,1141241-1253(2018)
[57] Vo、B.N。;Drovandi,C.C。;Pettitt,A.N。;Simpson,M.J.,使用近似贝叶斯计算量化集体细胞扩散随机模型参数估计的不确定性,数学。生物科学。,263, 133-142 (2005) ·Zbl 1371.92033号
[58] Wooldridge,J.M.,《横截面和面板数据的计量经济学分析》(2010年),麻省理工学院出版社·兹比尔1327.62009
[59] 徐,T。;瓦洛奇,A.J。;Ye,M。;Liang,F.,量化模型结构误差:使用代理和数据驱动误差模型对区域地下水流量模型进行有效的贝叶斯校准,Water Resour。第3号决议,4084-4105(2017年)
[60] 张,D。;布林纳斯,I.M。;粘合剂,B.J。;兰德曼,K.A。;Newgreen,D.F.,《肠神经系统形成的神经嵴区域化:先天性巨结肠和干细胞治疗的意义》,发展。《生物学》,339280-294(2010)
[61] Ziaeipoor,H。;泰勒,M。;Pandy,M。;Martelli,S.,实时预测股骨应变的新型无训练方法,J.Biomech。,86, 110-116 (2019)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。