Temel tehlike oranını elde etmek için R hayatta kalma paketindeki fonksiyon basehaz
.
Daha sonra , bulduğunuz katsayılara göre belirli tehlike oranını elde etmek için çeşitli için çarpmanız gerekir .eβ
Basit bir örnek yardımcı olabilir:
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")
kopyalanan Cox regresyon özet sonuçları yapıştırılan, , ima veri kadın tehlike oranları olduğu erkeklere göre daha düşük.% 18exp(−0.1929)=0.8245summary(reg_fit)
18%