第3章-广义线性模型PPT课件

2020/5/22,.,1,主编：费宇中国人民大学出版社,2020/5/22,.,2,第3章广义线性模型,3.1广义线性模型概述3.2Logistic模型3.3对数线性模型,2020/5/22,.,3,3.1广义线性模型概述,第2章我们研究了多元线性模型,该模型的一个重要假定是因变量是连续型的变量(通常假定服从正态分布),但在许多情况下,这种假定并不合理,例如下面这两种情况.(1)结果变量可能是类型变量.二值分类变量和多分类变量.(比如:是/否,差/一般/良好/优秀等)显然都不是连续型变量.,2020/5/22,.,4,3.1广义线性模型概述,(2)结果变量可能是计数型变量(比如:一周交通事故的数目)这类变量都是非负的有限值,而且它们的均值和方差通常是相关的(一般线性模型假定因变量是正态变量,而且相互独立).普通线性回归模型(2.3)假定因变量y服从正态分布,其均值满足关系式:=X,这表明因变量的条件均值是自变量的线性组合.本章介绍两种常见的广义线性模型:Logistic模型与对数线性模型.,2020/5/22,.,5,3.1广义线性模型概述,1.广义线性模型的定义：(1)随机成分:设y1,y2,yn是来自于指数分布族的随机样本,即yi的密度函数为其中ai(.),b(.),ci(.)是已知函数,参数i是典则参数,是散度参数.,2020/5/22,.,6,1.广义线性模型的定义：,(2)联结函数:设yi的均值为i而函数m(.)是单调可微的联接函数,使得其中是协变量,是未知参数向量.,2020/5/22,.,7,指数分布族,正态分布二项分布泊松分布,2020/5/22,.,8,2.正态线性回归模型,正态分布属于指数分布族,其密度函数为与(3.1)对照可知,2020/5/22,.,9,2.正态线性回归模型,只要取联结函数为,则正态线性回归模型满足广义线性模型的定义.类似的,容易验证,二项分布和泊松分布都属于指数分布族.下面介绍实际中应用广泛的两种广义线性模型:Logistic模型和对数线性模型.,2020/5/22,.,10,3.2Logistic模型,1.模型定义设yi服从参数为pi的二项分布,则i=E(yi)=pi采用逻辑联结函数,即这个广义线性模型称为Logistic模型.,2020/5/22,.,11,例3.1(数据文件为eg3.1),表3.1某地区45个家庭的调查数据,2020/5/22,.,12,2.模型的参数估计和检验,采用R软件中的广义线性模型过程glm()可以完成回归系数的估计,以及模型回归系数的显著性检验.程序如下：,#eg3.1广义线性模型:Logistic模型#打开数据文件eg3.1.xls,选取A1:B46区域,然后复制data3.1-read.table(clipboard,header=T)#将eg3.1.xls数据读入到data3.1中glm.logit-glm(yx,family=binomial,data=data3.1)#建立y关于x的logistic回归#模型,数据为data3.1summary(glm.logit)#模型汇总,给出模型回归系数的估计和显著性检验等ypypp.fit-exp(yp)/(1+exp(yp);p.fit#估计x=15时y=1的概率10.9665418,容易看出:回归模型的回归系数在5%水平上显著,于是得回归模型为当x=15时,估计y=1的概率约为0.97,即年收入为15万元的家庭有私家车的可能性约为97%.,2020/5/22,.,15,3.3对数线性模型,1.模型的定义设y服从参数为的泊松分布,则=E(y)=,采用对数联结函数,即这个广义线性模型称为泊松对数线性模型.,2020/5/22,.,16,例3.2(数据文件为eg3.2),表3.4Breslow癫痫数据,2020/5/22,.,17,例3.2(数据文件为eg3.2),这个数据是robust包中的Breslow癫痫数据(Breslow,1993).我们讨论在治疗初期的八周内,癫痫药物对癫痫发病数的影响,响应变量为八周内癫痫发病数(y),预测变量为前八周内的基础发病次数(x1),年龄(x2)和治疗条件(x3),其中治疗条件是二值变量,x3=0表示服用安慰剂,x3=1表示服用药物.根据这个数据建立泊松对数线性模型并对模型的系数进行显著性检验.,2020/5/22,.,18,表3.2Breslow癫痫数据,2020/5/22,.,19,2.模型的参数估计和检验,采用R软件中的广义线性模型过程glm()来建立泊松对数线性模型并对模型的系数进行显著性检验.程序如下：,#eg3.2广义线性模型:泊松对数线性模型#打开数据文件eg3.2.xls,选取A1:E60区域,然后复制data3.2-read.table(“clipboard”,header=T)#将eg3.2.xls数据读入到data3.2中glm.ln|z|)(Intercept)1.94882590.135619114.3702e-16*x