Como calcular as taxas de risco previstas de um modelo Cox PH?

11

Eu tenho o seguinte modelo Cox PH:

(Hora, Evento) ~ X + Y + Z

Gostaria de obter os perigos previstos taxas (eu estou falando sobre taxas de risco NÃO taxas de risco) dado valores específicos de X, Y, Z. Sei que o pacote muhaz R pode calcular as taxas de risco observadas, mas estou interessado no modelo previsto.

Existe uma maneira de fazer isso no R?

sem estoque
fonte
2
Parece que sua "equação" modela um resultado binário. Agora, ofertas modelo PH de Cox com dados de sobrevivência e que se encontra-se o lado esquerdo deve ser uma taxa de risco ... Eu acho que você deve esclarecer ...
Ocram
2
λ0 0(t)
1
Há alguma ambiguidade na sua pergunta. A função de risco é uma probabilidade, não uma taxa. É a função de risco que você deseja ou algo mais?
Fomite 27/03

Respostas:

16

A função no pacote de sobrevivência R para obter a taxa básica de risco é basehaz.

eβ

Um exemplo simples pode ajudar:

library(survival) #survival analysis
library(eha) #used for data 
data(oldmort) #create the data

# Create surv data set
mort <- Surv(time=oldmort$enter,time2=oldmort$exit,event=oldmort$event)

reg_fit <- coxph(formula=mort~oldmort$sex)
summary(reg_fit)

# Now get baseline curve
baseline <- basehaz(reg_fit)

# Draw baseline hazard (that's male)
plot(baseline$time, baseline$hazard, type='l',main="Hazard rates") 

# Draw female hazard
lines(baseline$time, exp(-0.1929)*baseline$hazard, col="blue") 

exp(-0,1929)=0,8245summary(reg_fit)18%

CarrKnight
fonte
13

A função basehaz(da resposta anterior) fornece o risco cumulativo , não a função de risco (a taxa). Eu acredito que essa pergunta era sobre a função de risco. A estimativa da função de risco exigiria a especificação do tipo de suavização (como na estimativa de densidade). O Muhazpacote R pode fazer isso para uma amostra de dados. Não conheço uma função que possa fazer isso para o risco de linha de base no modelo Cox. Eu também preciso disso. Acho que precisarei diminuir o risco cumulativo por conta própria.

user21413
fonte