数学建模华中赛B题优秀论文正稿

1、专业资料第八届华中地区大学生数学建模邀请赛承诺书我们仔细阅读了第八届华中地区大学生数学建模邀请赛的竞赛细则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们的参赛报名号为：参赛队员（签名）：队员1:队员2:队员3:武汉工业与应用数学学会

2、第八届华中地区大学生数学建模邀请赛组委会word完美格式专业资料第八届华中地区大学生数学建模邀请赛编号专用页选择的题号：B参赛的编号：（以下内容参赛队伍不需要填写）竞赛评阅编号：word完美格式专业资料第八届华中地区大学生数学建模邀请赛题目：基因调控网络的重构及病毒感染的致病机制【摘要】一个基因的表达受其他基因的影响，而这个基因又影响其他基因的表达,这种相互影响相互制约的关系构成了复杂的基因调控网络。基因调控网络的研究是从基因之间相互作用的角度揭示复杂的生命现象，是当前生物信息学研究的前沿。疾病的发病因素和原理，对于医疗领域有着十分重要的作用。这不仅仅能够让更多的患者免受病痛的困扰，还能促进人

3、类医学史的进步。所以根据基因数据谱来重构基因调控网络，以及某个疾病症状产生的原因的研究具有很大的意义。本文对基因调控网络的重构以及导致严重临床症状的蛋白质进行了研究和推测。由于所给的基因数据谱（附录一）十分庞大，所以首先要对数据进行降维处理。本题基于时间序列给出了272组基因数据，为了减小噪声以及缺失值对实验精度的干扰，在实验前对四组噪声较大或有缺失的数据进行剔除。具体的降维方式采用了多元统计法中的主成分分析和聚类分析：先对这一万多个数据做主成分分析，从这一万多个数据中，通过线性变化选出了1000个左右的重要变量来组成新的样本。既降低了数据的处理难度，又尽量保持了新数据和原数据相比，尽可能保持

4、原数据的信息。然后用spss两阶聚类法粗略地对要聚类的数目进行一个估计，根据此彳&计用K-means算法对数据进行处理，得到相应的30组数据。对这30组数据建立模型，来重构基因调控网络。本文中采用的模型是线性回归模型，并对它的合理性，以及相较贝叶斯网络作了对比。最后依据所得到的系数矩阵进行基因网络图的绘制与呈现。问题二在第一问的基础上，寻求导致产生严重临床症状的蛋白质。根据附录二给出的个体出现感染症状时间节点示意图，1代表此志愿者在该时间节点表现出了临床症状，0则表示没有，这是一个二分类。本题采用逻辑回归模型，利用LR分类器模型去寻找该重要蛋白质。用268组数据，其中每一个基因视为该组

5、数据的一个属性，对这些基因进行LR分类，并得到相应的系数矩阵。然后对系数矩阵进行分析，取出影响比较大的几个基因，然后对照基因表对基因作用的描述去寻求该重要蛋白。本题最终找出四个导致志愿者产生严重的临床症状的蛋白质。所有代码实现，以及每次得到的系数矩阵均在附录中给出。关键词：线性回归模型，基因调控网络重构，多元统计法，主成分分析，聚类分析，逻辑回归（LR分类器模型）word完美格式专业资料1 .问题重述通过基因之间的相互调控，生物体可以实现细胞的生长，器官的发育、以及免疫等各种生物机能。随着测序技术的发展，产生了越来越多的高通量实验数据。基于这些实验数据重建基因调控网络（Generegulato

6、rynetworks,GRNS,对于深入了解生物机能的实现过程具有重要作用。生物实验中，在17个健康志愿者鼻内接种流感病毒H3N2/Wisconsin,其中9个人出现了严重的感染症状，另外的8个人没有出现症状。接种后，每隔大约8h从血液中采集样本测量基因表达谱数据，实验数据一共有16个时间点（单位：h）,包括baseline（-24）,0,5,12,21,29,36,45,53,60,69,77, 84,93,101,108,共268个样本。基因表达谱数据见附件1,其中前8个为未出现严重感染症状的数据，后9个为出现严重感染症状的数据。（其中行代表探针号，对应着不同的基因；列为各个个体血液样本在

7、各个时间节点的数据）个体出现感染症状的时间节点示意图见附件2。问题：1）根据实验数据重构基因调控网络；2）通过比较出现感染症状的志愿者和健康志愿者的样本数据，试确定病毒感染人体后导致志愿者是否会出现严重临床症状的重要蛋白。2 .问题分析一个基因的表达受其他基因的影响，而这个基因又影响其他基因的表达，这种相互影响相互制约的关系构成了复杂的基因调控网络。更一般些，几乎所有的细胞活动都被基因网络所控制。生命是存储并加工信息的复杂系统，孤立地研究单个基因及其表达往往不能确切地反映生命现象本身的内在规律。因此，需要从复杂系统的角度研究基因网络。对于问题一，考察我们如何根据已有的基因表达谱（附录一）去重构

8、基因调控网络，从而推断调控网络各节点之间潜在的调控关系。考虑“反向分析法”来重构基因调控网络，常见的基因调控网络模型有布尔网络模型、线性组合模型和贝叶斯网络模型等等。然而题目所给的数据集十分庞大，如果直接将这一万个基因全部带入模型，那么计算量是惊人的。所以需要用到多元统计方法中的主成分分析和聚类分析去实现降维的操作。对于问题二，在已经重构好的基因网络的基础上寻找导致病毒感染人体以后导致志援者是否产生严重临床症状的蛋白质。首先我们要对数据进行分析，寻找与染病相关系数大的基因，然后依据附录一的sheet2中对于基因的描述去进一步确定关键蛋白质。word完美格式专业资料3 .模型假设针对本问题，建立

9、如下合理假设：(1)题目所给数据准确可靠；(2)假设不考虑个体差异性；(3)基因表达呈高斯分布；4 .符号说明Xn,m表示第n个基因基于时间序列的第m组数据;AK表示一个基因；Pi为回归系数；at，x代表基因X在时间点t具有的表达值；bibh为常数；鸟，为误差项。5 .问题一的建模与算法实现求解5.1 数据的分析问题一需要根据所给的基因表达谱数据来重构基因调控网络，附录一中的sheetl中给出了17个志愿者体内的10000种基因，随着注入病毒后的时间变化而出现的数值变化。由于数据集过大，所以第一步要做的就是对这一万种基因进行筛选降维操作。只选取部分具有代表性的数据代入模型，从而减少计算量。对于

10、数据的处理部分，采用多元统计中的常用方法，主成分分析和聚类分析。5.2 数据预处理5.2.1 数据处理方法选择由于这道题目的数据量庞大，所以，如何筛选数据就成了很重要的一步。我们这里采取先对10000组数据做主成分分析，形成1000组新变量，再对这些新变量进行聚类分析，进一步降维。5.2.2 主成分分析(1)主成分分析的基本思想：主成分分析的基本思想是通过构造10000个基因初始数据的适当的线性组合，以产生一系列互不相关的新变量，从中选出少数几个新变量并使它们尽可能多地包含原先所有基因的信息(降维)，从而使得用这几个新变量替代原变量分word完美格式专业资料析问题成为可能。即在尽可能少丢失信息

11、的前提下从所研究的m个变量中求出几个新变量，它们能综合原有变量的信息，相互之间又尽可能不含重复信息。(2)主成分分析的实现：设有n个样品，m个变量(指标)的数据矩阵。本题中n=10000,表示10000种基因；m=268表示基于时间序列的基因数据变化指标。,XuXi2IIIXim1X(i)Xn加=X2iX22IIIX2mX(2)IIIIIIHIIII*.XniXn2IIIXnmJ1，(n)寻找k个新变量yi,y2,|,yk(kEm),使得1、yi=aiiXi+ai2X2十川+amXm,(l=1,2H,k)2、yi,y2,lllyk彼此不相关主成分的系数向量ai=(aii,ai2,|l|,aim

12、)的分量aij刻划出第j个变量关于第l个主成分的重要性。可以证明，若x=(x1,X2,山,Xm)T为m维随机向量，它的协方差矩阵V的m个特征值为儿之界|之*之0,相应的标准正交化的特征向量为Ui,U2,|,Um,则xKxnXzJIhXm)的第i主成分为小=UiTx(i=i,21|,m)。mkm称九/Z%为主成分yi=ux(i=i,2,|,m)的贡献率，工%ZZ%为主成分jij+jmyi,y2,Hlyk的累计贡献率，它表达了前k个主成分中包含原变量为？2,|,Xm的信息量大小，通常取k使累计贡献率在85%Z上即可。当然这不是一个绝对不变的标准，可以根据实际效果作取舍，例如当后面几个主成分的贡献率

13、较接近时，只选取其中一个就不公平了，若都选入又达不到简化变量的目的，那时常常将它们一同割舍。计算步骤如下：i、由已知的原始数据矩阵Xn不计算样本均值向量*X=(无亚2,出,鼠);i/其中Xi=-Xij(i=i,2,|l|,m)nj£2、计算样本协方差矩阵?=(Sj)?2ij)n-in其中Sij="(Xii-Xi)(Xj-Xj)(i,j=i,2,|l|,m)iiword完美格式专业资料3、把原始数据标准化，即ijx-x；j,记Xn刈=(Xij)o形成样本相关矩阵:二jj4、求口的特征根兀之九2至IM之3±0及相应的标准正交化的特征向量u1,u2,|,um,可得主成分

14、为yi=UiTx(i=1,2,|,m)o(3)主成分分析降维结果用Matlab实现以上算法(代码见附录)，实现结果如下:-new_gene田1000x272double1234557&910111-136.1692*145.3295»135.0&72162,1874-140,3731-156.7242»142.4655-141-5586J42.4953355.0893-12孤恻T4G.M1函.4£14140.4264SJ33143.S54442.01963-123.8926-112.0174122.4528*1131081d07.8305112,6

15、692-113.7195-114.2983-1153714-117.343607.0175421.3S01期434303011817.000219.197祖招拈20.717S2106623.7M327.4964543.1902*71惨4041104&.2Z9852,956747377545,2021150,80514&515647.056B49AB32e-121.4M4-116.443Q-1194194-iso.ass-117.466-118.9976-118.6342-120-3057-120.7&2S-119X3176-17*13.171214J471-13.120

16、3-13.2759d3.7&35-14.01B7=14,M21-9,M75-142739d2.30OB-11.0011g眼*59507964693963.119462a4前64.379363.5787晓51M9-3.95364.5757-2.3M1-&.O731415210-4,846-7J27J-7.9276-3.6125-5.8H14-9JD3471a76.7269711.3323744W73M5976.546472黜5876.334675.774476564277.512411-2.2W64.5444-1.56J3-2.TO34-2.7257-13422-1,B935-2

17、.5770=1.9712需8569-53B95u-41.774fi-47.71&7-414341-4147-434749-45-013-4.04lQ05M-4?.«W-4243471315.126318.029316.90®513,905018.1W719.719617,506015,208518.B23819,291415.96291456.577759-121557,910059.1588总黑器61.W97W.821560634060.0712航59360515dll.0318-109.4654-1O9.&2B9-1WJ241-112.5562-112.

18、2411-111.4234di0.4424-113-9180-112.2415-11O47&9116S5.304830.7075M122731.S913施.045433701235-610032.047134J77133.98W3S.91911726,4542252B1426,857524355826,247225,791025J22723.W7225,821124,430625605313M.14拓59367064.63096350759.585562.742163.636163.057S64.1177加.1热762.122219-20.0239用。出150-20.1272-19.3

19、2-11-16.70S4-1S,OT76-17.6315-19.563319此强-15.2937-1722a42D-5位通-4B.1W9-462B5-44.073-4682?照7725-43.5730-46宙第A-43-6W21-11.7029-114219-12.1946-114352-12.5736-12.2213-111.65011-10.854M-12J1S2-11.9341-12.7B81靛一能5460-4Q.1BQ4*2722-,3775T&温5-37.M15-37J29O-39.2333-36l43W23-22.呢®5-21L9022-22.1901-21.49

20、90-22.&407-21.41B1-20.4147-21.7874-21L5722-22.S290-223085图5.1主成分分析结果如图可见是一个1000组新的变量，由于数据集比较大，在这里只截出一部分下面再对这1000组新变量做聚类分析处理。5.2.3 聚类分析(1)聚类分析的基本思想：聚类(clustering),简单的讲就是将一个给定的数据集分成若干个不同簇的过程。聚类算法中的簇指的是数据对象的集合，且这种数据对象集合必须满足条件：同一簇中的数据对象间具有较大的相似性，而不同簇中的数据对象间具有较小的相似性。聚类的主要指导思想就是尽可能使同一簇内对象相似度达到最大，且不同簇间

21、对象相异度达到最大。（2） K-means算法：word完美格式专业资料首先从含有n个数据对象的数据集中随机选择K个数据对象作为初始中心然后计算每个数据对象到各中心的距离，根据最近邻原则，所有数据对象将会被划分到离它最近的那个中心所代表的簇中，接着分别计算新生成的各簇中数据对象的均值作为各簇新的中心，比较新的中心和上一次得到的中心，如果新的中心没有发生变化，则算法收敛，输出结果，如果新的中心和上一次的中心相比发生变化，则要根据新的中心对所有数据对象重新进行划分。直到满足算法的收敛条件为止。（3） K-means算法的实现:从含有1000个数据对象的基因表达谱（附录一）中随机选择100个数据对象

22、作为初始的聚类中心；两基因数据的相似度可通过计算两个基因数据的欧式距离来得到，再根据最近邻原则将数据对象逐个划分到离其最近的聚类中心所代表的簇中，计算误差平方和准则函数E的值；更新聚类中心，即分别计算各个簇中所有数据对象的均值作为各个簇的新的中心，以新的聚类中心来计算误差平方和准则函数E的值；将步骤计算得到的E值和前一次计算得到的E值进行比较，若两者差值的绝对值小于等于预先设定的阈值，即聚类准则函数收敛，则转步骤（5）,否则转步骤；输出K个聚类。5.2.4数据降维后的结果显示将10000个基因数据通过主成分分析和聚类分析缩小到30组特征数据。图5.2各类中基因数目分布word完美格式专业资料5

23、.3 模型一的分析目前重构基因调控网络的模型主要有贝叶斯网络，线性组合模型等，本题我们采用的是线性回归模型，下面先对基因调控网络与线性回归模型的对应关系作分析，以证明模型的可行性：假设基因Ak受基因AK1,Ak2,.Akm调控，则认为相应的基因表达值，有以下线性关系：at,k一blat-1,k1b2at-1,k2.bmat-1,kmk式中at,k代表基因Ak在t时刻的表达水平，bv.bm为常数。若基因A,B,C之间的真实调控关系如图所示，其中A,B,C代表基因，而边代表调控关系，比如从A到B有一条有向边，代表了基因A对基因B有调控关系。其对应的线性模型为:at,c=b1at-1,ab2at-1

24、,b1at,b=b3at-1,a25.4 模型一的建立5.4.1 模型说明：采用多元线性回归模型设因变量y与自变量x1,x2，xp,的线性回归模型为y="1+吆2+.+，p+&（2-1）其中，1,久，.，灯是p个未知参数，称为回归系数。y称为因变量（被解释变量），而X1,X2,.,xp是p歌可以精确测量并可以控制的一般变量，称为自变量（解释变量）。上面的公式就称为多元线性回归模型。名为随机误差项,对随机误差项我们常假定word完美格式专业资料'E=02(2-2)var(s)=a称E(y)=PiXi+P2X2+.+PpXp+名为理论回归方程。如果有n个样本数据知n,&#

25、39;'）,i=i,2,.,n,则样本数据代入（2-3）式可表示为：y=4再i+四再2+吃邛+ay2=Ax2i+Pixn+Bp/+Jyn=+幺相2+用/即+4写成矩阵形式为:5.4.2 本题中的基因调控网络与线性回归模型的对应关系：17位志愿者，每位志愿者体内的P个基因数据，这里P=11961。根据时间序列排布，一共16个连续采样时间点，获得272组观测数据。其中有四组错误数据，排除后剩余268组数据xt1，xt2,.,xtp,这里t=1,2,.,268,其中，xti,表示t时刻时基因i（i代表一个基因）的表达值。对于基因i,i=1,2,.,p,若假设它与全部基因（包括它本身在内）存在

26、线性关系，则它的线性回归模型为4=AlVu+0Mq+乩为+岛（3-5）代入268组样本数据，得：X2i丹丙I+再2T卜PipXp+%；X3i=Al-t21+_*fiipX2p+62XTi+Pi2XT-2+%写成矩阵的形式表示为：Z=仇风+星与+6word完美格式专业资料也即式中,V二xi月=>fill，,Pip).与)矩阵X是一个（T-1）*P矩阵，称X为回归设计矩阵或者设计矩阵。对于所有P个基因的全模型，可以表示成统一的矩阵形式如下Y=XB+E其中，Y=yi>y2）>yp>B=%，02,，口p'，E=、1,82,，即。在处理基因调控网络重构问题时，实际上是逐一

27、寻找对基因i,i=1,2,p起调控作用的调控基因的集合，而不是从基因i,i=1,2,p调控哪些其他基因的角度出发。这从网络的拓扑结构上来讲，我们是在逐个寻找网络中每个顶点（基因）的入边的集合（即其他基因对该顶点对应基因的调控关系），找到了每个顶点的入边集合后，由顶点和边组成的整个网络的拓扑结构就确定了。以上线性回归算法用Matlab实现，重构出的基因调控网络如下所示：word完美格式专业资料图5.3基因网络重构效果图6.问题二的建模与算法实现求解6.1 模型选择问题二：在已经重构好的基因网络的基础上寻找导致病毒感染人体以后导致支援者是否产生严重临床症状的蛋白质。假设是否展现严重的临床症状为最后

28、的因变量，那么由于这个因变量是二分类，所以考虑用逻辑回归模型去对基因对是否产生症状的关系做一个分析。通过逻辑回归我们可以得到一个基因对于严重临床症状的相关性的大小，以及对产生临床症状是起着正向的作用还是抑制的作用。我们要对所得的数据进行分析，寻找与染病相关系数大的基因，然后去依据基因表达谱中的数据中对于基因的描述去进一步确定关键蛋白质。6.2 模型简述该模型的对输入没有严格的要求，而且模型简单直观，又容易解释.而且该模型并不容易产生过拟合。而当前题目中只有产生和没有产生临床症状两种不同的结果，所以采用此模型去寻求重要的蛋白质是一个可以考虑的方法。定义一种概率函数1一、昨Prl.XifXzF.X

29、n=Xn)=占Xi_4Xn(6-1)1.a*e1f对log()进行线性回归。对于上述模型可以米用最大似然估计方法对冗其回归参数进行估计。既用总体的分布密度或概率分布的表达式及样本所提供的信息求位置参数估计两的一种方法。1.3 模型建立注射后不同人随时问是否产生症状不同。回归的因变量是二分类。因此采用的是word完美格式专业资料LR分类器的方法去求解。将每一个志愿者的每一个时间节点作为一组数据，而每一个基因的值就是该组数据的一个属性。如果一个志愿者在某时间点产生了严重的临床症状，则因变量结果为1,相反没有产生严重的临床症状则为00本次实验有17个志愿者，每个志愿者有16个时间节点的采样数据。又由

30、于有四个数据有较大误差，故对268组数据进行LR分类。得到每一个基因与最终是否产生临床症状这个因变量的相关性。通过相关性的分析便可以得到能够推动产生临床症状的关键基因，从而推出关键蛋白。将通过逻辑回归所得到的系数矩阵绘制如图6.1各基因与临床症状的相关性1.4 模型估计结果设定一个阙值，将大于阙值的基因数据拿出来，对所拿出来的基因数据进行分析。满足要求的基因在基因表达谱中的编号如下：图6.2基因在基因谱中的位置对应基因表达谱，将上述位置的基因编号拿出，方便后续数据分析，如下:个教1指针编号10675_at21742_at33669_at464135_at527074_at图6.3基因以及基因的

31、指针编号对照基因表达谱中的数据得到了四个对于最终的因变量影响较大的基因，在谱中寻找所对应的蛋白质的信息如下表：word完美格式专业资料Gene_IDSymbolDescription1742DLG4discs,largehomolog4(Drosophila)3669ISG20interferonstimulatedexonucleasegene20kDa27074LAMP3lysosomal-associatedmembraneprotein364135IFIH1interferoninducedwithhelicaseCdomain1表6.1对临床症状产生重要的蛋白质LocationTyp

32、e(s)PlasmakinaseMembraneNucleusenzymePlasmaMembraneNucleusenzyme综上，结果认为导致志愿者产生严重的临床症状的蛋白质主要有四种，其标志依次为DLG4,ISG20,LAMP3,IFIH17 .模型的分析及改进7.1 对线性回归模型拟合优度的分析假设Ki为第I个类用线性回归重构基因网络的拟合优度，那么该网络的平均拟合优度i=mK(Ki/m)a(7-1)带入数据计算得到，平均拟合优度大约在0.8左右，所以该网络具有一定的误差，但是大致对于一个基因网络还是有一定真是的反应。7.2 改进在使用线性回归模型区重构基因网络的时候，虽然不容易产生拟

33、合过度的情况，但是线性回归和逻辑回归都只能对变量只见一个线性的关系具有良好的刻画。对于完全不相关的两组变量的时候数据的误差就会增大。所以可结合非线性动态贝叶斯网络来重构基因网络，并且在逻辑回归之前尽早的取出不相关变量，对变量进行分析。word完美格式专业资料8 .附录1287*6-13-1.403le-16-I.7525e-157.0617f-16T对913-5.1264dS-l61JOOOO1.Sl4e-14Kl602e-l5也乃7&S-16-3423fte-l?2.4149e-l3-1-51Me-16-3.429Ze-l62.591&e-169.8l?3e-1?MaF1B1

34、母3W156.7713e-1fi1.3463e-Hfi454&9Z5-Z2&17e-15L6J50e15-I.5437e-153.7135e-H60.5441-0134157-D.5M30.D0434K口112g-O,113ldQOM-D.03D40.D256-0,S24B-Q.1D430.0077F|-1-M72e-1b|-15951>-1162s9901e-16-4.3BWe-161.304.gff-165u73HTe-163J1SDe-16-147git-16-1.5166e-15Z2gflOe-15-3t6571e-16-3.9TD3e-llG-23166-16-

35、0.D61D-033440.0654.D.332O-0.M410,D124OJDfljS-D.OOB5DU0772-0.D2450.21159OJD064).007o.D&ai4).0645-D.3272-D.4819-O.D1I22-O.1DTQ0.00&4-D.O761I口加税-0.2127-Q.0M10.00721J257e-16-1.1W1r-166.5529&-16lJtM7De-16-2.14B7e-17-9.5415c-167452Be-172.3W7e-16BL6442e-1fi-7.K77e-1l71.0000-0.03004K1&250JQ3

36、22DL1IB31口.2如5，口皿口0.DD610.0042口郁而-0.01210.10620.04450.003*-1.562(k-169019m-5.079417，玷9河,171-155通-2.21We-174>M2Je-1Ei-6LQ769e-17*强4.1工,丸1英Z"了1N0CM152Tz257*16iJOOW701-4.4645*-16区整35力173341*1$4.3W6t.1715圣55.155.5071t.15M5M1512%U+，15和15125fl»16544fl1»-175,903gli乙433”7-1e725&fd61.6|

37、42tie5377741752-162435-164潮1,口姬4157.186*-17徇吐1$SUM*帕瑞卸*帼-fi,0071e-174Ji»Dv-15i1.15047#.153而ME-S.SOO-172J>134e-16-5.187&B-174,0764e-i6i37l2*-164,1717e-17-2.60318-165.0964e-i*9.9122e-171ra93«d54.7323a-16-1,i752fl-16<5BO6e-17-0.0319-0.17270,173502556-0.0a2由05蚓-0XKH5O.OSOJ-0.012SA112

38、S曲05施0.0590为-40613-0.3210-0.472&O.OCMO-0.0120910删OXIOS2460.037-0.0973-£.1l03e-16-.?S33r-l6TJNMrSlJa5Te-18S47S7e-14露1惘M15-1.00T5e-14-4lDffi5e-151星蜘-13如如"76-lJ733e-13-0.01134-0.0727ftO144口和310.107&-ai115e-M002可0.0239-DJ00190,D17D-0.00541QM75aozsi4口居N-Q.1F9-7D.O356DL19D60.366D-0.0023O

39、.OD&rcm湖-D.OCM*0.M2D-D.D1I34I0.11740.05473-M11e-16<1.5641e-16-4.E313e-17-2.2501e-1746775e-161.05DBe-16-2.73i8De-1644256fr.l币1.9351?-13.792fie-15-Z2a79e-15=ZH72e-H66.2434e-16>图8.1基因网络系数矩阵1132.0194e-16-5.ia7fe-17-1.S&2De-164.0764e-162.3712e-164.1717e-17-2.®51=-1651如9.51Z2e-17-3.ai53

40、e-1l5-1.7323C-16£1.1752*-164.58D6e-1714-OJQ319-0.172T0.0342CL17J50.2556-D.(M22Q.DM5DUO&&SQDH5D.04D3-CI.D12&0.1123口.05州150.0WO03195/D6居-D.321O-0.729DUD04Q-0.D12D-D.1DW0.DOSS-DUO0.D237小金四-0.097316-2.11D5e-16-1.5945ff-174.1395c-16-9L7539e-16-1.7222e-161.5851-1854757e-161.1D65t-15-1.(&#

41、187;75e-16-41D05e»151.5434e-155.6M2e-16-1月75翡书174U0134-Q.D727(k口1440jD731O.1D76-5.111l5e4MO.DO27ugg也网gD.0170-0.D054<1.04750.022116心-0.1797a0256D.13D60.2&50-0.00230.DD670.0591ASMSD.0420-0.D124(11174D.0547193.S611t-161.5641e-1fi4,B31Je-17-22M1e-174.6775*-161.a&DSc-1fi-2.714-IE-4-525&am

42、p;e-161.913511£-1&J.7926e-1l5-22£l79e-15-2.7172e-1620DJ023SQ.12e64J.D255-0.1293-O.1M4D.0O16-0.DD4&-DLM23口侬m-003010.D096-49.0640-0.039221OOO57O.D309ftDOfil-Qjoaw-O.D4573666e-D4-0.D012-DU01017.9'5e-D4-DJ0072O.DD23-OJ02O2-O.D09422OJOZ27.1212-0.0244-0.1I23A-Q1£220.00154HME-&am

43、p;JM05口侬2-4UI2BS0.KH911-aoas-0.0175221J0D002.141Be-16g回血-16工蚓了k173.5Q74e-17的比e-17-4.4716e-1i54.5d0c-171U175&-1152.01&4e-16Q.omE1.5621e-16240JQ216.11ea/*34<11as川751CL00154tMM4DlO姻0.D031o.DaeaLfl773-0.03602$4lugi72-0M34队口侬口加邦.1胆-DLQ012*00350.DO24琬阳自,口.DoegOJD610口.a2a4263.416e-1G1L的95酰15.1.2

44、24151.7743S-151-U7W-1&可鼻1*15OMSe-ISn-1.051410745b-151庭吠152700210Q.1139ME-anus1惭(MKI14由皿3他次0.5-00744=0.1347演00714e.丽5喀ng增-CL12720,0100owe7心电-Q.117B2906510揶取郭耶-MOKQ.Q244M141Mice我国104$4口惟汕-4S7$9r-1§1.W0016构7E64.1册加1$区1咖e16-2.S564«d6壮笈打.博Memo-isS10*155.7碑0口155,0191<-166J712赛图8.2基因网络系数矩阵

45、211.7490e-161.9364c-15dJ142td521期*15-4J09fl2td62.Z27&-15r1_22Me-15-1.22触T5心511i.saaieiGi13W5吐博4.BniOe-15-3.M51tdl52I.OODO1.2484«-15-4.4fi11e-161.242ie-151J8227&-16-1424De-16-1-5259&-16g.9553e-173.1120-16-54361e-163j62O1t-161.34gfie-15-Z.a3O6e-151D.00M9103口HOST-0.1721D.(E4fi-aocs口-DJ

46、0634-O.D197D2247-CJ.JM9-DL15404-瓦21gbi了-7,B56Be-16L9847ff-16-9.032&e-16M910ff-171.B434e-151,0000-12226e-16-1.4049¥-16-2.2321e-H5-2,DCWa?-1B1.20&6e-1550.0055注脸T-O.OOSJ0,1*37OM幻也1的M-O052J0,00770m的-0215S0,M29*丛14西$4855-ft39670.<J05谈M5MOOSV-0,0076-M«0-G.0186CL2126©2141*577>2

47、l1O5td744158,16用湖青17，1.453,16WglT.IT.5424edg方4制35小1已修21T1了£用99w174.PQKM-1&4.17Q9f150.0027o.mo-0.00230.0510-0LQ2?g0.DB13：口力2弼O.D03BD-03DQ0.D041-H1D610.144D0U&72792.7i373e-17-1.D45&e-16El27B3e-185.IM5&e-22-4.9M1&-17-5.1459C-166-8440e-17154i-171.3756e-17-1.4393e-16-7.56OOb-177.

48、5O1&-17&.6D46e-17ID-MS13e-17-2.&&5&e-1S126B7ff-16-I.WSIe-151J6293EE-1,S9OTC-15-4.2636e-164.D999«-166.2272e-16-KaWSe-TS-1.2003e-155.1767E-16in-2,lOT4e-l61XBM-1,&03&e-16-l.lSftle-lS-2.6303e-174T皿H工辅浪骨-±1兆於-MS.2:M-161,0246-15i21必6*-t柝1,ifi79e-15曲W-4,27*3*-16i4MW-1.

49、762&B-163.i3J4e1754453#14i-40M*-i51.iOO2e-l5£72第e-16135l535>3«-17电.7892316丸7aT5£.4Wr1615,we175.5W9t.171.CW«闻61263.5«9»-16旧知鹏博14D.00H0.2104-QjOOSa0.0542-DLO297Q.0864-CL0274Q.D040DU&31S口.0099-DL112SQ.1531DU077315O.(M54-0.3*95DjOC55-O.1D03DJSE49-0.1598DJ05D7-0.D

50、D74-D-05M-o.DiaiDJDB6-a.MSI-D.14I3O16Z19B4?-172.25744：-16-2L82B4?dfi1.6252C-16-2l17ige-1T1.&637e-15ZJiWOe-166.315&C-171.00008J26fe-17-13701e-15-1,DOO3e-151t6492e-1517O.OTIr也晒6-0X)0130.022B-XH2501.0364WlS010017CU>1M0.0W2-DLMI750.0645MS25IB0.2190-0.003.li0.0564-0X>3090.帕躺-aoaflrsao(M2O4&

51、gt;331MWS-0L11?40J59aCuOS(M195.”底164.24酶.像4目5*16-3,674161-OCXX)am通-17山&41.162.6951777本呜3.K0i*d6.5.041Of16图8.3基因网络系数矩阵3word完美格式专业资料135.5355e-17国7蚪5-1.1031e-17-1.3391-1日3.7523e-16-2.4648e-16-Z.42B2e-163.9634可5.59a9e-171.00003.5Q09e-HfiJ3.2201e-16)140.0025DLZ104-0.D03D0.05420.0364Q.D2740.0D40oojiaD.009S-0.112B00.0054-0.13910.D0550.1DMOJ0644-0L