1、全国大学生数学建模竞赛C题国家奖一等奖优秀论文脑卒中发病环境因素分析及干预摘要本文主要讨论脑卒中发病环境因素分析及干预问题。根据题中所给出的数据，利用SPSS20 软件进行相关性统计分析，分别对各气象因素进行单因素分析，进而建立后退法线性回归分析模型，得到脑卒中与气压、气温、相对湿度之间的关系。同时在广泛收集各种资料并综合考虑环境因素，对脑卒中高危人群提出预警和干预的建议方案。首先，利用SPSS20软件,从患病人群的性别、年龄、职业进行统计分析，得到2007-2010年男性患病人数高于女性，且男性所占比例有逐年下降趋势，女性则有上升趋势，因此，性别比例呈减小趋势。分析不同年龄段患病人数，得到患
2、病高峰期为75-77岁之间，且青少年比例逐年呈增长趋势，可见患病比例趋于年轻化。同时在不同的职业中，农民发病人数最多，教师，渔民，医务人员，职工，离退人员的发病人数较少。其次，由题中所给数据先进行单因素分析，剔除对脑卒中影响不显著的因素，得出气温、气压、相对湿度对脑卒中的影响程度大小，进而采用后退法线性回归分析建立模型，利用SPSS20对数据进行分析，求得脑卒中发病率与气温、气压、相对湿度之间的关系。即发病率与平均温度成正相关，与最高温度成负相关，发病率与平均气压成正相关，与最低气压成负相关，与平均相对湿度成负相关，与最小相对湿度成正相关。最后，通过查找资料发现，影响脑卒中的因素有两类，一类是
3、不可干预因素，如年龄、性别、家族史，另一类是可干预因素，如高血压、高血脂、糖尿病、肥胖、抽烟、酗酒等因素。分析这些因素，建立双变量因素分析模型，并结合问题1和问题2，对高危人群提出预警和干预的建议方案。关键词 脑卒中 单因素分析 后退法线性回归分析 双变量因素分析一 问题的重述脑卒中（俗称脑中风）是目前威胁人类生命的严重疾病之一，它的发生是一个漫长的过程，一旦得病就很难逆转。这种疾病的诱发已经被证实与环境因素，包括气温、湿度之间存在密切的关系。对脑卒中的发病环境因素进行分析，其目的是为了进行疾病的风险评估，对脑卒中高危人群能够及时采取干预措施，也让尚未得病的健康人，或者亚健康人了解自己得脑卒中
4、风险程度，进行自我保护。同时，通过数据模型的建立，掌握疾病发病率的规律，对于卫生行政部门和医疗机构合理调配医务力量、改善就诊治疗环境、配置床位和医疗药物等都具有实际的指导意义。数据（见Appendix-C1）来源于中国某城市各家医院2007年1月至2010年12月的脑卒中发病病例信息以及相应期间当地的逐日气象资料（Appendix-C2）。请你们根据题目提供的数据，回答以下问题：1根据病人基本信息，对发病人群进行统计描述。2建立数学模型研究脑卒中发病率与气温、气压、相对湿度间的关系。3查阅和搜集文献中有关脑卒中高危人群的重要特征和关键指标，结合1、2中所得结论，对高危人群提出预警和干预的建议方
5、案。二 问题分析2.1 问题的背景随着城市进程的加快、人口密度的加大、生活习性的改变、全球气候的变暖，一些严重威胁人们健康的疾病，如脑卒中的发病率呈现上升趋势。全球每年有460万人死于脑卒中，中国每年死于脑卒中的病人有160多万。脑卒中以其发病率高、死亡率高、致残率高的特点严重影响了人们的正常生活，成为仅次于癌症的“第二号杀手”。脑卒中一旦发生，多数患者治疗效果不能满意，完全恢复正常者只占少数，大多数患者会遗留严重的后遗症，给患者和家庭带来沉重的心理和经济负担。其成因包括环境、情绪、生活习惯、生理状况等多种因素，气象条件是脑卒中诱发的主要因素之一，根据气象条件变化对脑卒中发生的可能性进行预测，
6、可以指导人们在不同天气条件下对脑卒中采取不同的预防措施，从而预防、减少脑卒中的发生。如果能对脑卒中高危人群高位程度、患脑卒中的类型、患病风险、以及患病时间进行准确的评估，从而提出有针对性的防范措施，可以帮助医护人员对脑卒中的诊断，对于卫生行政部门和医疗机构合理调配医务力量、改善就诊治疗环境、配置床位和医疗药物等都具有实际的指导意义。2.2 对发病人群的统计描述根据题中所给出的数据，利用SPSS20 软件，分别从性别、年龄、职业这几个因素出发，对数据进行统计性分析。男女之间不同的生理结构，抵御环境的能力的不同，对脑卒中的影响也不同，分析男女在2007-2010年间患病的比例，进而可以更好的预防脑
7、卒中。随着的年龄的增长，患病人数在逐渐增长，分析不同年龄段患病人数的比例，得到脑卒中达到高峰时的年龄段。不同的职业工作的环境不同，受教育程度不同，个人习惯也不同，进而影响到其得病的概率。其次，不同的职业对脑卒中的预防与重视度不同。农民可能由于经济状况的原因，不及时治疗，使得患病人数所占比例较大，而其他职业则相对较小。2.3 脑卒中发病率与气温、气压、相对湿度间的关系 由于平均温度、最高温度、最低温度、平均气压、最低气压、最高气压、平均相对湿度、最小相对湿度对脑卒中的影响并不都非常显著，所以本文先通过单因素分析得出最低气温、最高气压对脑卒中的影响不显著。然后建立后退法线性回归分析模型，利用SPS
8、S20软件得出脑卒中发病率与气温、气压、相对湿度之间的关系。即发病率与平均温度成相关，与最高温度成负相关，发病率与平均气压成正相关，与最低气压成负相关，与平均相对湿度成负相关，与最小相对湿度成正相关2.4 高危人群提出预警和干预的建议方案广泛收集各种资料并综合考虑环境因素，发现影响脑卒中的因素有两类，一类是不可干预因素，如年龄、性别；另一类是可干预因素，如高血压、糖尿病、血脂异常、肥胖、饮酒等因素。分别分析这些因素对脑卒中的影响，建立双变量分析模型，进而提出高危人群预警和干预的建议方案。三 问题的假设1.假设不考虑该地区人数的变化。2.假设不考虑脑卒中的患病类型。3.假设影响脑卒中发病率的各因
9、素之间互不相关。4.假设不考虑影响遗传因素、个体差异对脑卒中的影响。四 符号说明平均温度最高温度最低温度平均气压最高气压最低气压平均相对湿度最小相对湿度统计量随机误差脑卒中患病率五 模型的建立与求解5.1 模型一的建立与求解随着社会的发展，统计学在医疗气象研究中的作用越来越重要，成为一种普遍应用的方法，并为人们所认可，发挥着极其重要的指导性作用。本文运用SPSS 20 软件分别从性别、年龄、职业对发病人群进行统计分析。表1 脑卒中不同性别逐年发病情况年份性别频率百分比有效百分比累积百分比2007男729355.255.255.2女592844.844.8100合计132211001002008
10、男1038054.554.654.6女864645.445.4100合计1902699.91002009男5171525252女47804848100合计99511001002010男1053753.553.553.5女916346.546.5100合计19700100100图一 脑卒中不同性别发病人数比例趋势图从上表可以看出，2007年-2010年脑卒中发病人群中，男性高于女性。分析可能男性普遍有吸烟饮酒的不良生活习惯，还与当前男性可能比女性面临更多的生活压力，导致高血压、糖尿病、高血脂等脑卒中危险因素的增加。同时从上图可以看出，男性所占比例有下降趋势，女性有上升趋势，因此，性别比例呈减小趋
11、势。表2 脑卒中不同职业逐年发病情况2007频率百分比有效百分比累积百分比有效农民716454.268.668.6工人8146.27.876.4退休人员190514.418.294.7教师380.30.495渔民180.10.295.2医务人员230.20.295.4职工620.50.696离退人员4163.14100合计1044079100缺失278121合计132211002008有效农民1008152.967.867.8工人12526.68.476.2退休人员250013.116.893教师630.30.493.5渔民380.20.393.7医务人员190.10.193.8职工1070.
12、60.794.6离退人员8094.25.4100合计1486978.1100缺失系统417221.9合计190411002009有效农民149915.142.442.4工人147914.941.984.3退休人员1061.1387.3教师710.7289.3渔民60.10.289.5医务人员350.4190.5职工3133.18.999.3离退人员230.20.7100合计353235.5100缺失系统641964.5合计99511002010有效农民1100055.872.172.1工人13026.68.580.7退休人员213110.81494.6教师440.20.394.9渔民40095
13、医务人员130.10.195职工2531.31.796.7离退人员5032.63.3100合计1525077.4100缺失系统445022.6合计19700100从上表可以看出：（1）农民的发病人数最多，其原因可能与农民高血压病例规则用药治疗所占的比率低，抽烟饮酒在人群中比例大，以及膳食结构不合理、摄入膳食结构不合理、摄入动物性优质蛋白少等因素有关。也可能由于经济状况，对脑卒中的不重视，不能及时治疗脑卒中。（2）退休人员发病人数较多，其原因可能是退休人员年龄较大，脑供血，脑供氧不足等原因引发脑卒中。（3）教师，渔民，医务人员，职工，离退人员的发病人数较少。其原因可能与受教育程度，工作环境等因素
14、有关。表3 2007-2010年不同年龄段不同性别患病情况2007性别男女总人数年龄人数比例人数比例本年龄本性别本年龄本性别18岁以下1152.38%0.15%1047.62%0.17%211840岁11659.18%1.61%8040.82%1.36%1964180岁598457.52%82.90%442042.48%75.17%1040480岁以上110744.69%15.34%137055.31%23.30%2477总人数721855.11%100.00%588044.89%100.00%130982008性别男女总人数年龄人数比例人数比例本年龄本性别本年龄本性别18岁以下1934.55
15、%0.18%3665.45%0.42%551840岁21069.77%2.03%9130.23%1.05%3014180岁843456.94%81.35%637843.06%73.88%1481280岁以上170544.48%16.44%212855.52%24.65%3833总人数1036854.57%100.00%863345.43%100.00%190012009性别男女总人数年龄人数比例人数比例本年龄本性别本年龄本性别18岁以下4133.33%0.79%8266.67%1.72%1231840岁23460.62%4.53%15239.38%3.18%3864180岁419253.65%
16、81.07%362246.35%75.82%781480岁以上70443.32%13.61%92156.68%19.28%1625总人数517151.98%100.00%477748.02%100.00%99482010性别男女总人数年龄人数比例人数比例本年龄本性别本年龄本性别18岁以下1062.50%0.09%637.50%0.07%161840岁18562.93%1.76%10937.07%1.19%2944180岁850355.85%80.70%672344.15%73.37%1522680岁以上183944.16%17.45%232555.84%25.37%4164总人数1053753
17、.49%100.00%916346.51%100.00%19700由上表可知：（1）18岁以下的青少年患病比例逐年基本呈增长趋势，可见患病年龄比例趋于年轻化。（2）。患病人数主要集中在41-80岁，其原因可能是年龄逐渐增长，血液循环所需能量减少，代谢减慢。（3）80岁以上的患病人群中，女性高于男性图二2007年脑卒中不同年龄段发病人数趋势图图三2008年脑卒中不同年龄段发病人数趋势图图四2009年脑卒中不同年龄段发病人数趋势图图五2010年脑卒中不同年龄段发病人数趋势图从上图中可以看出，脑卒中发病高峰年龄段2007年为75-77岁之间，2008年为74-76岁之间，2009年76-78岁之间，
18、2010年为75-78岁之间。由此可知，脑卒中发病率高峰年龄段为75-77岁之间。5.2 模型二的建立与求解通过题中所给数据，首先通过单因素分析得出气温、气压、相对湿度对脑卒中发病情况的影响关系大小，其次建立多元线性回归分析模型，最后利用SPSS20软件对数据进行分析，得到气温、气压、相对湿度与发病情况的关系。5.2.1 单因素模型的建立假设a为给定的显著性水平，则列出：气温、气压、相对湿度分别对脑卒中的影响程度分析，即单因素分析，公式如下：若，表明均值之间的差异显著，该因素对观测值有显著影响。若,表明均值之间的差异不显著，该因素对观测值没有显著影响。5.2.2 模型的求解利用SPSS20 得
19、出平均温度、最高温度、平均气压、最低气压、平均相对湿度、最小相对湿度对脑卒中的影响显著，而最低气温、最高气压对脑卒中的影响不显著。5.2.3 多元化线性回归模型的建立在大多数的实际问题中，影响因变量的因素不止一个而是多个，因此对这类的问题本文采用多元化的线性回归分析。根据脑卒中发病率与气温、气压、相对湿度之间的关系建立多元化线性回归模型，如下：5.2.4 多元线性回归模型的基本假定假设1，解释变量是非随机的或固定的，且各x之间互不相关（无多重共线性）假设2，随机误差项具有零均值、同方差及不序列相关性，即： 假设3，解释变量与随机项不相关，即： 假设4，随机项满足正态分布，即：在医学病学研究中经
20、常会遇到某一疾病的发生与变化取决于几个影响因素的情况，也就是一个因变量和几个自变量有依存关系的情况，而且有时几个影响因素主次难以区分，或者有的因素虽属次要，但也不能略去其作用，这时采用一元回归分析预测法是难以奏效的，需要采用多元回归分析预测法。列出多元线性回归模型的一般表现形式为：将n个统计数据代入上述模型，则问题转化为：5.2.5 模型的求解利用SPSS20 软件对数据进行分析，得到气温、气压、相对湿度之间的关系如下：表4 回归效果检验表Anovad模型平方和df均方FSig.1回归.0038.0003.866.002a残差.00339.000总计.006472回归.0037.0004.51
21、1.001b残差.00340.000总计.006473回归.0026.0004.656.001c残差.00441.000总计.00647a. 预测变量: (常量), 最小相对湿度, 最低气压, 平均相对湿度, 最高气压, 最高温度, 最低温度, 平均气压, 平均温度。b. 预测变量: (常量), 最小相对湿度, 最低气压, 平均相对湿度, 最高气压, 最高温度, 平均气压, 平均温度。c. 预测变量: (常量), 最小相对湿度, 最低气压, 平均相对湿度, 最高温度, 平均气压, 平均温度。d. 因变量: 患病比例由以上结论可知：p0.05，气温、气压、相对湿度具有统计学意义。表5 脑卒中与气
22、象资料回归分析SPSS20回归分析系数a模型非标准化系数标准系数tSig.相关性共线性统计量B标准 误差试用版零阶偏部分容差VIF1(常量)0.0510.6180.0820.935平均气压0.0030.0012.052.3620.023-0.0150.3540.2830.01952.645最高气压00.001-0.754-1.0050.321-0.042-0.159-0.120.02539.383最低气压-0.0020.001-1.267-3.6950.001-0.105-0.509-0.4420.1228.226平均温度0.0020.0011.3361.1940.240.0480.1880.
23、1430.01187.492最高温度-0.0020.001-1.336-2.4030.0210.004-0.359-0.2870.04621.611最低温度00.0010.2870.2820.7790.0750.0450.0340.01472.452平均相对湿度-0.0010-0.716-4.1810-0.312-0.556-0.50.4872.053最小相对湿度0.00100.4992.2350.0310.0110.3370.2670.2873.4892(常量)0.0950.5910.160.873平均气压0.0030.0012.1592.8110.008-0.0150.4060.3320.
24、02442.206最高气压-0.0010.001-0.894-1.6030.117-0.042-0.246-0.1890.04522.263最低气压-0.0020.001-1.293-3.9490-0.105-0.53-0.4670.137.668平均温度0.0020.0011.6052.7760.0080.0480.4020.3280.04223.921最高温度-0.0020.001-1.388-2.680.0110.004-0.39-0.3170.05219.203平均相对湿度-0.0010-0.719-4.2480-0.312-0.558-0.5020.4882.049最小相对湿度0.0
25、0100.5262.6380.0120.0110.3850.3120.3512.8473(常量)0.0710.6020.1170.907平均气压0.0020.0011.2382.3850.022-0.0150.3490.2870.05418.559最低气压-0.0020-1.093-3.5440.001-0.105-0.484-0.4270.1536.554平均温度0.0020.0011.7913.1040.0030.0480.4360.3740.04422.957最高温度-0.0020.001-1.479-2.820.0070.004-0.403-0.340.05318.972平均相对湿度-
26、0.0010-0.718-4.1640-0.312-0.545-0.5020.4882.049最小相对湿度0.00100.5572.7520.0090.0110.3950.3310.3542.821a. 因变量: 患病比例图六 回归标准化残差的P-P图分析以上结果得到最优化的回归方程：图七 2007-2010年月平均气压与患病比例的双轴变化图图八 2007-2010年月平均相对湿度与患病比例的双轴变化图图九 2007-2010年月平均温度与患病比例的双轴变化图由回归方程结合折线图可得出各因素与脑卒中发病率的影响如下：（1）脑卒中发病率与平均温度成正相关，即平均温度越高，发病率越高。与最高温度成负相关，即最高温度越低，发病率越高。（2）脑卒中发病率与
