r - 匹配对的条件逻辑回归

Question

我正在尝试对类似于以下示例的数据运行条件逻辑回归:

table.10.3 <-  data.frame(pair=rep(1:144,rep(2,144)),
                         MI=rep(c(0,1),144),
                         diabetes=c(rep(c(1,1),9),
                                    rep(c(1,0),16),
                                    rep(c(0,1),37),
                                    rep(c(0,0),82))
                          )
# head(table.10.3)
# pair MI diabetes
#    1  0        1
#    1  1        1
#    2  0        1
#    2  1        1
#    3  0        1
#    3  1        1

library("survival")
fit.CLR <- clogit(MI ~ diabetes + strata(pair), method="exact", data=table.10.3)

summary(fit.CLR)

我得到了总结。我的问题是如何以图形方式表示结果？我需要举行，因为我对 R 绘图很陌生。我试过 vcd 包。我能够获得其他一些虚拟数据的马赛克图。但我想绘制 clogit 模型的结果。

Answer 1

我认为您无法通过 clogit 获得生存曲线。条件逻辑回归不会自动考虑生存时间；它只处理包含匹配案例的层的成员资格，并控制 Cox 模型处理生存时间的方式(因此它出现在“生存”包中)。

如果您对估计基线生存和风险函数感兴趣，我认为您需要使用像 coxph 这样的生存模型，或者您可以直接使用 survfit 估计生存曲线。