कॉक्स PH मॉडल से अनुमानित खतरनाक दरों की गणना कैसे करें?

Question

मेरे पास निम्नलिखित कॉक्स PH मॉडल है:

(समय, घटना) ~ एक्स + वाई + जेड

मैं अनुमानित खतरे की दरें प्राप्त करना चाहूंगा (मैं खतरनाक दरों के बारे में बात कर रहा हूं न कि जोखिम अनुपात) X, के विशिष्ट मूल्य दिए गए हैं Y, और Z। मुझे पता है कि मुहाज़ आर पैकेज मनाया खतरों की दर की गणना कर सकता है, लेकिन मुझे पूर्वानुमानित मॉडल में दिलचस्पी है।

क्या आर में ऐसा करने का कोई तरीका है?

Comments

ऐसा लगता है कि आपका "समीकरण" एक द्विआधारी परिणाम का मॉडल है। अब, कॉक्स का PH मॉडल उत्तरजीविता डेटा के साथ काम करता है और जो झूठ बोलता है वह बाएं हाथ की ओर एक खतरनाक दर होना चाहिए ... मुझे लगता है कि आपको स्पष्ट करना चाहिए ...

— ocram

कॉक्स मॉडल आधारभूत खतरे के कार्य का अनुमान नहीं लगाता है (आमतौर पर रूप में संकेत दिया गया है )। एक मॉडल से घटना दर का अनुमान लगाने के लिए, आपको इसके बजाय पैरामीट्रिक उत्तरजीविता मॉडल का उपयोग करना चाहिए। $\lambda_0(t)$

— बोसोविच

आपके प्रश्न में थोड़ी अस्पष्टता है। खतरा फ़ंक्शन एक संभावना है, दर नहीं। क्या यह खतरनाक कार्य है जिसे आप चाहते हैं, या कुछ और?

— फोमाइट

Answer 1

आधारभूत खतरे की दर प्राप्त करने के लिए 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") 

$\exp(-0.1929)=0.8245$summary(reg_fit)$18\%$

Answer 2

फ़ंक्शन basehaz(पिछले उत्तर से) संचयी खतरा प्रदान करता है , न कि खतरनाक फ़ंक्शन (दर)। मेरा मानना ​​है कि प्रश्न खतरे के कार्य के बारे में था। खतरे के कार्य का अनुमान लगाने के लिए चौरसाई के प्रकार (घनत्व के आकलन में) के विनिर्देश की आवश्यकता होगी। Muhazआर पैकेज एक नमूना डेटा के लिए ऐसा कर सकते हैं। मैं एक फ़ंक्शन के बारे में नहीं जानता जो कॉक्स मॉडल में बेसलाइन के खतरे के लिए ऐसा कर सकता है। मुझे भी इसकी जरूरत है। मुझे लगता है कि मुझे अपने आप ही संचयी खतरे को सुचारू करने की आवश्यकता होगी।