[发明专利]父母子三元亲属结构的候选区域基因关联检测系统及方法

权利要求书

1.父母子三元亲属结构的候选区域基因关联检测系统，其特征在于：所述系统包括：信息获取模块，主成分提取模块，家庭内部相关性统计量构造模块，全部家庭数据关联检测统计量获取模块，基因关联检测模块；

所述信息获取模块用于获取父母子三元亲属结构家庭成员基因组上的标记位点信息x_ijl以及性状信息y_ijl；

其中，x_ijl和y_ijl分别表示第i个子群体第j个家庭中第l个个体的基因型得分和数量性状值；

所述主成分提取模块用于提取信息获取模块获取的信息主成分t_ijl＝(t_ijl1,…,t_ijlf)；

其中，f是主成分的总数量；

所述家庭内部相关性统计量构造模块用于根据x_ijl和y_ijl和t_ijl构造家庭内部相关性统计量U_ij：

步骤一、建立主成分t_ijl＝(t_ijl1,…,t_ijlf)与y_ijl的广义线性函数和主成分t_ijl＝(t_ijl1,…,t_ijlf)与候选标记点位信息x_1ij^l的广义线性函数：

y_ijl＝g₁(t_ijl)+ψ_ijl

x_1ijl＝g₂(t_ijl)+ε_ijl

其中，ψ_ijl是随机变量，ε_ijl是随机变量，g₁(t_ijl)和g₂(t_ijl)是模型核心函数；

步骤二、利用步骤一获取的广义线性函数获取家庭内部的相关性统计量U_ij；

所述全部家庭数据关联检测统计量获取模块用于根据U_ij获取全部家庭数据关联检测统计量；

所述基因关联检测模块用于根据获取全部家庭数据关联检测统计量的统计P值，并根据统计P值确定基因与性状信息是否有关联。

2.根据权利要求1所述的父母子三元亲属结构的候选区域基因关联检测系统，其特征在于：

其中，β₁＝(β₁₀,β₁₁…,β_1f)是模型核心函数g₁(t_ijl)中主成分的系数；β₂＝(β₂₀,β₂₁…,β_2f)是模型核心函数g₂(t_ijl)中主成分的系数。

3.根据权利要求2所述的父母子三元亲属结构的候选区域基因关联检测系统，其特征在于：E(ψ_ijl)＝E(ε_ijl)＝0；

其中，E(ψ_ijl)和E(ε_ijl)是随机变量ψ_ijl和ε_ijl的期望。

4.根据权利要求3所述的父母子三元亲属结构的候选区域基因关联检测系统，其特征在于：所述步骤二中的利用步骤一获取的广义线性函数获取家庭内部的相关性统计量U_ij：

步骤二一、基于广义估计方程理论对g₁(t_ijl)和g₂(t_ijl)进行估计获得的

步骤二二、将步骤二一获得的代入步骤一获得广义线性函数中获得和并利用和配置残差变量；

步骤二三、利用步骤二二获得的残差变量构造家庭内部的相关性统计量U_ij。