R语言orscoreci函数参数疑问：为何用病例/对照总数而非未暴露数

研究背景与数据 #

在一项回顾性病例对照研究中，分析肺癌与吸烟的关联，选取709例病例和709例对照，暴露因素为是否曾吸烟的二分类变量，数据如下：

Predictor  Disease1 Disease2 Total
  Exposed1      688      650  1338
  Exposed2       21       59    80
  Total         709      709  1418

手动计算OR及95%置信区间 #

手动计算过程如下：

ln(odds ratio)+- z alpha/2(SE):
ln(2.97)+-1.96(0.2599)

-> 1.089+1.96(0.2599) = 1.5984
-> 1.089-1.96(0.2599) = 0.5795

• 置信区间：(exp(0.5795), exp(1.5984)) = (e^0.5795, e^1.5984) = (1.78 , 4.94)

R语言实现代码 #

指导老师要求使用orscoreci函数完成计算，代码如下：

library('PropCIs')
orscoreci(688, 709, 650, 709, conf.level = 0.95)

核心疑问 #

手动计算时未用到21、59这些未暴露数，为何orscoreci函数要使用病例/对照总数709、709而非21、59？且手动计算与R运行结果均为(1.78 , 4.94)。

内容的提问来源于stack exchange，提问作者Alexandra García