相互作用区室中反应网络模型的参数敏感性分析

David F. Anderson, Aidan S. Howells
{"title":"相互作用区室中反应网络模型的参数敏感性分析","authors":"David F. Anderson, Aidan S. Howells","doi":"arxiv-2408.09208","DOIUrl":null,"url":null,"abstract":"Models of reaction networks within interacting compartments (RNIC) are a\ngeneralization of stochastic reaction networks. It is most natural to think of\nthe interacting compartments as ``cells'' that can appear, degrade, split, and\neven merge, with each cell containing an evolving copy of the underlying\nstochastic reaction network. Such models have a number of parameters, including\nthose associated with the internal chemical model and those associated with the\ncompartment interactions, and it is natural to want efficient computational\nmethods for the numerical estimation of sensitivities of model statistics with\nrespect to these parameters. Motivated by the extensive work on computational\nmethods for parametric sensitivity analysis in the context of stochastic\nreaction networks over the past few decades, we provide a number of methods in\nthe RNIC setting. Provided methods include the (unbiased) Girsanov\ntransformation method (also called the Likelihood Ratio method) and a number of\ncoupling methods for the implementation of finite differences. We provide\nseveral numerical examples and conclude that the method associated with the\n``Split Coupling'' provides the most efficient algorithm. This finding is in\nline with the conclusions from the work related to sensitivity analysis of\nstandard stochastic reaction networks. We have made all of the Matlab code used\nto implement the various methods freely available for download.","PeriodicalId":501266,"journal":{"name":"arXiv - QuanBio - Quantitative Methods","volume":null,"pages":null},"PeriodicalIF":0.0000,"publicationDate":"2024-08-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Parametric Sensitivity Analysis for Models of Reaction Networks within Interacting Compartments\",\"authors\":\"David F. Anderson, Aidan S. Howells\",\"doi\":\"arxiv-2408.09208\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"Models of reaction networks within interacting compartments (RNIC) are a\\ngeneralization of stochastic reaction networks. It is most natural to think of\\nthe interacting compartments as ``cells'' that can appear, degrade, split, and\\neven merge, with each cell containing an evolving copy of the underlying\\nstochastic reaction network. Such models have a number of parameters, including\\nthose associated with the internal chemical model and those associated with the\\ncompartment interactions, and it is natural to want efficient computational\\nmethods for the numerical estimation of sensitivities of model statistics with\\nrespect to these parameters. Motivated by the extensive work on computational\\nmethods for parametric sensitivity analysis in the context of stochastic\\nreaction networks over the past few decades, we provide a number of methods in\\nthe RNIC setting. Provided methods include the (unbiased) Girsanov\\ntransformation method (also called the Likelihood Ratio method) and a number of\\ncoupling methods for the implementation of finite differences. We provide\\nseveral numerical examples and conclude that the method associated with the\\n``Split Coupling'' provides the most efficient algorithm. This finding is in\\nline with the conclusions from the work related to sensitivity analysis of\\nstandard stochastic reaction networks. We have made all of the Matlab code used\\nto implement the various methods freely available for download.\",\"PeriodicalId\":501266,\"journal\":{\"name\":\"arXiv - QuanBio - Quantitative Methods\",\"volume\":null,\"pages\":null},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2024-08-17\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"arXiv - QuanBio - Quantitative Methods\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/arxiv-2408.09208\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - QuanBio - Quantitative Methods","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2408.09208","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

摘要

相互作用区内的反应网络(RNIC)模型是随机反应网络的一般化。最自然的做法是将相互作用的隔室视为 "细胞",它们可以出现、退化、分裂甚至合并,每个细胞都包含一个不断演化的基本随机反应网络的副本。这种模型有许多参数,包括与内部化学模型相关的参数和与区室相互作用相关的参数,因此自然需要高效的计算方法来数值估计模型统计量对这些参数的敏感性。过去几十年来,在随机反应网络背景下,针对参数敏感性分析的计算方法开展了大量工作,受此激励,我们提供了一些 RNIC 环境下的方法。所提供的方法包括(无偏)吉尔萨诺夫变换方法(也称似然比方法)和一些用于实现有限差分的耦合方法。我们提供了几个数值示例,并得出结论:与 "拆分耦合 "相关的方法提供了最有效的算法。这一结论与标准随机反应网络灵敏度分析相关工作的结论一致。我们已将用于实现各种方法的所有 Matlab 代码免费提供下载。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
Parametric Sensitivity Analysis for Models of Reaction Networks within Interacting Compartments
Models of reaction networks within interacting compartments (RNIC) are a generalization of stochastic reaction networks. It is most natural to think of the interacting compartments as ``cells'' that can appear, degrade, split, and even merge, with each cell containing an evolving copy of the underlying stochastic reaction network. Such models have a number of parameters, including those associated with the internal chemical model and those associated with the compartment interactions, and it is natural to want efficient computational methods for the numerical estimation of sensitivities of model statistics with respect to these parameters. Motivated by the extensive work on computational methods for parametric sensitivity analysis in the context of stochastic reaction networks over the past few decades, we provide a number of methods in the RNIC setting. Provided methods include the (unbiased) Girsanov transformation method (also called the Likelihood Ratio method) and a number of coupling methods for the implementation of finite differences. We provide several numerical examples and conclude that the method associated with the ``Split Coupling'' provides the most efficient algorithm. This finding is in line with the conclusions from the work related to sensitivity analysis of standard stochastic reaction networks. We have made all of the Matlab code used to implement the various methods freely available for download.
求助全文
通过发布文献求助,成功后即可免费获取论文全文。 去求助
来源期刊
自引率
0.00%
发文量
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信