关于疾病研究问题的数学模型

《关于疾病研究问题的数学模型》由会员分享，可在线阅读，更多相关《关于疾病研究问题的数学模型（23页珍藏版）》请在装配图网上搜索。

1、关于疾病研究问题的数学模型摘 要现在与50年前相比，生出一个活泼健康的孩子越来像一个复杂的系统工程，需要精心运作每一个细节，到底是什么原因威胁胎儿的健康也是医学上一个非常复杂的问题。本模型就题目给出的鼻炎家族史、主（被）动吸烟等12个因素的相关数据，建立多分类条件Logistic回归分析数学模型，利用SPSS软件包拟合出Sig-显著性（P）的值，再利用逐步回归中的后退法，把所有的变量引入模型，然后逐次把Sig值最大并且Siga（a=0.05）的变量剔除出模型，每次只剔除一个变量，每次剔除一个变量后重新拟合模型，按照上述剔除标准继续剔除变量，直至模型中的所有变量的Sig0.05为止。最后所剩Si

2、g值0.05所对应的因素:鼻炎家族史，孕期用药，先兆早产，剖宫产和孕期营养即为诱发胎儿疾病的主要影响因素。关键词：疾病 Logistic回归分析 SPSS软件包 Sig(P)值 逐步回归后退法一、 问题重述世界卫生组织最近发布的资料显示，与50年前相比，人们的生殖能力明显降低了，生一个活泼健康的孩子，越来越像一个复杂的系统工程，需要精心运作每一个细节，不敢稍有松懈麻痹，到底是什么原因威胁胎儿的健康也是医学上一个非常复杂的问题。本题就胎儿的某种疾病影响因素做了详细的调查与记录见附表1，根据附表1中的数据，利用统计学知识判断到底该疾病的诱发因素有哪些？二、 模型分析考虑到题目给出了病例组与对照组，

3、且给出的多个影响因素为分类变量，不能使用线性回归模型，因此，考虑使用Logistic回归分析数学模型进行了统计分析。在考虑用逐步后退法分析出该疾病的诱发因素。Logistic回归分析数学模型是专门用来分析疾病与危险因素间联系的一种统计方法。假设影响因素对疾病的影响是独立的，利用SPSS软件包输入数据，拟合出Sig值，Sig值为结果可信程度的一个递减指标（即将观察结果认为有效即具有总体代表性的犯错概率）。本模型中取a=0.05（在许多研究领域，0.05的Sig值通常被认为是可接受错误的边界水平）。当Siga时，就认为此Sig值对应的因素即为该疾病的诱发因素。三、符号说明y ：患病情况（0=未患病

4、 1=患病）x1 ：鼻炎家族史(0=无 1=有)x2 ：主(被)动吸烟(0=无 1=有)x3 ：孕早期病毒感染(0=无 1=有)x4 ：孕期用药(0=无 1=有)x5 ：先兆早产(0=无 1=有)x6 ：流产史(0=无 1=有)x7 ：早产(0=无 1=有)x8 ：剖宫产(0=无 1=有)x9 ：孕期营养（0=一般 1=好）x10：孕期活动量（0=小 1=一般）x11：出生体重（2.5=1 2.54=2 4=3）x12：妊娠反应（轻1 中2 重3）P：在m个自变量的作用下患病的概率 : 回归系数（j=0，112）Si：Sig的检验值Smax：Si的最大值S: Sig的临界值Sig=0.05Si

5、g： “显著性”，当0.01Sig0.05,则为差异显著，如果Sig0.01,则差异极显著。赞同 四、模型假设1、影响因素对疾病的影响是独立的。2、各个因素的有无影响程度大小不考虑。3、除这12种因素外有可能引起该病的因素不考虑。五、模型建立1、y=1 患病；y=0 未患病；自变量 在m个自变量的作用下患病的概率记作： （1）即： （2）将P作logit转换为，有 （3）2、建立似然函数模型： （4）通过，求出，代入（3）得logistic回归方程同时也可求出Sig值。m=12,j=(1,2200)3、逐步回归的后退法： 根据原始数据由SPSS软件得出多元线性回归结果； 找到最大的Sig值Sm

6、ax与临界值S=0.05比较，如果SmaxS，剔除Smax对应变量，重新根据余下的变量数据由SPSS软件再得出多元线性回归结果，再返回（1），否则进入（3）。 如果SmaxS，则逐步回归结束，即剩下的变量统计显著水平高。根据上面步骤，可以简单绘出流程如下：原始数据SmaxS剔除回归模型SPSS计算查找判断所剩因子六、模型求解：借助于SPSS统计软件完成具体的过程如下：Data View：（数据窗口）（见附表1） Variable View：（变量窗口）再如下：AnalyzeRegressionLinear再如下：Lineardependent：y independent：x1-x12然后单击O

7、K,得：以上的Coefficient，为回归参数表。再根据逐步回归的后退法，得：孕早期病毒感染的Sig值最大0.9050.05,剔除。孕期活动量Sig值最大为0.8580.05，剔除主（被）动吸烟Sig值最大为0.6160.05，剔除。早产的Sig值最大为0.3970.05，剔除。流产史的Sig值最大为0.1920.05，剔除妊娠反应的Sig值最大为0.1860.05，剔除出身体重Sig值最大为0.1800.05，剔除所有变量的Sig值0.05，故最终回归模型Logit（P）=0.111+0.144x1+0.184x4+0.400x5+0.186x8+0.188x9综上知：:鼻炎家族史，孕期用

8、药，先兆早产，剖宫产和孕期营养为诱发胎儿疾病的主要影响因素。七、模型检验方法1：根据t的检验方法，试验证，如下：用SPSS软件得：由于Sig的值都小于0.05，从而验证正确。方法2：判别分析法与回归分析的最优判别函数是一致的。用Fisher判别法的判别函数经计算得：y=0.0051x1+0.0009x2+0.0002x3+0.0052x4+0.0116x5+0.0025x6+0.0033x7+0.0051x8+0.0047x9-0.0004x10-0.0035x11+0.0016x12判别系数大小前五的为：x1 x4 x5 x8 x9与上模型计算结果一致。八、模型应用可以用来预测在不同的自变量

9、情况下，发生某病或某种情况的主要影响因素。适用于与本题类似的医学研究问题，也适用于其他情况。在生物学中生物呈s增长的相关研究问题也适用。如基于Logistic回归服装品牌忠诚度影响因素分析等。 九、模型评价在影响结果的混杂因素比较多的时候，难以满足均衡可比的要求，且像本题给出各因素的定性数据，不能直接利用线性回归进行分析，因而Logistic回归模型更为适合。利用SPSS软件包计算方便、精准，易于判断、理解，但是，适用本模型的前提是自变量对疾病的影响是独立的，实际情况有时并不如此，就需利用其他方法进行降维处理，再利用本模型进行统计分析。由这种模型，根据Sig的值，也可得到最重要的影响因素，依次

10、次要的因素等。如本题：先兆早产Sig值最小，为最重要的，次要的为孕期用药，剖宫产和孕期营养，再为鼻炎家族史十、参考文献十一、附录附表1：病例组：患儿代码鼻炎家族史主(被)动吸烟 孕早期病毒感染孕期用药先兆早产流产史早产剖宫产孕期营养孕期活动量出生体重 （kg）妊娠反应无=0 有=1无=0 有=1无=0 有=1无=0 有=1无=0有=1无=0有=1无=0有=1无=0有=1一般0好1小0 一般12.5=1 2.54=2 4=3轻1 中2 重31111101010121200000101112230110000101214011100000121501001101112161101010110237

11、010001000123801011011102290100000001211000011000102111100010011021120011110110221311000100013114100110011023151001000011211611101001102217011100000122180001000001211901110000113320000011000121210101100111212210000001012123000000010123241001110111212501111001002226001101010123271001000011212810110001

12、112229000110101021301011100110213101100101112332111101011021331101001111213411010011112135000010011013361011010001213701000001012138010010011121391111000000234000000011112141100110000132421101010100234311111101012144011010011023451011000010214600010111001247101111010031481011100111214900000100002350

13、010101011121511001000101225200110110102353111010010121541101000100225500010001112156000010001121570001000101225811011001112159110000000021600011010101316100111101112362001001011121630100101000216401001001112165000100010122661101100010226711000101112168110100011123690110110011217001010000112171011010

14、011121720101000011237311000100112174011100001121751100100100217600110101012277010000001121781100111110217910110001012280010011001021811110000110228200010000012183010010011123841100010111218500110100002186010111011123871111000111218811110000012189010001011121901010010100219110010001012192000010001123

15、93001101001123941111100001219511110000012196010000010123971010000110219810010101112199000001011133100001101001121对照组：儿童代码鼻炎家族史主(被)动吸烟 孕早期病毒感染孕期用药先兆早产流产史早产剖宫产孕期营养孕期活动量出生体重 （kg）妊娠反应无=0 有=1无=0 有=1无=0 有=1无=0 有=1无=0 有=1无=0 有=1无=0 有=1无=0 有=1一般0 好1小0 一般12.5=1 2.54=2 4=3轻1 中2 重310100000101212001100001121300

16、111010002341000010011325101100010121601000100013170000000011218110100000121901000001012110011001000121111010000101211201000000012213011001010021140111000001211500000000112116010100000122171001010110311810010000012219000000011121201000000001312101100101113222010101001123230000000001222400000101012225

17、010100000121260101000101112701000000012228011000000121290100000001213000000100102131010100010121320000000001213301010000112234011001000131351000000101213601100000002137000100010121380100000001223901100100111140000000000121410100000001224211101001012343000001001021440101000001224500000100013146101000

18、100121470101000101224800100100012149100100000133500100000101215110000000012152011000001122530001010001315411000001012155011101000121560000000010225701010000012158101000010131590110010011226000010000012161101100011123620100000001216301000000112264101101010121650100010101216601010000002167000000001121

19、6810010100011269000000010121700111010110217100000000013272110100000121730110010011237400000001013275000100000022760110010111217701100000013178100000001122790001000011228001110101013181000000001121820100010110228301010000012184011100011122851100010001218600000001102187001100000121881100010101328901000000012190011000000121910000010111219211010000013193000000100123940100000101219500111100102196110000000121970101000101319800100000112299010001000121100101100000121附表2：