[发明专利]一种多重熔解曲线的软件分析方法

权利要求书

1.一种多重熔解曲线的软件分析方法，其特征在于，所述方法包括：

S1，获取或接收包括多条多重熔解曲线的数据集；

S2，对所述数据集中的每条多重熔解曲线进行校正，获得归一化的多重熔解曲线；

S3，根据归一化的多重熔解曲线，计算得到目标峰；

S4，对归一化的目标峰进行聚类分析，确定目标基因的类别。

2.根据权利要求1所述的软件分析方法，其特征在于，每条所述多重熔解曲线包括熔解前区域、熔解区及熔解后区域，所述熔解区包括至少一个熔解段，且所述多个熔解段中，其中一个为参考基因的熔解段，其它的均为目标基因的熔解段。

3.根据权利要求1所述的软件分析方法，其特征在于，所述S2包括：

S21，统计所述数据集中的最大降幅，按最大降幅比例的经验值设定阴阳性阈值，并根据所述阴阳性阈值去除所述数据集中的阴性曲线；

S22，定位去除阴性曲线后的数据集的熔解前区域和熔解后区域；

S23，根据所述数据集的熔解前区域与熔解后区域，去除每条多重熔解曲线的背景曲线，获得无背景曲线的第一多重熔解曲线；

S24，对所述第一多重熔解曲线进行归一化，获得归一化后的第二多重熔解曲线。

4.根据权利要求3所述的软件分析方法，其特征在于，所述S22包括：

S221，对每条多重熔解曲线进行开窗，得到多个窗区，依次计算每个窗区的线性相关系数，并根据所述线性相关系数判断当前窗区是否为目标区域，其中，每个窗区的开窗长度为n个数据点，n为大于等于1的整数；

S222，统计去除阴性曲线后的数据集内所有多重熔解曲线的熔解前区域的最小起始温度，根据所述最小起始温度确定数据集的熔解前区域；及统计S1数据集内所有多重熔解曲线的熔解后区域的最大结束温度，根据所述最大结束温度确定数据集的熔解后区域。

5.根据权利要求4所述的软件分析方法，其特征在于，所述S221中，n的范围由多重熔解曲线的x轴的温度间距决定，其计算公式为：

其中，step为温度间距。

6.根据权利要求4所述的软件分析方法，其特征在于，所述S221中，判断当前窗区是否为目标区域的过程包括：若当前窗区的线性相关系数比前两个窗区的线性相关系数大，且比后一个窗区的线性相关系数大，则判定当前窗区为目标区域。

7.根据权利要求3所述的软件分析方法，其特征在于，所述S23包括：

S231，根据数据集的熔解前区域与熔解后区域，计算每条多重熔解曲线的背景曲线，计算公式如下：

y＝c×e^a(x-T0)，

其中，y表示计算所得的背景曲线；

k0为熔解前区域的拟合直线的斜率，

k1为熔解后区域的拟合直线的斜率，

T0为熔解前区域的平均温度，

T1为熔解后区域的平均温度；

S232，用多重熔解曲线减去背景曲线，得到无背景曲线的所述第一多重熔解曲线。

8.根据权利要求3所述的软件分析方法，其特征在于，所述S24中，归一化的公式如下：

其中，F2为归一化后的第二多重熔解曲线，F1为第一多重熔解曲线，Max，Min分别为F1的最大值与最小值。

9.根据权利要求1所述的软件分析方法，其特征在于，所述S3包括：

S31，对每条归一化的第二多重熔解曲线进行负一阶导的计算，得到峰值曲线；

S32，对每条所述峰值曲线进行开窗，且在峰值曲线上搜索目标峰；

S33，根据参考基因的熔解温度，确定参考峰；

S34，将峰值曲线归一化，得到归一化峰值曲线上的目标熔解峰。

10.根据权利要求1所述的软件分析方法，其特征在于，所述S4包括：

S41，根据标准品的类别，确定目标熔解峰的聚类分析的类别；

S42，将所有目标熔解峰根据峰高进行聚类，并将每个目标熔解峰类别映射到标准品的相应类别中；

S43，统计每个目标熔解峰类别的均值与方差，并根据所述均值与方差估计当前目标熔解峰类别在其映射的标准品相应类别中的概率；

S44，根据所述概率及预设的概率阈值，确定目标熔解峰分类正确与否。