1.一种基于跨平台基因表达数据的基因调控网络重建方法,其特征在于,包括以下步骤:S1)、从GEO基因表达数据库中获取来自p个测序平台基因表达数据,其中,每个测序平台包括多个基因样本,分别对p个测序平台基因表达数据进行预处理,删除每个测序平台中基因表达量缺失过多的基因样本,并将预处理后的p个测序平台基因表达数据合并得到基因样本量为m的基因样本集X={x1,x2,...xm},然后分别从基因样本集X={x1,x2,...xm}的每个基因样本中抽取n个基因表达量,从而得到基因调控网络样本集G={g11,g12,...g1n,g21,g22,...g2n,…gm1,gm2,...gmn};
S2)、初始化基因调控网络,用因果网络图表示基因调控网络,其中,一个基因表达量对应图中的一个节点,平台变量用一个对所有基因表达量都有影响的外生节点d表示,将基因调控网络初始化为空图;
S3)、用混合型条件独立性测试检验基因表达量gi与其他各基因表达量gj之间的条件独立性,若条件独立性测试结果表明gi和gj没有被其他基因表达量d-分隔(有向分隔),则将gj加入到gi的父子节点集PC(gi)中,即PC(gi)={gj};
S4)、重复步骤S3),得到基因调控网络中每个基因表达量节点的父子节点集PC(gi);
S5)、根据各基因表达量的父子节点集PC(gi),若基因调控网络中的任意两个基因表达量gi和gj互相存在于对方的父子节点集中,则将gi和gj用无向边连接起来,得到基因调控网络图骨架;
S6)、逐个检查网络图骨架中具有gi-gj-gk结构的任意三个基因表达量节点,根据混合型条件独立性测试的结果,确定是否存在v-结构,若混合型条件独立性测试的结果为gi⊥gk, 则认为gi,gj,gk之间构成了一个v-结构,则将图中这两条无向边的方向标记为gi→gj←gk;
S7)、根据不产生多余的v-结构、无环约束规则对基因调控网络图中未标记方向的边进行最大化方向标记,直到不能标记更多无向边的方向为止;
不产生多余的v-结构的规则为:网络图中具有gi→gj-gk结构的任意三个基因表达量节点,将gj-gk之间的边标记为gj→gk;
无环约束的规则为:网络图中两个相连节点gi-gj,若存在gi指向gj的有向路径,则将gi-gj标记为gi→gj;
S8)、输出跨平台基因调控网络。
2.根据权利要求1所述的一种基于跨平台基因表达数据的基因调控网络重建方法,其特征在于:步骤S3)中,所述的混合型条件独立性测试包括以下步骤:S301)、判断基因表达量gi,gj在给定基因表达量gK和平台时的条件独立性,用d表示平台变量,将gi,gj和gK的样本数据抽取出来;
S302)、对每个平台测序下对应的gi,gj,gK的样本数据,计算偏相关系数ρi,j|K,偏相关系数ρi,j|K可以通过线性回归计算得到,或者递归计算下面的公式得到,其中h∈K:S303)、重复S302步骤,分别得到p个测序平台对应的偏相关系数集Pi,j|K={(ρi,j|K)1,(ρi,j|K)2,...,(ρi,j|K)p};
S304)、若这p个偏相关系数总体上为零,则gi⊥gj|[gK,d]成立,然后对偏相关系数集Pi,j|K={(ρi,j|K)1,(ρi,j|K)2,...,(ρi,j|K)p}的p个偏相关系数进行费雪z-转换,得到Zi,j|K={z1,z2,...zp},其转换的公式为:S305)、若 成立,则基因表达量gi,gj在给
定基因变量gK和平台时条件独立,则表明基因表达量gi,gj被基因变量gKd-分隔(有向分隔);
其中,Nq为第q个平台的样本数,Z为gK的维度,Φ-1(1-α/2,p)表示均值为0,均方差为p的正态分布累计函数的反函数。