如何根据 rcs 逻辑回归创建参考 OR=1 的 OR 图

Question

我正在尝试使用 rms 包中的脊柱创建一个显示 lrm 模型的优势比的图。我创建了对数赔率和赔率的图，但对于显示赔率比的图，参考从赔率比 0 开始，而不是赔率比 1。

这是我在生存包中使用 cgd 数据集的说明。

summary(cgd)

dd = datadist(cgd)
options(datadist="dd")
attach(cgd)

fit <- lrm(inherit ~ rcs(age, 4), data=cgd)
fit
summary(fit)
an <- anova(fit)
f1 <- ggplot(Predict(fit), ylab="Log-Odds")+ theme_grey()
f2 <- ggplot(Predict(fit, fun=exp), anova=an, pval = TRUE, ylab="Odds ratio") + theme_grey()

cowplot::plot_grid(f1, f2, nrow = 1, ncol = 2, scale = .9)

这些是生成的图：

我的代码有问题吗？我该怎么做才能使比值比从 1 开始？

Answer 1

在第二张图中，您没有绘制几率比，而是绘制 x 连锁与常染色体本身的几率。 优势比，即年龄 = 10 的预测优势与年龄 = 11 的预测优势之间的比率。 基本上，您必须除以赔率（年龄 = 11）/赔率（年龄 = 10）才能获得预测变量中年龄增加 1 分的 OR。