[发明专利]一种基于双向XGBoost的基因调控网络构建方法

权利要求书

1.一种基于双向XGBoost的基因调控网络构建方法，其特征在于，包括以下步骤：

步骤1、根据时序基因表达数据，依次选取每一个基因作为特定基因，针对每一个特定基因分别构建一个双向模型，一个双向模型包括一个local-in模型和一个local-out模型，在local-in模型中，特定基因为调控目标基因，其他基因为其候选调控基因；在local-out模型中，特定基因为调控基因，其他基因为其候选调控目标基因；

local-in模型表示为以下回归方程：

其中，表示除特定基因G_i以外的其他基因在t-p至t-1时刻的表达值集合，N表示时序基因表达数据中基因的个数，fⁱ( )表示特定基因G_i在t时刻的表达值与G^p之间存在的函数关系，p表示时滞系数，T为时序基因表达数据中时间片的个数，∈_i表示随机误差；该回归方程以其他基因在t-p至t-1时刻的表达值为自变量，以特定基因G_i在t时刻的表达值为因变量，表示其他基因在t-p至t-1时刻的表达值对于特定基因在t时刻的表达值的影响程度，用于为特定基因选择调控基因，在基因调控网络中产生调控基因-G_i这样的有向边，-表示调控关系；

local-out模型表示为以下回归方程：

其中，表示除特定基因G_i以外的其他基因在t+1至t+p时刻的表达值集合；该回归方程以其他基因在t+1至t+p时刻的表达值为自变量，以特定基因G_i在t时刻的表达值为因变量，表示特定基因在t时刻的表达值对于其他基因在t+1至t+p时刻的表达值的影响程度，用于为特定基因选择调控目标基因，在基因调控网络中产生G_j-调控目标基因这样的有向边；t+p＜T；

步骤2、使用XGBoost对回归方程建立梯度提升回归树；在回归树的建立过程中，以回归方程的因变量，即特定基因G_i在t时刻的表达值的样本空间为根节点，以回归方程的自变量，即其他基因的表达值作为特征，在每个节点上，选择带惩罚的方差降低最多的特征作为分裂特征，根据分裂特征的布尔函数对节点进行划分；对每个特征统计其在建立回归树的过程中被选择作为分裂特征的次数，作为该特征的重要性打分，并将其作为特定基因G_i与该特征对应的候选基因之间的调控关系的初步权重值；

步骤3、对每个独立的回归方程计算残差值，并对残差值进行归一化处理；将归一化处理后的残差值作为评价每个独立的回归方程的权重，将每个独立的回归方程的权重乘以每个回归方程中调控关系的初步权重值作为每个回归方程中的调控关系的最后打分；最终结合双向模型中同一对调控关系的两次打分进行最终打分，得到该对调控关系的全局权重；

步骤4、对调控关系的全局权重由大到小进行排序，选取排名靠前的M个调控关系，在相应基因之间建立一条边，从而形成基因调控网络。

2.根据权利要求1所述的基于双向XGBoost的基因调控网络构建方法，其特征在于，所述步骤3中残差值为均方误差，回归方程中的调控关系的最后打分通过以下方法计算：

对任意一个基因G_i建立的回归方程计算均方误差SSE_i，并进行归一化，得到nSSE_i：

其中，表示特定基因G_i在t时刻的表达值表示由回归树得到的估计值；

并采用以下方式计算基因G_i与基因G_j之间的调控关系的最后打分：

I(i,j)＝nSSE_i·nw_i,j

其中，nw_i,j表示调控关系的初步权重值w_i,j归一化后的值，nw_i,j＝w_i,j/∑_alljw_i,j。