如何在R语言中将已计算的Odds Ratio及置信区间整合到新表格？

步骤1：批量提取各模型的OR与置信区间 #

先写一个小函数自动处理每个模型，避免重复代码：

# 定义函数：提取模型的OR、95%置信区间（自动排除截距项）
get_or_ci <- function(model) {
  # 计算OR和置信区间
  or_ci <- exp(cbind(OR = coef(model), confint(model)))
  # 移除截距行（如果不需要截距的OR）
  or_ci <- or_ci[rownames(or_ci) != "(Intercept)", , drop = FALSE]
  # 转成数据框方便后续合并
  as.data.frame(or_ci) %>% 
    rename(`95% CI下限` = 2, `95% CI上限` = 3)
}

# 处理4个模型，得到各自的OR结果表
or_results <- list(
  get_or_ci(table1),
  get_or_ci(table2),
  get_or_ci(table3),
  get_or_ci(table4)
)

步骤2：与已有data表格合并 #

假设你的data表格是按table1-table4的顺序排列的4行数据，直接将OR结果按行合并即可：

# 将OR结果转为矩阵，匹配data的行顺序
or_matrix <- do.call(rbind, or_results)

# 合并到原data表格中
final_table <- cbind(data, or_matrix)

可选：优化表格可读性 #

如果想把OR和CI合并成一个更直观的列，可以加这一步：

library(dplyr) # 需要dplyr包

final_table <- final_table %>% 
  mutate(`OR (95% CI)` = sprintf("%.2f (%.2f-%.2f)", OR, `95% CI下限`, `95% CI上限`)) %>% 
  select(-OR, -`95% CI下限`, -`95% CI上限`) # 可选：移除单独的OR和CI列

注意事项 #

• 如果confint()函数报错，可尝试增加迭代次数：confint(model, nsim = 1000)
• 确认data的行顺序和table1-table4的顺序完全对应，避免合并错位
• 如果需要保留截距项的OR，删除函数里rownames(or_ci) != "(Intercept)"的判断即可

内容的提问来源于stack exchange，提问作者irene123