编辑代码

# 加载 epiDisplay 包
install.packages("epiDisplay") # 若尚未安装，请先运行此行
library(epiDisplay)

# 加载数据集 Ectopic
data(Ectopic)

# 查看数据集的结构
str(Ectopic)
head(Ectopic)

# 数据预处理
# 将 EP (宫外孕) 定义为因变量（1 = case组，0 = control组）
Ectopic$case <- ifelse(Ectopic$EP == 1, 1, 0)

# 定义自变量：人工流产 (IA) 和未分娩的孕妇 (Deli)
# IA 和 Deli 已在数据集中，可直接使用

# Logistic 回归分析
model <- glm(case ~ IA + Deli, data = Ectopic, family = binomial)

# 显示回归结果
summary(model)

# 计算 Odds Ratio（OR）及置信区间
exp(cbind(OddsRatio = coef(model), confint(model)))

# 结果可视化
logistic.display(model)