R语言-logistics回归分析

xujunzju • 2023年6月18日 12:36 • 统计分析, 编程语言 • 阅读 181

使用R语言包AER中的Affairs作为例子

data(Affairs, package = "AER")
View(Affairs)
summary(Affairs)


table((Affairs$affairs))


#分类变量，变成因子
Affairs$ynaffair[Affairs$affairs > 0] <- 1
Affairs$ynaffair[Affairs$affairs == 0] <- 0
Affairs$ynaffair <- factor(Affairs$ynaffair,
                           levels= c(0,1),
                           labels = c("no","yes"))

table((Affairs$ynaffair))

#建模，因子全纳入
fit.full <- glm(ynaffair ~ gender + age + yearsmarried + children +
                  religiousness + education + occupation + rating ,
                data =  Affairs, family = binomial())
summary(fit.full)
coef(fit.full)
#OR值
exp(coef(fit.full))

confint(fit.full)
#95%的置信区间
exp(confint(fit.full))

这个就是OR值