【数学/计算问题求教】 (转载)# 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里面没有对数增长过程。
大家有可行的计算/理论方法吗?
发信人: 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里面没有对数增长过程。
大家有可行的计算/理论方法吗?