0
1)如何在以下示例中将y轴更改为“比值比”,“死亡概率”和“死亡率”?如何绘制带有危险无线电,死亡概率或y轴死亡率的受限立方样条?
2)在以下示例中,如何将y轴更改为“fit2”的“风险比”?
library(Hmisc)
library(survival)
library(rms)
data(pbc)
d <- pbc
rm(pbc)
d$died <- ifelse(d$status == 2, 1, 0)
d$status <- ifelse(d$status != 0, 1, 0)
ddist <- datadist(d)
options(datadist='ddist')
fit <- lrm(status ~ rcs(age, 4), data=d)
(an <- anova(fit))
plot(Predict(fit), anova=an, pval=TRUE)
fit2 <- cph(Surv(time, status) ~ rcs(age, 4), data=d)
(an2 <- anova(fit2))
plot(Predict(fit2), anova=an, pval=TRUE)
我期待着您的帮助!
更新1 继BondedDust的答案,我又写道如下:
# probability
getProbability <- function(x) {
exp(x)/(1+exp(x))*100
}
fit <- lrm(status ~ rcs(age, 4), data=d)
(an <- anova(fit))
plot(Predict(fit, fun=getProbability), anova=an, pval=TRUE, ylab="Probability of death [%]")
# overall probability to die
table(d$status)
round(table(d$status)[[2]]/sum(table(d$status))*100, digits=1) # = 44.5%
由于总体概率死了,是44.5%的预测概率的计算所得的情节似乎是正确的,以我为非统计学家,不是吗?