Logistic模型统计分析的R软件实现.doc

Logistic模型统计分析的R软件实现 2011年12月 第4期 连云港师范高等专科学校 JournalofLianyungangTeachersCollege December,2011 No.4 Logistic模型统计分析的R软件实现 光琳 (江苏联合职业技术学院连云港财经分院,江苏连云港222003) [摘要]随着回归分析理论的不断发展,Logistic回归模型已成为目前重要的研究课题,并在医学,工程技术等领 域得到了广泛应用.它的数据计算依赖于统计软件,文章结合实例研究了如何在R软件中实现对Logistic模型的统计分 析. [关键词]Logistic模型;R软件;统计诊断 [中图分类号]o212.1[文献标识码]A[文章编号]1009—7740(2011)04—0106—03 1引言 Logistic模型最早由比利时生物学家Verhulst于 1838年导出,该模型形式简单但在生物,经济,环境科 学等领域都有着广泛应用,特别是在人口数量预测, 耐用消费品销售量预测等方面有着重要的应用.它主 要研究二分变量与诸多自变量的相互关系.随着理论 研究的不断深入,Logistic模型也生发出多种形式,其 中如下形式:… ,n, logit(P)=logIJ=flo+X1+X2+…pp—l’上, X(*) 在工程以及医学,流行病学,生物学上都有广泛的应 用价值,13益受到人们的重视. 由于Logistic回归模型涉及较复杂的数学理论, 数据统计分析的计算往往较为复杂,一般无法用手工 计算,L2只能在计算机上实现.目前的统计软件有很 多,如SAS,SPSS,STAT,S—PLUS等,但大多是商业性 的,购买需支付昂贵的费用.而R软件是一款免费的 统计软件,它提供了有弹性的,互动的环境来分析,展 示数据,_3且提供若干统计程序包以及一些集成统计 计算工具和函数,使得用户可以灵活机动地进行数据 分析,简化了数据分析过程.它可以完成大多数模型 的统计计算,并帮助用户分析计算结果.本文将结合 实例,展示如何在R软件中实现对Logistic模型的统 计分析. 2Logistic模型的拟合 回归模型(*)的拟合主要是求模型中的参数估 计值,Logistic模型的参数估计一般采用极大似然估 计,对于得到的似然方程xT(Y—P)10:=0,可采用 Gauss—Newton迭代法口进行计算.鉴于Logistic模型 是基于二项分布族的广义线性模型,因此在R软件中 可通过glm语句建立回归关系,再用summary语句得 到其详细结果. 3Logistic模型的诊断 对得到的模型拟合结果,R软件还可以对其进行 诊断,可用influence.measures语句进行影响分析.通 过该语句可以直接得到像cook距离,帽子矩阵对角 元值hii等基本统计量.再通过函数运算可得到其他统 计量,进而进行诊断分析. 4数值实例 资料为5O位急性白血病人的生存时间Y,入院后 取得的淋巴浸润等级x以及出院后有无巩固治疗x2 的关系.b该数据曾被多次引用,周雁曾利用S— plus和Matlab软件进行拟合,而谢书培J也只用了30 例数据,得到的结果不够详细. 根据数据,用R软件拟合Logistic模型,具体过程 如下: Xlt;一c(0,2,2,2,2,0,0,2,2,0,2,0,0,0,0,0,0, 0,0,2,0,0,0,0,0,2,0,0,0,0,0,2,0,2,0,0,2,0,2,0, 0,0,2,0,0,2,O,0,0,0) x2lt;一c(0,0,0,0,0,0,1,0,0,0,1,0,0,0,1,0,1, 1,l,1,1,0,1,1,1,0,0,0,0,0,0,0,1,0,1,1,1,0,1,0, 0,1,1,1,1,1,1,1,1,1) Ylt;一C(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1, [收稿日期]2011—09—08 [作者简介]光琳(1980一),男,河南灵宝人,江苏联合职业技术学院连云港财经分院讲师,扬州大学数学科学学院硕士 生. 1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 1,1,1,1,1,1,1,1,1,1) lifelt;一data.fr~ane(X1,x2,Y) loglt;一glm(Y一?Xl+x2,family=binomial,data= life) summary(1og) 根据拟合的结:果,得到 l0g()=一1.64l9—0.7070X1+2.7844X2, 即 . exp(一1.一6419—0.7070X1+2.7844X2) 一 1+exp(一1.6419—0.7070X1+2.7844x2)’ R软件可以用来进行预测分析:-3] pre