首页 理论教育多元统计分析:对数线性模型的应用

多元统计分析:对数线性模型的应用

【摘要】:对于广义线性模型,除了以上介绍的Logistic回归模型外,还有其他的模型,如Poisson模型,这里就不作详细介绍.以下简要介绍R 软件中“glm( )”关于这些模型的使用方法.Poisson分布族模型和拟Poisson分布族模型的使用方法如下:fm<-glmfm<-glm其直观意义是ln[E]=β0+β1x1+β2x2+…+βp xp,即E=exp(β0+β1 x1+β2 x2+…

对于广义线性模型,除了以上介绍的Logistic回归模型外,还有其他的模型,如Poisson模型,这里就不作详细介绍.以下简要介绍R 软件中“glm( )”关于这些模型的使用方法.

Poisson分布族模型和拟Poisson分布族模型的使用方法如下:

fm<-glm(formula,family=poisson(link=log),data.frame)

fm<-glm(formula,family=quasipoisson(link=log),data.frame)

其直观意义是

ln[E(y)]=β0+β1x1+β2x2+…+βp xp,

E(y)=exp(β0+β1 x1+β2 x2+…+βp xp).

Poisson分布族模型和拟Poisson分布族模型的唯一差别就是:Poisson分布族模型要求响应变量y 是整数,而拟Poisson分布族模型则没有这个要求.(www.chuimin.cn)

对于联列表还可以用(多项分布)对数线性模型来描述.以二维联列表为例,只有主效应的对数线性模型为

ln mij=αi+βj+εij

这相当于只有主效应αi 和βj,而这两个变量的效应是简单可加的.但是有时两个变量在一起时会产生交叉效应,此时相应的对数线性模型为

ln mij=αi+βj+(αβ)ij+εij

对于表中数目代表一个观测数目时,就要考虑是否用Poisson对数线性模型.例如,如果有两个定性变量、一个定量变量的Poisson对数线性模型可以表示为

lnλ=μ+αi+βj+γx+εij,

其中,μ 为常数项;αi 和βj 为两个定性变量的主效应;x 为连续变量;而γ 为其系数;εij 为残差项.这里之所以对Poisson分布的参数λ 取对数,是为了使模型左边的取值范围为整个实数轴