CN109637586A - 测序深度的矫正方法及装置 - Google Patents
测序深度的矫正方法及装置 Download PDFInfo
- Publication number
- CN109637586A CN109637586A CN201811616187.2A CN201811616187A CN109637586A CN 109637586 A CN109637586 A CN 109637586A CN 201811616187 A CN201811616187 A CN 201811616187A CN 109637586 A CN109637586 A CN 109637586A
- Authority
- CN
- China
- Prior art keywords
- section
- site
- depth
- content
- sites
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明提供了一种测序深度的矫正方法及装置,其中,上述装置包括:获取模块,用于获取基因组上每个位点的GC含量,以及所述基因组的所有位点的测序深度;划分模块,用于根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,其中,N为正整数;确定模块,用于根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布;处理模块,用于对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度。
Description
技术领域
本发明涉及通信领域,具体而言,涉及一种测序深度的矫正方法及装置。
背景技术
一些研究表明,通过扩增和测序引入的GC偏差对染色体非整倍性的敏感性产生了限制作用。在不同测序条件下,例如试剂成分不同、测序温度不同或者簇密度不同,GC偏差可能在样品制备和测序过程中引入,这会导致GC含量过高或GC含量过低的基因组区域的测序数据出现偏差。
相关技术中开发了消除GC偏差的计算方法。HC Fan和SR Quake开发了一种通过局部多项式回归拟合(LOESS回归)的计算方法消除GC偏差,该方法通过基于局部基因组GC含量对每个GC密度给予权重,以通过乘以相应权重矫正每个区域的Reads数目。但是,该方法在计算时可能造成数据轻微失真,会干扰检测的准确性。
针对相关技术中,通过局部多项式回归拟合的计算方法来消除测序深度的偏差的方法,容易出现数据失真,结果准确性差等问题,尚未提出有效的技术方案。
发明内容
本发明实施例提供了一种测序深度的矫正方法及装置,以至少解决相关技术中通过局部多项式回归拟合的计算方法来消除测序深度的偏差的方法,容易出现数据失真,结果准确性差等问题。
根据本发明的一个实施例,提供了一种测序深度的矫正装置,包括:获取模块,用于获取基因组上每个位点的GC含量,以及所述基因组的所有位点的测序深度;划分模块,用于根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,其中,N为正整数;确定模块,用于根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布;处理模块,用于对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度。
在本发明实施例中,所述划分模块,用于将GC含量落在第一范围内的位点化分在N个不同区间的两端区间内;以及将GC含量落在第二范围内的位点化分在中间区间内,其中,中间区间为所述N个不同区间中,除所述两端区间的其他区间,所述第一范围内的位点数目小于所述第二范围内的位点数目。
在本发明实施例中,所述确定模块,还用于根据以下公式确定所述泊松分布:p(x=n|λ)=eλλn/n!,其中,λ=λall-λi,λi为所述区间平均深度,n的取值包括:0、1、2、3、4、5,λall为所述N个区间的平均深度。
在本发明实施例中,所述确定模块,还用于通过以下方式确定所述区间平均深度:将区间内所有位点的总深度除以该区间内位点的个数,得到所述区间平均深度。
在本发明实施例中,所述确定模块,还用于通过以下方式确定所述N个区间的平均深度:将所述基因组内所有位点的总深度除以所述基因组内位点的个数,得到所述N个区间的平均深度。
在本发明实施例中,所述处理模块,还用于将矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值的和,作为矫正后的测序深度。
在本发明实施例中,所述装置还包括:矫正模块,用于根据矫正后的每个位点的测序深度和所述每个位点的GC含量矫正所述每个位点的GC偏好。
根据本发明的另一个实施例,还提供了一种测序深度的矫正方法,包括:获取基因组上每个位点的GC含量,以及所述基因组的所有位点的测序深度;根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,其中,N为正整数;根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布;对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度。
在本发明实施例中,根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,包括:将GC含量落在第一范围内的位点化分在N个不同区间的两端区间内;将GC含量落在第二范围内的位点化分在中间区间内,其中,中间区间为所述N个不同区间中,除所述两端区间的其他区间,所述第一范围内的位点数目小于所述第二范围内的位点数目。.
在本发明实施例中,根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布,包括:p(x=n|λ)=eλλn/n!,其中,λ=λall-λi,λi为所述区间平均深度,n的取值包括:0、1、2、3、4、5,λall为所述N个区间的平均深度。
根据本发明的另一个实施例,还提供了一种存储介质,所述存储介质中存储有计算机程序,其中,所述计算机程序被设置为运行时执行以上任一项所述的测序深度的矫正方法。
通过本发明,根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布,对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度,采用上述技术方案,解决了相关技术中通过局部多项式回归拟合的计算方法来消除测序深度的偏差的方法,容易出现数据失真,结果准确性差等问题,进而提高了一种准确性高的方法来消除测序深度的偏差。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为根据本发明实施例的测序深度的矫正方法的流程图;
图2是根据本发明实施例的测序深度的矫正装置的结构框图;
图3是根据本发明实施例的测序深度的矫正装置的另一结构框图;
图4是根据本发明实施例的测序深度的矫正前效果的示意图(一);
图5是根据本发明实施例的测序深度的矫正后效果的示意图(一);
图6是根据本发明实施例的测序深度的矫正前效果的示意图(二);
图7是根据本发明实施例的测序深度的矫正后效果的示意图(二)。
具体实施方式
下文中将参考附图并结合实施例来详细说明本发明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。
相关技术中,在对测序深度进行矫正的过程中,当数据量很大时,要对每个待预测数据拟合一次,运行时间较长,计算量大,本发明实施例将泊松分布引入到了测序深度的矫正过程中,避免了欠拟合和过拟合引起的计算过程,具体详见以下实施例以及优选实施例的技术方案。
实施例1
本发明实施例提供了一种测序深度的矫正方法,图1为根据本发明实施例的测序深度的矫正方法的流程图,如图1所示,包括以下步骤:
步骤S102:获取基因组上每个位点的GC含量,以及所述基因组的所有位点的测序深度;
步骤S104:根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,其中,N为正整数;
步骤S106:根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布;
步骤S108:对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度。
通过上述步骤,根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布,对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度,采用上述技术方案,解决了相关技术中通过局部多项式回归拟合的计算方法来消除测序深度的偏差的方法,容易出现数据失真,结果准确性差等问题,进而提高了一种准确性高的方法来消除测序深度的偏差,且采用泊松分布的计算方式,运算速度快。
上述步骤S104可以有多种实现方式,在一个可选实施例中,可以通过以下技术方案实现:将GC含量落在第一范围内的位点化分在N个不同区间的两端区间内;将GC含量落在第二范围内的位点化分在中间区间内,其中,中间区间为所述N个不同区间中,除所述两端区间的其他区间,所述第一范围内的位点数目小于所述第二范围内的位点数目,本发明实施例中的第二范围应该指的是在含有位点数目比较多的GC区间,而第一范围应该指的是含有位点数目比较少的GC区间,对于基因组,通常GC含量在20%~60%的位点比较多,也就是说GC在20%~60%的位点为第二范围的位点,其余的位点位第一范围的位点。
在本发明实施例中,确定每个区间的区间平均深度,具体可以通过该区间的所有位点总深度除以总位点数来得到;确定N个区间的平均深度,具体可以基因组的所有位点总深度除以总位点数来得到。
在本发明实施例中,根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布,包括:p(x=n|λ)=eλλn/n!,其中,λ=λall-λi,λi为所述区间平均深度,n的取值包括:0、1、2、3、4、5,λall为所述N个区间的平均深度。
可选地,在步骤S108后,还可以执行以下技术方案:根据矫正后的每个位点的测序深度和所述每个位点的GC含量矫正所述每个位点的GC偏好,也就是说,在本发明实施例中,步骤S108得到的测序深度可以理解为一种中间状态的信息,其中,测序深度与GC含量的相关性,称为GC偏好(英文简称GC bias),GC含量是在DNA4种碱基中,鸟嘌呤和胞嘧啶所占的比率称为GC含量。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到根据上述实施例的方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中,包括若干指令用以使得一台终端设备(可以是手机,计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
实施例2
在本实施例中还提供了一种测序深度的矫正装置,该装置用于实现上述实施例及优选实施方式,已经进行过说明的不再赘述。如以下所使用的,术语“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图2是根据本发明实施例的测序深度的矫正装置的结构框图,如图2所示,该装置包括:
获取模块20,用于获取基因组上每个位点的GC含量,以及所述基因组的所有位点的测序深度;
划分模块22,用于根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,其中,N为正整数;
确定模块24,用于根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布;
处理模块26,用于对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度。
通过本发明,根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布,对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度,采用上述技术方案,解决了相关技术中通过局部多项式回归拟合的计算方法来消除测序深度的偏差的方法,容易出现数据失真,结果准确性差等问题,进而提高了一种准确性高的方法来消除测序深度的偏差,且采用泊松分布的计算方式,运算速度快。
在本发明实施例中,所述划分模块22,用于将GC含量落在第一范围内的位点化分在N个不同区间的两端区间内;以及将GC含量落在第二范围内的位点化分在中间区间内,其中,中间区间为所述N个不同区间中,除所述两端区间的其他区间,所述第一范围内的位点数目小于所述第二范围内的位点数目,本发明实施例中的第二范围应该指的是在含有位点数目比较多的GC区间,而第一范围应该指的是含有位点数目比较少的GC区间,对于基因组,通常GC含量在20%~60%的位点比较多,也就是说GC在20%~60%的位点为第二范围的位点,其余的位点位第一范围的位点。
在本发明实施例中,所述确定模块24,还用于根据以下公式确定所述泊松分布:p(x=n|λ)=eλλn/n!,其中,λ=λall-λi,λi为所述区间平均深度,n的取值包括:0、1、2、3、4、5,λall为所述N个区间的平均深度。
图3是根据本发明实施例的测序深度的矫正装置的另一结构框图,如图3所示,在本发明实施例中,所述装置,还包括:矫正模块28,用于根据矫正后的每个位点的测序深度和所述每个位点的GC含量矫正所述每个位点的GC偏好。
需要说明的是,上述实施例1-实施例2的技术方案可以结合使用,也可以单独使用,本发明实施例对此不作限定。
以下结合一优选实施例对上述测序深度的矫正流程进行说明,但不用于限定本发明实施例的技术方案,主要包括以下步骤:
步骤1、通过高通量测序获得样本的全基因组序列。
步骤2、将测序得到的序列与基因组hg19进行比对,比对结果去除比对到重复序列和比对到多个位置上的序列,去除非完全比对的序列。
步骤3、计算基因组hg19上每个位点的GC含量,将所有位点根据其GC含量分成43个Bin,中间Bin的区间为1%的GC含量,两端Bin的区间为20%的GC含量,其中,将基因组根据GC含量分成43个Bin,每个Bin上所有位点的深度服从泊松分布。
步骤4、利用Samtools软件(也可以是其他可执行软件)对基因组上的所有位点进行深度分析,统计所有位点的测序深度信息。
步骤5、对于每一个GC Bin,所有位点的测序深度近乎服从泊松分布。分别计算每一个Bin的平均深度λi(i为1,2...n,代表第1,2...n个Bin)和所有Bin的平均深度λall。
步骤6、对于每一个Bin i,以λ=λall-λi模拟一个新的泊松分布,将该Bin里面所有位点校正前的深度加上新的泊松分布各位点的值即为矫正后所有位点的测序深度。
步骤7、以100kb为一个滑动窗口Bin,统计基因组上不同窗口Bin的平均深度信息。
步骤8、计算每个Bin的GC含量,以GC含量为横轴,平均深度为纵轴,绘制同一样品矫正前后的GC分布图,比较矫正前后的深度信息是否有GC偏好。
基于上述步骤1-步骤8的技术方案,如图4所示,对平均深度随GC含量偏上的样品进行校正,校正前后的结果如图4-5所示,横轴是每一个Bin的GC含量,纵轴是每一个Bin的平均深度。
如图6-7所示,对平均深度随GC含量偏下的样品进行校正,校正前后的结果如图6-7所示,横轴是每一个Bin的GC含量,纵轴是每一个Bin的平均深度。
可选地,在本实施例中,上述存储介质可以包括但不限于:U盘、只读存储器(Read-Only Memory,简称为ROM)、随机存取存储器(Random Access Memory,简称为RAM)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
可选地,本实施例中的具体示例可以参考上述实施例及可选实施方式中所描述的示例,本实施例在此不再赘述。
显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种测序深度的矫正装置,其特征在于,包括:
获取模块,用于获取基因组上每个位点的GC含量,以及所述基因组的所有位点的测序深度;
划分模块,用于根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,其中,N为正整数;
确定模块,用于根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布;
处理模块,用于对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度。
2.根据权利要求1所述的装置,其特征在于,所述划分模块,用于将GC含量落在第一范围内的位点化分在N个不同区间的两端区间内;以及将GC含量落在第二范围内的位点化分在中间区间内,其中,中间区间为所述N个不同区间中,除所述两端区间的其他区间,所述第一范围内的位点数目小于所述第二范围内的位点数目。
3.根据权利要求1所述的装置,其特征在于,所述确定模块,还用于根据以下公式确定所述泊松分布:
p(x=n|λ)=eλλn/n!,其中,λ=λall-λi,λi为第i个区间的区间平均深度,n的取值包括:0、1、2、3、4、5,λall为所述N个区间的平均深度。
4.根据权利要求3所述的装置,其特征在于,所述确定模块,还用于通过以下方式确定所述区间平均深度:将区间内所有位点的总深度除以该区间内位点的个数,得到所述区间平均深度。
5.根据权利要求3所述的装置,其特征在于,所述确定模块,还用于通过以下方式确定所述N个区间的平均深度:将所述基因组内所有位点的总深度除以所述基因组内位点的个数,得到所述N个区间的平均深度。
6.根据权利要求1所述的装置,其特征在于,所述处理模块,还用于将矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值的和,作为矫正后的测序深度。
7.根据权利要求1所述的装置,其特征在于,所述装置还包括:
矫正模块,用于根据矫正后的每个位点的测序深度和所述每个位点的GC含量矫正所述每个位点的GC偏好。
8.一种测序深度的矫正方法,其特征在于,包括:
获取基因组上每个位点的GC含量,以及所述基因组的所有位点的测序深度;
根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,其中,N为正整数;
根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布;
对于所述每个区间,均根据矫正前的每个区间内位点的区间平均深度和所述泊松分布中各位点的值得到矫正后的测序深度。
9.根据权利要求8所述的方法,其特征在于,根据每个位点的GC含量将所述所有位点按照指定规则分成N个不同区间,包括:
将GC含量落在第一范围内的位点化分在N个不同区间的两端区间内;
将GC含量落在第二范围内的位点化分在中间区间内,其中,中间区间为所述N个不同区间中,除所述两端区间的其他区间,所述第一范围内的位点数目小于所述第二范围内的位点数目。
10.根据权利要求8所述的方法,其特征在于,根据每个区间内位点的区间平均深度和N个区间的平均测序深度确定泊松分布,包括:
p(x=n|λ)=eλλn/n!,其中,λ=λall-λi,λi为所述区间平均深度,n的取值包括:0、1、2、3、4、5,λall为所述N个区间的平均深度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811616187.2A CN109637586B (zh) | 2018-12-27 | 2018-12-27 | 测序深度的矫正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811616187.2A CN109637586B (zh) | 2018-12-27 | 2018-12-27 | 测序深度的矫正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109637586A true CN109637586A (zh) | 2019-04-16 |
CN109637586B CN109637586B (zh) | 2020-11-17 |
Family
ID=66078488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811616187.2A Active CN109637586B (zh) | 2018-12-27 | 2018-12-27 | 测序深度的矫正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109637586B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2009293354A1 (en) * | 2008-09-20 | 2010-03-25 | The Board Of Trustees Of The Leland Stanford Junior University | Noninvasive diagnosis of fetal aneuploidy by sequencing |
US20100216153A1 (en) * | 2004-02-27 | 2010-08-26 | Helicos Biosciences Corporation | Methods for detecting fetal nucleic acids and diagnosing fetal abnormalities |
CN103201744A (zh) * | 2010-10-13 | 2013-07-10 | 考利达基因组股份有限公司 | 用于估算全基因组拷贝数变异的方法 |
CN104428425A (zh) * | 2012-05-04 | 2015-03-18 | 考利达基因组股份有限公司 | 测定复杂肿瘤全基因组绝对拷贝数变异的方法 |
CN107480470A (zh) * | 2016-06-08 | 2017-12-15 | 广州华大基因医学检验所有限公司 | 基于贝叶斯与泊松分布检验的已知变异检出方法和装置 |
CN108229101A (zh) * | 2017-12-29 | 2018-06-29 | 北京科迅生物技术有限公司 | 基于ngs的靶向测序数据模拟方法和装置 |
-
2018
- 2018-12-27 CN CN201811616187.2A patent/CN109637586B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100216153A1 (en) * | 2004-02-27 | 2010-08-26 | Helicos Biosciences Corporation | Methods for detecting fetal nucleic acids and diagnosing fetal abnormalities |
AU2009293354A1 (en) * | 2008-09-20 | 2010-03-25 | The Board Of Trustees Of The Leland Stanford Junior University | Noninvasive diagnosis of fetal aneuploidy by sequencing |
CN103201744A (zh) * | 2010-10-13 | 2013-07-10 | 考利达基因组股份有限公司 | 用于估算全基因组拷贝数变异的方法 |
CN104428425A (zh) * | 2012-05-04 | 2015-03-18 | 考利达基因组股份有限公司 | 测定复杂肿瘤全基因组绝对拷贝数变异的方法 |
CN107480470A (zh) * | 2016-06-08 | 2017-12-15 | 广州华大基因医学检验所有限公司 | 基于贝叶斯与泊松分布检验的已知变异检出方法和装置 |
CN108229101A (zh) * | 2017-12-29 | 2018-06-29 | 北京科迅生物技术有限公司 | 基于ngs的靶向测序数据模拟方法和装置 |
Non-Patent Citations (1)
Title |
---|
张高川等: "一株出芽短梗霉全基因组序列测定及其分析", 《基因组学与应用生物学》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109637586B (zh) | 2020-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fan et al. | MuSE: accounting for tumor heterogeneity using a sample-specific error model improves sensitivity and specificity in mutation calling from sequencing data | |
McKenzie et al. | Distributions of cherries for two models of trees | |
US9792827B2 (en) | Providing questions to entity groups | |
CN110222029A (zh) | 一种大数据多维分析计算效率提升方法及系统 | |
CN108647261A (zh) | 基于气象数据离散点网格化处理的全球等值线绘制方法 | |
US10592617B2 (en) | Computer implemented method, processor device and computer program product for designing intervention into real complex systems | |
WO2012142052A1 (en) | Generation of occupant activities based on recorded occupant behavior | |
Yuan et al. | A dose distribution evaluation technique using the k‐d tree for nearest neighbor searching | |
CN109616154A (zh) | 测序深度的矫正方法及装置 | |
Salmon et al. | Asteroseismology of β Cephei stars: The stellar inferences tested in hare and hound exercises | |
CN109637586A (zh) | 测序深度的矫正方法及装置 | |
KR101635791B1 (ko) | P파 도달시간을 이용한 진원 및 진원시 결정방법 | |
CN114037531A (zh) | 一种区块链分片系统的分片账户调整方法及其相关装置 | |
CN109727638A (zh) | 测序深度的矫正方法及装置 | |
CN109637585A (zh) | 测序深度的矫正方法及装置 | |
CN109584959A (zh) | 测序深度的矫正方法及装置 | |
CN111569424B (zh) | 游戏场景内的物理结算方法和装置 | |
CN109767813A (zh) | 测序深度的矫正方法及装置 | |
CN106373192B (zh) | 一种非拓扑一致性三维网格块体追踪算法 | |
Münchmeyer | PyOcto: A high-throughput seismic phase associator | |
CN109347680B (zh) | 一种网络拓扑重构方法、装置及终端设备 | |
CN108133234A (zh) | 基于稀疏子集选择算法的社区检测方法、装置及设备 | |
JP4842494B2 (ja) | イオン分布計算方法およびプログラム | |
CN115759250A (zh) | 归因分析方法、装置、电子设备以及存储介质 | |
Caramia et al. | CHECKCOL: Improved local search for graph coloring |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |