这个例子就是这样的
多元的
用到广义线性回归 glm()函数
# get summary statistics
data(Affairs, package = "AER")
summary(Affairs)
table(Affairs$affairs)
# create binary outcome variable
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 model
fit.full <- glm(ynaffair ~ gender + age + yearsmarried +
children + religiousness + education + occupation + rating,
data = Affairs, family = binomial())
summary(fit.full)
# fit reduced model
fit.reduced <- glm(ynaffair ~ age + yearsmarried +
religiousness + rating, data = Affairs, family = binomial())
summary(fit.reduced)
|