Cox回归对象中分组观测数的存储位置及提取咨询

• Cox回归对象（coxph()输出的结果）本身不存储分组观测数，它仅保留与模型拟合相关的统计量（如系数、p值、总样本量等），分组样本量这类信息不会被保存在回归结果里。
• gtsummary包的add_n()函数能显示分组观测数，是因为它会回溯到你输入回归的原始数据集df，对分组变量score的每个类别直接统计观测数，而非从res这个回归对象中读取。

手动提取分组观测数的方法 #

如果需要单独提取分组样本量，可以直接对原始数据集操作，比如用dplyr：

library(dplyr)
df %>% 
  group_by(score) %>% 
  summarise(分组观测数 = n())

从gtsummary表格对象中提取 #

如果已经生成了带n的汇总表格，也可以用gtsummary的extract_n()函数提取分组样本量信息：

tbl <- res %>% tbl_regression(exp = T) %>% add_n(location = 'level')
extract_n(tbl)

补充：coxph对象里的res$n返回的是总观测数，这是回归过程中会用到的统计量，所以被保留了。

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