จะคำนวณอัตราอันตรายที่คาดการณ์จากแบบจำลอง Cox PH ได้อย่างไร

ฉันมีโมเดล Cox PH ต่อไปนี้:

(เวลา, กิจกรรม) ~ X + Y + Z

ฉันต้องการที่จะได้รับอันตรายจากการคาดการณ์อัตรา (ฉันพูดคุยเกี่ยวกับอัตราการเกิดอันตรายไม่อัตราส่วนอันตราย) ได้รับค่าเฉพาะของX, ,Y Zฉันรู้ว่าแพ็คเกจ muhaz Rสามารถคำนวณอัตราอันตรายที่สังเกตได้ แต่ฉันสนใจในแบบจำลองที่ทำนายไว้

มีวิธีทำเช่นนี้ใน R หรือไม่?

ฟังก์ชั่นในแพคเกจ R basehazอยู่รอดที่จะได้รับอัตราอันตรายพื้นฐานคือ

จากนั้นคุณต้องเพิ่มทวีคูณสำหรับเพื่อรับอัตราอันตรายเฉพาะตามค่าสัมประสิทธิ์ที่คุณพบ$e^{\beta}$

ตัวอย่างง่ายๆอาจช่วย:

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") 

จะถูกคัดลอกและวางจากผลสรุปของการถดถอย Cox, หมายความว่าหญิงในข้อมูลที่มีอัตราการเกิดอันตรายต่ำกว่าเพศชาย$\exp(-0.1929)=0.8245$summary(reg_fit)$18\%$

ฟังก์ชั่นbasehaz(จากคำตอบก่อนหน้านี้) ให้อันตรายสะสมไม่ใช่ฟังก์ชั่นอันตราย (อัตรา) ฉันเชื่อว่าคำถามนั้นเกี่ยวกับฟังก์ชั่นอันตราย การประมาณฟังก์ชั่นอันตรายจะต้องใช้คุณสมบัติของการปรับให้เรียบ (เช่นในการประมาณความหนาแน่น) Muhazแพคเกจ R สามารถทำเช่นนี้สำหรับข้อมูลตัวอย่างหนึ่ง ฉันไม่ได้ตระหนักถึงฟังก์ชั่นที่สามารถทำสิ่งนี้เพื่ออันตรายพื้นฐานในรูปแบบ Cox ฉันต้องการสิ่งนี้ ฉันคิดว่าฉันจะต้องทำให้อันตรายที่สะสมได้ด้วยตัวเองราบรื่นขึ้น