1.本发明涉及生物信息学领域，特别是涉及一种微卫星稳定性的检测方法、检测装置、计算机设备及存储介质。
背景技术：2.微卫星不稳定性(microsatellite instability，msi)是dna复制及损伤过程中出现的碱基错配、未配或多配造成的现象，表现为同一微卫星位点在不同个体之间的差异，或者同一个体的正常组织与某些异常组织之间重复单位的差异。
3.由一系列特异性修复dna碱基错配的酶组成的dna错配修复(mismatch repair，mmr)系统能够查出msi并进行修复，保证复制的准确性。如果mmr修复系统出现故障，如参与mmr修复的基因发生了突变，则有可能导致msi和整个基因组的不稳定性，最终可能造成细胞增殖分化异常和肿瘤的发生。
4.通过msi检测，可以评价mmr系统功能，从而预测癌症患者是否对免疫疗法敏感。目前msi相关的检测方法主要有3种：
5.1)免疫组化(immunohistochemistry，ihc)：通过使用相应抗体检测四种常见错配修复基因(mlh1、msh2、msh6和pms2)的表达情况，明确是否存在错配修复缺陷。存在1种及以上蛋白表达阴性即为错配修复缺陷(deficient mismatch repair，dmmr)，否则为错配修复蛋白完整(proficient mismatch repair，pmmr)。其优势在于应用性广泛，并且能确定哪些mmr蛋白在肿瘤细胞中表达缺失。但是ihc本身存在主观性，同时受抗体质量和实验因素等影响，有时无法检出某些蛋白的变化，导致mmr结果偶有报错。
6.2)基于pcr(polymerase chain reaction)的方法：主要采用多重荧光pcr结合毛细管电泳的方法，通过pcr扩增特定的微卫星序列，然后通过毛细管电泳比较肿瘤组织与正常组织微卫星序列长度的差异来判断该位点是否存在msi现象。这种检测方法是公认的msi检测的金标准，也是使用最广泛的方法。pcr检测方法不仅弥补了ihc在因非截断式错义突变导致的msi无法检出的漏洞，同时还具备良好的可重复性。但是，其检测的基因(panel)的位点较少、通量较低、无法提供具体的基因突变信息，而且实验周期较长。
7.3)基于ngs(next-generation sequencing)的方法：与pcr相比，ngs方法通量大，涉及基因范围广、灵敏度和特异性更高，可与靶点的突变检测和肿瘤突变负荷(tumor mutation burden，tmb)检测等共用一份测序数据。在目前已发表的ngs方法中，一般都是以pcr检测结果作为金标准，通过比较二者结果一致性作为评价ngs检测性能的标准。但目前基于ngs检测方法都需要以正常样本为参照，增加了取样步骤。
技术实现要素：8.基于此，有必要提供一种微卫星稳定性的检测方法，以能够实现对样本微卫星稳定性的无参且较为准确的检测。
9.此外，还提供一种微卫星稳定性的检测装置、一种计算机设备和一种计算机存储
介质。
10.一种微卫星稳定性的检测方法，包括以下步骤：
11.获取待测样本的测序数据；
12.根据待测样本的测序数据，得到待测样本的突变基因状态矩阵其中，t为所有待测样本的待分析的基因数量，q＝(b
1t
，b
2t
，b
3t
，...，b
nk
)
t
，b为各待测样本的目标位点的突变状态，n为待测样本的样本量，k＝1，2，3，...，t；及
13.根据预设的系数β和公式得到yv，yv＝(d1，d2，d3，...，dn)
t
，d代表得到的各待测样本的微卫星稳定性状态，n为待测样本的样本量。
14.在其中一个实施例中，上述微卫星稳定性的检测方法还包括如下确定预设的系数β的步骤：
15.获取已确定微卫星稳定性状态的训练样本的测序数据，作为训练集测序数据；
16.根据训练集测序数据，得到训练样本的突变基因状态矩阵其中，i为所有训练样本的待分析的突变基因数量，p＝(a
1j
，a
2j
，a
3j
，...，a
mj
)
t
，a为各训练样本的目标位点的突变状态，m为训练样本的样本量，j＝1，2，3，...，i；及
17.根据由弹性网络回归得到的关系式得到预设的系数β，其中，y＝(c1，c2，c3，...，cm)
t
，c代表各训练样本的微卫星稳定性状态，m为训练样本的样本量。
18.在其中一个实施例中，λ1和λ2的计算包括步骤：令λ＝λ1＝λ2，采用5-折交叉验证估计得到λ。
19.在其中一个实施例中，上述检测方法中待测样本的测序数据和训练集测序数据是通过二代测序得到的。
20.一种微卫星稳定性的检测装置，包括以下模块：
21.测序数据获取模块，用于获取待测样本的测序数据；
22.突变状态获取模块，用于根据待测样本的测序数据，得到待测样本的突变基因状态矩阵其中，t为所有待测样本的待分析的基因数量，q=(b
1t
，b
2t
，b
3t
，...，b
nk
)
t
，b为各待测样本的目标位点的突变状态，n为待测样本的样本量，k＝1，2，3，...，t；及
23.微卫星稳定性状态检测模块，用于根据预设的系数β和公式得到yv，yv＝(d1，d2，d3，...，dn)
t
，d代表得到的各待测样本的微卫星稳定性状态，n为所述待测样本的样本量。
24.在其中一个实施例中，上述微卫星稳定性状态的检测装置还包括确定预设的系数β的模块，确定所述预设的系数β的模块还包括数据获取模块、状态获取模块和系数确定模块；其中数据获取模块用于获取已确定微卫星稳定性状态的训练样本的测序数据，作为训练集测序数据；状态获取模块用于根据所述训练集测序数据，得到训练样本的突变基因状
态矩阵其中，i为所有训练样本的待分析的突变基因数量，p＝(a
1j
，a
2j
，a
3j
，...，a
mj
)
t
，a为各训练样本的目标位点的突变状态，m为训练样本的样本量，j＝1，2，3，...，i；系数确定模块用于根据由弹性网络回归得到的关系式得到预设的系数β，其中，y＝(c1，c2，t3，...，cm)
t
，c代表各训练样本的微卫星稳定性状态，m为训练样本的样本量。
25.在其中一个实施例中，上述系数确定模块中的λ1和λ2的计算包括步骤：令λ＝λ1＝λ2，采用5-折交叉验证估计得到λ。
26.在其中一个实施例中，上述检测装置中待测样本的测序数据和训练集测序数据是通过二代测序得到的。
27.一种计算机设备，具有处理器和存储器，存储器上存储有计算机程序，所述处理器执行所述计算机程序时实现如上述任一实施例中的微卫星稳定性的检测方法的步骤。
28.一种计算机存储介质，其上存储有计算机程序，计算机程序被执行时实现如上述任一实施例中的微卫星稳定性的检测方法的步骤。
29.上述微卫星稳定性的检测方法和检测装置，结合统计模型和相关基因位点，能够实现在没有阴性参照样本的情况下对微卫星稳定性状态进行评估，较传统的检测方法，提高了效率，降低了检测成本，提高了准确性。
具体实施方式
30.为使本发明的上述目的、特征和优点能够更加明显易懂，下面对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似改进，因此本发明不受下面公开的具体实施例的限制。
31.除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
32.本技术一实施方式提供了一种微卫星稳定性的检测方法，包括步骤s11、步骤s12和步骤s13。
33.具体地，步骤s11：获取待测样本的测序数据。
34.在一个可选的具体示例中，通过二代测序(next generation sequencing，ngs)得到待测样本测序数据。具体地，采用二代测序能够方便快捷地获得大量样本的大量基因的测序数据，由于涉及的基因数量多，能够提高微卫星稳定性检测的灵敏度和准确性。可以理解的是，也可以采用其他测序方法得到待测样本测序数据。
35.步骤s12：根据待测样本的测序数据，得到待测样本的突变基因状态矩阵其中，t为所有待测样本的待分析的基因数量，q＝(b
1t
，b
2t
，b
3t
，...，b
nk
)
t
，b为各待测样本的目标位点的突变状态，n为待测样本的样本量，k＝1，2，3，...，t。
36.具体地，根据待测样本的测序数据，得到待测样本的突变基因状态矩阵，即是将待
测样本的测序数据分析得到各待测样本的各目标位点的突变状态，形成突变基因状态矩阵。例如，使用gatk(the genome analysis toolkit)或lianti等软件可以获得各待测样本的各目标位点的突变状态(突变/未突变)，如果有十个待测样本，每个待测样本包含有cdh1、egfr、kmt2d、nfe2l2、rb1、spta1和top1七个基因，每个基因检测一个目标位点，可以得到每个待测样本的基因突变结果如表1所示：
37.表1
[0038] cdh1egfrkmt2dnfe2l2rb1spta1top1样本01未突变未突变未突变未突变未突变未突变未突变样本02未突变未突变突变未突变突变未突变突变样本03未突变未突变未突变未突变未突变未突变未突变样本04突变未突变未突变突变未突变突变未突变样本05未突变未突变未突变未突变未突变未突变未突变样本06未突变未突变突变未突变未突变未突变未突变样本07未突变突变未突变未突变未突变未突变突变样本08未突变未突变未突变未突变未突变未突变未突变样本09突变未突变未突变突变未突变未突变未突变样本010未突变未突变未突变未突变未突变突变未突变
[0039]
根据表1中的突变结果，即能得到如下突变基因状态矩阵：
[0040][0041]
步骤s13：根据预设的系数β和公式得到yv，yv＝(d1，d2，d3，...，dn)
t
，d代表得到的各待测样本的微卫星稳定性状态，n为待测样本的样本量。
[0042]
具体地，h为计算得到的各待测样本微卫星稳定性的or值(odds ratio)，微卫星稳定性状态通过d的数值判断，当d等于1时，表示该待测样本为微卫星高不稳定状态(msi-h)；当d等于0时，表示该待测样本为非微卫星高不稳定状态，即包括微卫星低不稳定状态(msi-l)和微卫星稳定状态(mss)。由于微卫星高不稳定状态与对免疫疗法的预后有较强相关性，而微卫星低不稳定状态与免疫疗法的预后相关性较弱，因此本技术中将微卫星高不稳定状态与微卫星低不稳定状态区分开。
[0043]
在一些实施例中，上述微卫星稳定性的检测方法还包括步骤s131、步骤s132和步骤s133。
[0044]
具体地，步骤s131：获取已确定微卫星稳定性状态的训练样本的测序数据，作为训
练集测序数据。
[0045]
在一个可选的具体示例中，采用与步骤s11相同的测序方法进行测序，具体地，采用相同的测序方法能够获得相同位点的测序数据，并保证测序误差的一致性，能够使检测结果更准确。
[0046]
在一个可选的具体示例中，可以采用多重荧光pcr-毛细管电泳法确定训练样本的微卫星稳定性状态。可以理解的是，也可以采用其他方法确定训练样本的微卫星稳定性状态，例如免疫组化法。
[0047]
在一个可选的具体示例中，获取的微卫星高不稳定状态的训练样本与非微卫星高不稳定状态的训练样本的数量较大且相近，保持两者样本数量相近可以提高检测的准确度，避免训练得到的模型对其中一类样本产生偏好，可以有效降低假阳性和假阴性结果的出现。
[0048]
步骤s132：根据训练集测序数据，得到上述训练样本的突变基因状态矩阵其中，i为所有上述训练样本的待分析的突变基因数量，p＝(a
1j
，a
2j
，a
3j
，...，a
mj
)
t
，a为各训练样本的目标位点的突变状态，m为训练样本的样本量，j＝1，2，3，...，i。
[0049]
具体地，该步骤中得到训练样本的突变基因状态矩阵的方法与步骤s11相同。
[0050]
步骤s133：根据由弹性网络回归得到的关系式得到预设的系数β，其中，y＝(c1，c2，c3，...，cm)
t
，c代表各训练样本的微卫星稳定性状态，m为训练样本的样本量。
[0051]
具体地，c为0或1，是根据训练样本已知的微卫星稳定性状态设定，若训练样本为微卫星高不稳定状态则记为1，否则记为0。
[0052]
在一个可选的具体示例中，λ1和λ2的计算包括步骤：令λ＝λ1＝λ2，采用5-折交叉验证估计得到λ。可以理解的是，λ1和λ2为正则化系数，在其他一些实施例中，可以采用其他方法计算得出λ1和λ2的值，能实现最终计算出β的值即可。
[0053]
可以理解的是，在经过步骤s131～步骤s133后，得到的β能直接用于步骤s13中微卫星稳定性的检测，不需要每次检测之前都重新确定β的值。因此，能够实现只取待测样本进行微卫星稳定性的检测，而不需要取正常样本作为阴性对照，实现了无参的微卫星稳定性检测。
[0054]
此外，经过步骤s131～步骤s133，同时也实现了特征筛选，即筛选出了与微卫星不稳定相关的基因，也能够将筛选得到的这些基因用在与微卫星不稳定相关的其他研究上。
[0055]
基于与上述微卫星稳定性的检测方法相同的思想，本技术一实施方式还提供了一种微卫星稳定性的检测装置，包括模块21、模块22和模块23。
[0056]
具体地，模块21：测序数据获取模块，用于获取待测样本的测序数据。
[0057]
模块22：突变状态获取模块，用于根据待测样本的测序数据，得到待测样本的突变基因状态矩阵其中，t为所有待测样本的待分析的基因数量，q＝(b
1t
，b
2t
，b
3t
，...，b
nk
)
t
，b为各待测样本的目标位点的突变状态，n为待测样本的样本量，k＝1，2，3，...，t。
[0058]
模块23：微卫星稳定性状态检测模块，用于根据预设的系数β和公式得到yv，yv＝(d1，d2，d3，...，dn)
t
，d代表得到的各待测样本的微卫星稳定性状态，n为待测样本的样本量。
[0059]
在一些实施例中，上述微卫星稳定性的检测装置还包括模块231、模块232和模块233。
[0060]
具体地，模块231：数据获取模块，用于获取已确定微卫星稳定性状态的训练样本的测序数据，作为训练集测序数据。
[0061]
模块232：状态获取模块，用于根据训练集测序数据，得到训练样本的突变基因状态矩阵其中，i为所有训练样本的待分析的突变基因数量，p＝(a
1j
，a
2j
，a
3j
，...，a
mj
)
t
，a为各训练样本的目标位点的突变状态，m为训练样本的样本量，j＝1，2，3，...，i。
[0062]
模块233：系数确定模块，用于根据由弹性网络回归得到的关系式得到预设的系数β，其中，y=(c1，c2，c3，...，cm)
t
，c代表各训练样本的微卫星稳定性状态，m为训练样本的样本量。
[0063]
可以理解的是，在采用已确定微卫星稳定性状态的训练样本进行训练时需要调用模块231～模块233，在训练完成后可直接将β的值输入微卫星稳定性的检测装置的模块23中，可直接用于对待测样本的微卫星稳定性的检测，不需要每次检测都调用模块231～模块233进行训练。
[0064]
基于如上所述的各实施例，本发明还提供了一种计算机设备，具有处理器和存储器，存储器上存储有计算机程序，所述处理器执行所述计算机程序时实现如上述任一实施例和具体示例中的微卫星稳定性的检测方法的步骤。
[0065]
基于如上所述的各实施例，本发明还提供了一种计算机存储介质，其上存储有计算机程序，计算机程序被执行时实现如上述任一实施例和具体示例中的微卫星稳定性的检测方法的步骤。
[0066]
本领域普通技术人员可以理解实现上述方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性的计算机可读取存储介质中，该计算机程序在执行时，可包括如上各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其他介质的任何引用，均包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(read-only memory，rom)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(random access memory，ram)或外部高速缓冲存储器。作为说明而非局限，ram可以是多种形式，比如静态随机存取存储器(static random access memory，sram)或动态随机存取存储器(dynamic random access memory，dram)等。
[0067]
上述微卫星稳定性的检测方法和检测装置，结合统计模型和相关基因位点，能够实现在没有阴性参照样本的情况下对微卫星稳定性状态进行评估，较传统的检测方法，提高了效率，降低了检测成本，提高了准确性。
[0068]
具体实施例
[0069]
以下采用具体实施例说明本发明的微卫星稳定性的检测方法产生的技术效果。
[0070]
实施例1
[0071]
1.获取训练样本
[0072]
选取530个已确定微卫星稳定性状态的肿瘤样本作为训练样本，依次记为样本1～样本530。其中共270例微卫星高不稳定状态(msi-h)样本，260例非微卫星高不稳定状态(msi-l和mss)样本，并将微卫星高不稳定状态记为1，否则记为0，并记成y＝(c1，c2，c3，...，c
530
)
t
的形式，其中，c为各样本的微卫星稳定性状态，即(由于样本数量大，以下矩阵y中没有列出全部样本的微卫星未定型状态，只列举部分以示例)：
[0073][0074]
2.获取训练样本突变基因的状态
[0075]
采用二代测序对训练样本进行测序，获取测序数据，根据测序结果，将各目标位点的突变状态分为突变与未突变，将突变记为1，为突变记为0，如表2所示，并记成突变基因状态矩阵其中，i为所有训练样本的待分析的突变基因数量，i=17105，p＝(a
1j
，a
2j
，a
3j
，...，a
mj
)
t
，a为各训练样本的目标位点的突变状态，m为样本量，m=530，j＝1，2，3，...，i。由于样本和待分析的突变基因数量大，表2以及矩阵x中没有全部列出，只列举部分以示例。
[0076]
表2
[0077] cdk8cebpaepha7
…
fgf21gata4jak3样本1100
…
000样本2001
…
010样本3000
…
000
……………………
样本528010
…
000样本529100
…
001样本530000
…
100
[0078]
将表2的内容记成突变基因状态矩阵为：
[0079][0080]
3.通过训练确定λ并计算出β
[0081]
令λ＝λ1＝λ2，采用5-折交叉验证可估计得到λ。之后根据由弹性网络回归得到的关
系式计算得出β的值。其中，y和x分别为以上步骤1和步骤2中得到的。
[0082]
4.验证本发明检测方法的准确性
[0083]
将步骤1中的训练样本采用本发明的方法检测，即将步骤2中的矩阵x代入得到微卫星稳定性检测结果，与原来已确定的微卫星稳定性结果进行比对，如表3所示。
[0084]
表3
[0085][0086]
表3中，表3中，tp：真阳性样本数；tn：真阴性样本数；fp：假阳性样本数；fn：假阴性样本数。
[0087]
从表3可以看出，本发明的微卫星稳定性状态的检测方法的准确率、敏感性和特异性都很高，说明该方法能够用于微卫星稳定性状态的检测，在没有正常样本做阴性对照的条件下也能准确地确定样本的微卫星稳定性状态，且检测通量大，效率高。
[0088]
5.检测待测样本的微卫星不稳定状态
[0089]
获取40个待测样本，依次命名为待测样本1～待测样本40。采用二代测序对待测样本进行测序，获取测序数据，根据测序结果，将各目标位点的突变状态分为突变与未突变，将突变记为1，为突变记为0，如表4所示，并记为突变基因状态矩阵其中，t为所有待测样本的待分析的基因数量，t＝17105，q＝(b
1t
，b
2t
，b
3t
，...，b
nk
)
t
，b为各待测样本的目标位点的突变状态，n为待测样本的样本量，n＝40，k＝1，2，3，...，t。由于样本量和基因数量较大，以下列举出部分结果以示例：
[0090]
表4
[0091] cdk8cebpaepha7
…
fgf21gata4jak3待测样本1000
…
000待测样本2110
…
000待测样本3000
…
001
……………………
待测样本38001
…
100待测样本39100
…
000待测样本40010
…
000
[0092]
将表4的内容记成突变基因状态矩阵为：
[0093][0094]
将步骤3中计算得到的β代入中计算得到yv，yv＝(d1，d2，d3，...，dn)
t
，d代表得到的各待测样本的微卫星稳定性状态，微卫星稳定性状态通过d的数值判断，当d等于1时，表示该待测样本为微卫星高不稳定状态(msi-h)；当d等于0时，表示该待测样本为非微卫星高不稳定状态，即包括微卫星低不稳定状态(msi-l)和微卫星稳定状态(mss)。以下列举出部分结果以示例：
[0095][0096]
根据yv，各待测样本对应的微卫星稳定性状态如表5所示(只列举部分以示例)。
[0097]
表5
[0098]
样本名称微卫星稳定性状态待测样本1msi-l或mss待测样本2msi-l或mss待测样本3msi-l或mss
……
待测样本38msi-l或mss待测样本39msi-l或mss待测样本40msi-h
[0099]
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0100]
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。应当理解的是，在本领域技术人员在本发明提供的技术方案的基础上，通过合乎逻辑的分析、推理或有限的试验得到的技术方案，均在本发明所附权利要求的保护范围内。因此，本发明专利的保护范围应以所附权利要求为准，说明书可以用于解释权利要求的内容。
技术特征：1.一种微卫星稳定性的检测方法，其特征在于，包括以下步骤：获取待测样本的测序数据；根据所述待测样本的测序数据，得到所述待测样本的突变基因状态矩阵，其中，t为所有所述待测样本的待分析的基因数量，q＝(b
1t
，b
2t
，b
3t
，...，b
nk
)
t
，b为各所述待测样本的目标位点的突变状态，n为所述待测样本的样本量，k＝1，2，3，...，t；及根据预设的系数β和公式得到y
v
，y
v
＝(d1，d2，d3，...，d
n
)
t
，d代表得到的各所述待测样本的微卫星稳定性状态，n为所述待测样本的样本量。2.根据权利要求1所述的检测方法，其特征在于，还包括如下确定所述预设的系数β的步骤：获取已确定微卫星稳定性状态的训练样本的测序数据，作为训练集测序数据；根据所述训练集测序数据，得到所述训练样本的突变基因状态矩阵其中，i为所有所述训练样本的待分析的突变基因数量，p＝(a
1j
，a
2j
，a
3j
，...，a
mj
)
t
，a为各所述训练样本的目标位点的突变状态，m为所述训练样本的样本量，j＝1，2，3，...，i；及根据由弹性网络回归得到的关系式得到所述预设的系数β，其中，y＝(c1，c2，c3，...，c
m
)
t
，c代表各所述训练样本的微卫星稳定性状态，m为所述训练样本的样本量。3.根据权利要求2所述的检测方法，其特征在于，λ1和λ2的计算包括步骤：令λ＝λ1＝λ2，采用5-折交叉验证估计得到λ。4.根据权利要求2或3所述的检测方法，其特征在于，所述待测样本的测序数据和所述训练集测序数据是通过二代测序得到的。5.一种微卫星稳定性的检测装置，其特征在于，包括以下模块：测序数据获取模块，用于获取待测样本的测序数据；突变状态获取模块，用于根据所述待测样本的测序数据，得到所述待测样本的突变基因状态矩阵其中，t为所有所述待测样本的待分析的基因数量，q＝(b
1t
，b
2t
，b
3t
，...，b
nk
)
t
，b为各所述待测样本的目标位点的突变状态，n为所述待测样本的样本量，k＝1，2，3，...，t；及微卫星稳定性状态检测模块，用于根据预设的系数β和公式得到y
v
，y
v
＝(d1，d2，d3，...，d
n
)
t
，d代表得到的各所述待测样本的微卫星稳定性状态，n为所述待测样本的样本量。6.根据权利要求5所述的检测装置，其特征在于，还包括确定所述预设的系数β的模块，所述确定所述预设的系数β的模块还包括数据获取模块、状态获取模块和系数确定模块；所述数据获取模块用于获取已确定微卫星稳定性状态的训练样本的测序数据，作为训练集测序数据；所述状态获取模块用于根据所述训练集测序数据，得到所述训练样本的突变基因
状态矩阵其中，i为所有所述训练样本的待分析的突变基因数量，p＝(a
1j
，a
2j
，a
3j
，...，a
mj
)
t
，a为各所述训练样本的目标位点的突变状态，m为所述训练样本的样本量，j＝1，2，3，...，i；所述系数确定模块用于根据由弹性网络回归得到的关系式得到所述预设的系数β，其中，y＝(c1，c2，c3，...，c
m
)
t
，c代表各所述训练样本的微卫星稳定性状态，m为所述训练样本的样本量。7.根据权利要求6所述的确定装置，其特征在于，λ1和λ2的计算包括步骤：令λ＝λ1＝λ2，采用5-折交叉验证估计得到λ。8.根据权利要求6或7所述的微卫星稳定性的检测装置，其特征在于，所述待测样本的测序数据和所述训练集测序数据是通过二代测序得到的。9.一种计算机设备，其特征在于，具有处理器和存储器，所述存储器上存储有计算机程序，所述处理器执行所述计算机程序时实现如权利要求1～4任一项所述的微卫星稳定性的检测方法的步骤。10.一种计算机存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被执行时实现如权利要求1～4任一项所述的微卫星稳定性的检测方法的步骤。
技术总结本发明涉及一种微卫星稳定性的检测方法、检测装置、计算机设备及存储介质。上述微卫星稳定性的检测方法和检测装置，结合统计模型和相关基因位点，能够实现在没有阴性参照样本的情况下对微卫星稳定性状态进行评估，较传统的检测方法，提高了效率，降低了检测成本，提高了准确性。准确性。
技术研发人员：王维锋 刘成琨 郑新 姚继成
受保护的技术使用者：上海至本医学检验所有限公司
技术研发日：2022.03.15
技术公布日：2022/7/5
