[發(fā)明專利]一種多重熔解曲線的軟件分析方法

在審

• 申請?zhí)枺?/td>: 201911371056.7
• 申請日: 2019-12-26
• 公開（公告）號: CN111145835A
• 公開（公告）日: 2020-05-12
• 發(fā)明（設計）人: 苗保剛;李政;羅淑芬;孫瑤;李明;彭年才
• 申請（專利權）人: 蘇州天隆生物科技有限公司
• 主分類號: G16B50/00
• 分類號: G16B50/00
• 代理公司: 蘇州集律知識產(chǎn)權代理事務所(普通合伙) 32269
• 代理人: 安紀平

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

S1，獲取或接收包括多條多重熔解曲線的數(shù)據(jù)集；

S2，對所述數(shù)據(jù)集中的每條多重熔解曲線進行校正，獲得歸一化的多重熔解曲線；

S3，根據(jù)歸一化的多重熔解曲線，計算得到目標峰；

S4，對歸一化的目標峰進行聚類分析，確定目標基因的類別。

2.根據(jù)權利要求1所述的軟件分析方法，其特征在于，每條所述多重熔解曲線包括熔解前區(qū)域、熔解區(qū)及熔解后區(qū)域，所述熔解區(qū)包括至少一個熔解段，且所述多個熔解段中，其中一個為參考基因的熔解段，其它的均為目標基因的熔解段。

3.根據(jù)權利要求1所述的軟件分析方法，其特征在于，所述S2包括：

S21，統(tǒng)計所述數(shù)據(jù)集中的最大降幅，按最大降幅比例的經(jīng)驗值設定陰陽性閾值，并根據(jù)所述陰陽性閾值去除所述數(shù)據(jù)集中的陰性曲線；

S22，定位去除陰性曲線后的數(shù)據(jù)集的熔解前區(qū)域和熔解后區(qū)域；

S23，根據(jù)所述數(shù)據(jù)集的熔解前區(qū)域與熔解后區(qū)域，去除每條多重熔解曲線的背景曲線，獲得無背景曲線的第一多重熔解曲線；

S24，對所述第一多重熔解曲線進行歸一化，獲得歸一化后的第二多重熔解曲線。

4.根據(jù)權利要求3所述的軟件分析方法，其特征在于，所述S22包括：

S221，對每條多重熔解曲線進行開窗，得到多個窗區(qū)，依次計算每個窗區(qū)的線性相關系數(shù)，并根據(jù)所述線性相關系數(shù)判斷當前窗區(qū)是否為目標區(qū)域，其中，每個窗區(qū)的開窗長度為n個數(shù)據(jù)點，n為大于等于1的整數(shù)；

S222，統(tǒng)計去除陰性曲線后的數(shù)據(jù)集內(nèi)所有多重熔解曲線的熔解前區(qū)域的最小起始溫度，根據(jù)所述最小起始溫度確定數(shù)據(jù)集的熔解前區(qū)域；及統(tǒng)計S1數(shù)據(jù)集內(nèi)所有多重熔解曲線的熔解后區(qū)域的最大結(jié)束溫度，根據(jù)所述最大結(jié)束溫度確定數(shù)據(jù)集的熔解后區(qū)域。

5.根據(jù)權利要求4所述的軟件分析方法，其特征在于，所述S221中，n的范圍由多重熔解曲線的x軸的溫度間距決定，其計算公式為：

其中，step為溫度間距。

6.根據(jù)權利要求4所述的軟件分析方法，其特征在于，所述S221中，判斷當前窗區(qū)是否為目標區(qū)域的過程包括：若當前窗區(qū)的線性相關系數(shù)比前兩個窗區(qū)的線性相關系數(shù)大，且比后一個窗區(qū)的線性相關系數(shù)大，則判定當前窗區(qū)為目標區(qū)域。

7.根據(jù)權利要求3所述的軟件分析方法，其特征在于，所述S23包括：

S231，根據(jù)數(shù)據(jù)集的熔解前區(qū)域與熔解后區(qū)域，計算每條多重熔解曲線的背景曲線，計算公式如下：

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

其中，y表示計算所得的背景曲線；

k0為熔解前區(qū)域的擬合直線的斜率，

k1為熔解后區(qū)域的擬合直線的斜率，

T0為熔解前區(qū)域的平均溫度，

T1為熔解后區(qū)域的平均溫度；

S232，用多重熔解曲線減去背景曲線，得到無背景曲線的所述第一多重熔解曲線。

8.根據(jù)權利要求3所述的軟件分析方法，其特征在于，所述S24中，歸一化的公式如下：

其中，F(xiàn)2為歸一化后的第二多重熔解曲線，F(xiàn)1為第一多重熔解曲線，Max，Min分別為F1的最大值與最小值。

9.根據(jù)權利要求1所述的軟件分析方法，其特征在于，所述S3包括：

S31，對每條歸一化的第二多重熔解曲線進行負一階導的計算，得到峰值曲線；

S32，對每條所述峰值曲線進行開窗，且在峰值曲線上搜索目標峰；

S33，根據(jù)參考基因的熔解溫度，確定參考峰；

S34，將峰值曲線歸一化，得到歸一化峰值曲線上的目標熔解峰。

10.根據(jù)權利要求1所述的軟件分析方法，其特征在于，所述S4包括：

S41，根據(jù)標準品的類別，確定目標熔解峰的聚類分析的類別；

S42，將所有目標熔解峰根據(jù)峰高進行聚類，并將每個目標熔解峰類別映射到標準品的相應類別中；

S43，統(tǒng)計每個目標熔解峰類別的均值與方差，并根據(jù)所述均值與方差估計當前目標熔解峰類別在其映射的標準品相應類別中的概率；

S44，根據(jù)所述概率及預設的概率閾值，確定目標熔解峰分類正確與否。