Redian新闻
>
【数学/计算问题求教】 (转载)
avatar
【数学/计算问题求教】 (转载)# Joke - 肚皮舞运动
H*g
1
【 以下文字转载自 Military 讨论区 】
发信人: xiao86 (xiao86), 信区: Military
标 题: 【数学/计算问题求教】
发信站: BBS 未名空间站 (Thu Oct 19 14:36:54 2017, 美东)
下面这个问题,大家有好的方法简化拟合吗:
从一个细胞开始,里面5个拷贝质粒,其中一个是突变子。然后开始扩增,每一代,每
一个细胞分裂成两个,然后质粒复制一套(不考虑新突变)、之后随机分到两个子细胞
里,但是每个细胞总是保持5个拷贝的质粒。比如第一次分裂之后,有两个新细胞,8个
拷贝的没有突变的质粒,2个突变质粒。那么两个细胞可以是各取一个突变质粒,也可
以是一个细胞拿两个而另一个细胞全部都是没有突变的。
然后问题是,N代扩增以后,突变子在种群(N~10-30)里的分布:五个质粒拷贝里,有
0,1,2,3,4,5个突变拷贝的分布.
我能想到的就是拟合每一个细胞和每一次分裂;但是因为是分裂是以2为底数的指数过
程,计算量也是随N指数扩增,当N=30的时候,要计算5亿次才能得出下一代的结果,
computationally impossible。
这和urn process有点类似,但是urn里面没有对数增长过程。
大家有可行的计算/理论方法吗?
avatar
g*2
2
如果只想知道分布的话,只要对前一代取样就可以了 (比如当细胞数到了某个容忍上
线之后,只随机取上一代一半继续往下增长)。
avatar
J*e
3
这个有解析解吧?不用simulation的。

【在 H********g 的大作中提到】
: 【 以下文字转载自 Military 讨论区 】
: 发信人: xiao86 (xiao86), 信区: Military
: 标 题: 【数学/计算问题求教】
: 发信站: BBS 未名空间站 (Thu Oct 19 14:36:54 2017, 美东)
: 下面这个问题,大家有好的方法简化拟合吗:
: 从一个细胞开始,里面5个拷贝质粒,其中一个是突变子。然后开始扩增,每一代,每
: 一个细胞分裂成两个,然后质粒复制一套(不考虑新突变)、之后随机分到两个子细胞
: 里,但是每个细胞总是保持5个拷贝的质粒。比如第一次分裂之后,有两个新细胞,8个
: 拷贝的没有突变的质粒,2个突变质粒。那么两个细胞可以是各取一个突变质粒,也可
: 以是一个细胞拿两个而另一个细胞全部都是没有突变的。

avatar
j*g
4
vector (E1, E2, E3, E4, E5)
transition matrix (E1E1, E1E2, E1E3, ..., E2E1, E2E2, E2E3, ...)
v0 = (1, 0, 0, 0, 0)
相关阅读
logo
联系我们隐私协议©2024 redian.news
Redian新闻
Redian.news刊载任何文章,不代表同意其说法或描述,仅为提供更多信息,也不构成任何建议。文章信息的合法性及真实性由其作者负责,与Redian.news及其运营公司无关。欢迎投稿,如发现稿件侵权,或作者不愿在本网发表文章,请版权拥有者通知本网处理。