介绍:逻辑回归是一种广义线性模型,主要用于解决二分类问题。当因变量为生存与死亡。是否发病等二分类的类别变量时,一般的线性回归方程是不合适的。逻辑回归是将因变量等于1的概率作logit变换后,再与自变量建立线性回归方程。
假设:残差独立;因变量服从伯努利分布
模型:
,其中
即
用极大似然法估计β
所用函数:glm(formula,family,data)
formula为回归方程的公式
family为广义线性模型中的概率分布(连接函数)
逻辑回归中需设置family=binomial
与glm()
连用的常用函数如下:summary()
展示拟合模型的细节coefficients()、coef()
拟合模型的参数confint()
模型参数的置信区间residuals()
拟合模型的残差值plot()
评价模型的诊断图predict()
用拟合模型进行预测
代码示例:
x1 <- rnorm(20,mean = 3,sd = 1)
x2 <- rnorm(20,mean = 5,sd = 1)
y <- rbinom(20,size=1,prob = 0.5)
data1 <- data.frame('x1'=x1,'x2'=x2,'y'=y)
#拟合逻辑回归模型
fit1<-glm(y~x1+x2-1,family = binomial,data=data1)
#查看模型细节(参数估计、拟合效果)
summary(fit1)