[发明专利]一种大规模连锁不平衡遗传位点快速分区方法

权利要求书

1.一种大规模连锁不平衡遗传位点快速分区方法，其特征在于，包括以下步骤：

S1、根据位点的物理位置选定待分区间，得到待分区间；

S2、根据待分区间选定滑动处理窗口，得到滑动处理窗口；

S3、根据滑动处理窗口对待分区间进行初分区处理，得到初分区块；

S4、将初分区块大小与预设值进行比较并根据比较结果确定连锁不平衡区块；

S5、根据预设规则移动滑动处理窗口继续对待分区间进行初分区处理并返回步骤S4，直至滑动处理窗口移至待分区间的末端跳出循环，分区完成；

所述滑动处理窗口设为50个位点大小，所述根据滑动处理窗口对待分区间进行初分区处理，得到初分区块这一步骤，其具体包括以下步骤：

S301、获取滑动处理窗口内的相关系数信息；

S302、将滑动处理窗口内的相关系数信息与预设阈值进行比较，当滑动窗口内的相关系数大于预设阈值的个数超过预设数量，滑动处理窗口向右移，当滑动窗口内的相关系数大于预设阈值的个数没有超过预设数量，滑动处理窗口向下移；

S303、更新移动后的滑动窗口内的相关系数信息并返回步骤S302，直至滑动窗口向下移动后窗口首个位点的行坐标与列坐标相等，跳出循环；

S304、以滑动处理窗口跳出循环前的搜索边界作为区块边界，得到一个初分区块；

所述将初分区块大小与预设值进行比较并根据比较结果确定连锁不平衡区块这一步骤，其具体包括：

S401、获取初分区块尺寸大小；

S402、将初分区块尺寸大小与预设值进行比较，若初分区块尺寸大小大于预设值，调整滑动处理窗口位置和阈值并对粗略分区块进行重分区，得到连锁不平衡区块，若初分区块尺寸大小小于预设值，确认该初分区块为连锁不平衡区块。

2.根据权利要求1所述一种大规模连锁不平衡遗传位点快速分区方法，其特征在于，所述根据位点的物理位置选定待分区间，得到待分区间这一步骤之前，还包括：

S0、对位点进行预处理，过滤掉无法计算相关系数的位点。

3.根据权利要求2所述一种大规模连锁不平衡遗传位点快速分区方法，其特征在于，所述待分区间的大小为第一个位点与最后一个位点距离小于5MB。

4.根据权利要求3所述一种大规模连锁不平衡遗传位点快速分区方法，其特征在于，所述调整滑动处理窗口位置和阈值并对初分区块进行重分区，得到连锁不平衡区块，其具体为：

调整滑动处理窗口的首位置到上一个不平衡区块的首个位点，增大滑动处理窗口的相关系数预设阈值并以步骤S3的方法对初分区块进行重新分区，直至得到的分区块大小小于预设值，得到连锁不平衡区块。

5.根据权利要求4所述一种大规模连锁不平衡遗传位点快速分区方法，其特征在于，所述根据预设规则移动滑动处理窗口继续对待分区间进行初分区处理，其具体包括：

在当前获得的不平衡连锁区块中的最后一个区块的首位点作为滑动处理窗口的起始位置进行下一个分区处理。