匹配对的条件逻辑回归

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

如果您想绘制 S 曲线，下面的代码即可实现。

x <- -7:7
y <- 1 / (1 + exp(-x)) #Sigmoid Logistic function
plot(x,y,col="DarkGreen",pch=19,lwd=1)