统计软件介绍以及R语言数学建模实例推荐.ppt

男性发病人数占总和发病人数的54.18%，女性发病人数占总和发病人数的45.82%，男性发病人数要高于女性发病人数。 2．发病人数与年龄关系的条形图 datasetage-read.csv(age.csv) #读取数据 library(RColorBrewer) #使用分类调色板，将程序包调入内存 death-t(datasetage) #将数据集进行转置 barplot(death, col = brewer.pal(3, Set1), beside = TRUE ,legend = TRUE) #选用“Set1”作为颜色，beside=TRUE 表示条形图并列排列，legend表示显示相关旁注 由图2看出，60-89岁年龄段的人发病人数最多，而这段的人都是老年人。可见老年人是脑卒中发病的主要群体。而44岁以下和90岁以上的人当中发病人数占很小的比例。由于90岁以上人口原本相对较少，所以我们可以认为44岁以下即年轻人发病的比率不是很高。 从性别上看，74岁以下的发病人数中，男性的发病人数多于女性。而75岁以上的发病人数中，男女的比率相当。 3．发病人数与职业关系的条形图 datasetoccupation-read.csv(occupation.csv) par(mfrow=c(2,1)) #并列图，上下并列排列 dataset2007-read.table(clipboard) barplot(t(dataset2007),col=blue,main=2007) #条形图，颜色为蓝色，主标题为“2007” dataset2008-read.table(clipboard) barplot(t(dataset2008),col=blue,main=2008) par(mfrow=c(2,1)) dataset2009-read.table(clipboard) barplot(t(dataset2009),col=blue,main=2009) dataset2010-read.table(clipboard) barplot(t(dataset2010),col=blue,main=2010) datasettotal-read.table(clipboard) barplot(t(datasettotal),col=blue,main=2007-2010) 图3，从各年的发病人数以及07-10的总发病人数中看出，的农民的发病人数远远高于其他职业，而教师，渔民，医护人员中的发病人数非常少。其他职业当中发病人数由高到低依次为退休人员、工人、离退人员、雇员。 为研究脑卒中发病率与气温、气压和相对湿度之间的关系，运用R软件进行操作，建立多元回归模型 代码如下： lm.reg=lm(num~averpres+highpres+lowpres+avertemp+hightemp+lowtemp+comrh,data=dataset2) summary(lm.reg) 通过表1 R语言输出的结果可以看出，只有相对湿度（comrh）的p值小于0.05，通过了显著水平的检验。所以上述的多远回归模型还不能还不能完全解释发病数与其他相关因素的关系。 代码如下： library(DAAG) #将需要的程序包载入内存 vif（lm.reg) #各变量的方差膨胀因子 通过表2的结果可以看出，变量（averpres、highpres、lowpres、avertemp、hightemp、lowtemp）的方差膨胀因子均远大于10，则它们之间可能存在着多重线性关系。所以，需要排除多重线性的关系，然后再进行多元回归。