[发明专利]一种PCR检测方法的数据处理方法

权利要求书

1.一种PCR检测方法的数据处理方法，其特征在于，包括：

S1：采用光谱仪采集循环数0-40荧光强度A_i，A_i，i∈[0，40]；

S2：将采集的循环数0-40荧光强度A_i作归一化处理，得到归一化处理后的数据D_i，D_i＝C_i-y_i，i∈[1，40]；y_i为基线，C_i为平滑处理后的数据；

采集的循环数0-40荧光强度A_i归一化处理方法具体包括：

S11：扣除基底：A₁-A₄₀作为信号，减去基底A₀；

B_i＝A_i-A₀，i∈[1，40]，其中，B_i为去除基底的数据；

S12：平滑数据：对去除基底的数据B_i进行平滑处理，得到平滑处理后的数据C_i；

B_i-C_i

S13：确定基线：选择C₃-C₁₀，根据最小二乘法拟合出直线作为基线y_i，

y_i＝a₁χ_i+b₁，i∈[1，40]，其中，a₁为一次项系数，b₁为常数项；

S14：数据扣除基线，得到D_i＝C_i-y_i，i∈[1，40]；

S3：对步骤S2得到归一化处理后的数据D_i依次作一阶差分、二阶差分和三阶差分处理；

S31：一阶差分

其中，D′_i为第i个点的1阶差分值，D_i+1为归一化处理后的第i+1个点，D_i-1为归一化处理后的第i-1个点；

S32：二阶差分：

其中D″_i为第i个点二阶差分值，D′_i+2为一阶差分值中的第i+2个点的值，D′_i-2为一阶差分值中的第i-2个点的值；

S33：三阶差分：

D″′_i＝D″_i+1-D″_i，i∈[4，36]，其中D″′_i为第i个点的三阶差分值，D″_i+1为二阶差分中的第i+1个点的值；

S4：取二阶差分最大值，选出三阶差分中符合“正负负负”或“正负负”或“正负”规律的最大点，且要求对应循环数的一阶差分大于0，该坐标+1即为二阶差分最大值对应坐标，另±1共计3点；

S5：多项式拟合，将上述步骤S4得到的3点按二次函数y＝ax²+bx+c拟合，取

其中a是二次项系数，b是一次项系数，c是常数；

S6：切线求Ct值，对D_i进行曲线拟合，并过X处做D_i拟合曲线的切线，其与基线交点的横坐标即为Ct值，进而可得阈值，

S7：阴/阳性判别：

若Ct值≥38，判断为阴性；反之，则为阳性。

2.如权利要求1所述PCR检测方法的数据处理方法，其特征在于，还包括步骤S8：绘图，将D_i按循环数绘制曲线。

3.如权利要求1所述PCR检测方法的数据处理方法，其特征在于，以Savitzky-Golay平滑算法函数sg_smooth对去除基底的数据B_i进行平滑处理，其中参数为窗口大小：3，阶数：6。