{"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}
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.