疾病研究问题.docVIP

影响胎儿健康的因素分析 摘要：本文利用主成分分析法建立了影响胎儿健康的因素的数学模型，并结合实际提出了通用可行的算法并利用MATLAB软件进行实践求解。对于这个问题首先把所给数据用的矩阵表示，通过计算得到矩阵中每个因素的平均值与方差，借用平均值与方差分别表示数学期望与方差的近似值，将原先的矩阵标准化，再求出这个矩阵的相关系数矩阵，求出它的特征值、相应特征向量以及贡献率和累计贡献率，累计贡献率大于等于85%的值有9个，这样就将本题中原有12个影响因素减少为9个，分别为：鼻炎家族史、主(被)动吸烟、孕早期病毒感染、孕期用药、先兆早产、流产史、早产、剖宫产、孕期营养，通过特征值及主成分贡献率的图表可知，前9个特征值的的累积贡献率已达86.52%，说明前9个影响因素已经基本包含了全部影响因素的信息，所以，只需取前9个影响因素进行计算即可。最后利用主成分值比较分析得出影响胎儿健康的因素有六个，其中包括鼻咽家族史，孕期用药，先兆早产，流产史，剖宫产，孕期营养。本文最后还对模型进行了推广并分析了模型的优缺点。 关键字：主成分分析法 MATLAB软件 累计贡献率  1 问题重述 世界卫生组织最近发布的资料显示，与50年前相比，人们的生殖能力明显降低了，生一个活泼健康的孩子，越来越像一个复杂的系统工程，需要精心运作每一个细节，不敢稍有松懈麻痹，到底是什么原因威胁胎儿的健康也是医学上一个非常复杂的问题。本题就胎儿的某种疾病影响因素做了详细的调查与记录（见附录表一、二）。 问题：根据附录表一、二中的数据，利用统计学知识判断到底该疾病的诱发因素有哪些？ 2 问题分析 本题是分析关于影响胎儿健康的因素的问题，题中给出了12种因素，要判断这12种因素与胎儿的健康是否有关，既要排除无关的因素，也要找出有显著影响的因素，题中给了两组数据即病例组和对照组，显然，我们要以对照组为基准去分析病例组的异常情况，如果用统计的方法研究这个多变量的问题，量太大，会增加计算量和分析问题的复杂性，因为主成分分析法可以将变量减少，因此，我们将利用主成分分析法，减少影响胎儿健康的因素，得到诱发胎儿致病的因素。 3 符号说明 Xi 影响因素 第i个影响因素对应的第j个患者 每个因素Xi的平均值 每个因素Xi的平方差 第i个特征值 R 相关系数矩阵 Fi 第i个主成分 4模型假设 1 假设各数据有效完整 2 假设影响胎儿健康的因素只有：鼻炎家族史、主(被)动吸烟、孕早期病毒感染、孕期用药、先兆早产、流产史、早产、剖宫产、孕期营养、孕期活动量、出生体重 （kg）、妊娠反应这12个因素 3 各因素之间不相互作用 5模型建立与求解 表一中的病例组共有个数，每个因素Xi的平均值与方差按下述公式计算： 经计算得： Si为标准差，如果把Zi作为随机变量，则常用，分别表示Zi的数学期望与方差的近似值。 2.为了方便计算，我们将表一中的病例组数据用一个矩阵来表示： 即 即 将X通过以下的标准化变换为矩阵Z的过程，称为标准化变换。 系数相关矩阵为： 即 经MATLAB软件计算得出相关系数矩阵R对应的特征值及相应的特征向量，并把特征值从大到小排列，见表一: 表一 特征值 特征向量 1.7453 （0.3943，-0.1083，0.1528，-0.2062，0.2251，0.0633，0.2405，-0.2014，-0.0617，0.3883，-0.1347，0.278） 1.3926 （-0.0127，0.4289，-0.6098，0.8224，0.7306，0.7078，-0.3049，-0.0811，-0.4734，0.2354，-0.3979，-0.0452） 1.3171 （-0.3781，0.0655，0.2728，0.0877，0.2782，-0.1913，0.4393，0.4671，-0.4337，0.2323，0.2482，-0.0879） 1.2066 （0.1292，0.1886，-0.2154，0.1282 0.3517 -0.0011 0.1811 0.0088 0.0811 -0.0212 -0.0880 0.2570） 1.0408 0.3696 -0.1071 0.3176 0.1845 -0.1467 0.0884 0.1352 -0.0703 0.2959 0.1375 -0.1722 -0.0325 0.9669 -0.0358 -0.4371 -0.0891 0.7228 0.3651 -1