🔧 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)

2.3.4 泊松回归模型