Stata17自定义表格添加失效数N_fail报r(111)错误怎么解决？

• _r_b、_r_p、_r_se、_r_ci属于系数层级的内置返回结果，collect命令可直接识别这类前缀的系数级指标，自动匹配到对应的变量、因子水平下。
• N_fail是模型整体层级的返回标量，存储在Cox模型运行后的返回值e(N_fail)中，不属于系数级返回结果，因此不能直接和系数级指标并列放在collect的首行收集列表中，否则会触发未找到的r(111)报错。

完成系数级指标收集后，单独调用collect get命令将模型级的失效数添加到对应collection即可，完整可运行代码如下：

* 先清空同名collection避免冲突
collect drop model3

* 第一步：收集系数层级的核心指标
collect _r_b _r_p _r_se _r_ci,  ///
    name(model3) ///
    tag(model[(1)]) ///     
    : stcox i.tertile age_at_baseline initial_dose female_sex

* 第二步：单独添加模型级失效数指标，自定义设置标签
collect get e(N_fail), name(model3) tag(model[(1)]) label("Number of failures")

后续可以通过collect layout命令调整表格布局，将N_fail放在表格底部的汇总行展示即可。

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