Prever probabilidade do modelo Cox PH
Estou tentando usar o modelo cox para prever a probabilidade de falha após o tempo (que é chamado de parada) 3.
bladder1 <- bladder[bladder$enum < 5, ]
coxmodel = coxph(Surv(stop, event) ~ (rx + size + number) +
cluster(id), bladder1)
range(predict(coxmodel, bladder1, type = "lp"))
range(predict(coxmodel, bladder1, type = "risk"))
range(predict(coxmodel, bladder1, type = "terms"))
range(predict(coxmodel, bladder1, type = "expected"))
No entanto, as saídas da função de previsão não estão no intervalo de 0 a 1. Existe alguma função ou como posso usar a função de previsão lp e risco de linha de base para calcular a probabilidade?