🔧
tool
相关软件(R代码)
```r # Logistic回归 # install.packages("MASS") library(MASS) data(birthwt) # 建立Logistic回归模型 model_logit <- glm(low ~ age + lwt + race + smoke + ht, data = birthwt, family = binomial(...
📖 定义
# Logistic回归
# install.packages("MASS")
library(MASS)
data(birthwt)
# 建立Logistic回归模型
model_logit <- glm(low ~ age + lwt + race + smoke + ht,
data = birthwt, family = binomial(link = "logit"))
summary(model_logit)
# 提取OR及其置信区间
exp(cbind(coef(model_logit), confint(model_logit)))
# 模型评估:ROC曲线
# install.packages("pROC")
library(pROC)
prob <- predict(model_logit, type = "response")
roc_obj <- roc(birthwt$low, prob)
plot(roc_obj, main = paste("AUC =", round(auc(roc_obj), 3)))
# Probit模型
model_probit <- glm(low ~ age + lwt + race + smoke + ht,
data = birthwt, family = binomial(link = "probit"))
summary(model_probit)