आर का उपयोग करके समय निर्भर सहसंयोजकों के साथ उत्तरजीविता डेटा कैसे उत्पन्न करें


9

मैं एक कॉक्स आनुपातिक खतरों मॉडल से बचने का समय उत्पन्न करना चाहता हूं जिसमें समय पर निर्भर कोवरिएट होता है। मॉडल है

(टी|एक्समैं)=0(टी)exp(γएक्समैं+αमैं(टी))

जहाँ Binomial (1,0.5) और ।एक्समैंमैं(टी)=β0+β1एक्समैं+β2एक्समैंटी

सही पैरामीटर मानों का उपयोगγ=1.5,β0=0,β1=1,β2=1.5,h0(t)=1

समय-स्वतंत्र सहसंयोजक के लिए (यानी मैं इस प्रकार उत्पन्नh(t|Xi)=h0(t)exp(γXi)

#For time independent case
# h_0(t) = 1
gamma <- -1
u <- runif(n=100,min=0,max=1)
Xi <- rbinom(n=100,size=1,prob=0.5)
T <- -log(u)/exp(gamma*Xi)

किसी को भी समय-अलग-अलग सहसंयोजक के साथ उत्तरजीविता डेटा बनाने के लिए मेरी मदद कर सकते हैं।


किस प्रकार का कार्य है मैं(टी)? क्या यह निरंतर है? टुकड़ों में स्थिर? एक अलग एल्गोरिथ्म शायद उसी के अनुसार की आवश्यकता होगी।
ट्रिस्टन

mi(t)एक समय निर्भर सहसंयोजक है, सादगी के लिए आप समय के साथ आनुपातिक संबंध पर विचार कर सकते हैं।
शेख

मैंने अपने प्रश्न को संपादित किया है, जिसके एक फंक्शन पर विचार किया जा रहा है mi(t)
शेख

उपरोक्त समीकरण से आपने R कोड कैसे किया? इसका मतलब है कि एक ही आईडी के भीतर प्रत्येक मृत्यु के समय कार्यक्रम को यह पता लगाने की आवश्यकता होती है कि कोवरिएट्स हर किसी के लिए हैं जो या तो x 1 या 0. के बराबर है यदि सभी 1 कम्बल के बराबर है। उसके बाद उत्तरजीविता फ़ंक्शन की गणना करें। यह प्रत्येक विषय के लिए सही लाइन लेने देता है।
क्स आमल

जैसा कि जेड झांग बताते हैं तो इस लेख पर एक नज़र डालते हैं । इसके अलावा, आप मेरे प्रश्न का उत्तर देख सकते हैं, जहाँ मैं दिखाता हूँ कि उन लोगों के लिए कैसे अनुकरण किया जाएएक्समैं=1आर। में समूह
बेंजामिन क्रिस्टोफ़र्सन

जवाबों:


9

अपने आर कोड से ठीक है आप अपने बेसलाइन खतरे के लिए एक घातांक वितरण (निरंतर खतरा) मान रहे हैं। आपके खतरनाक कार्य इसलिए हैं:

(टी|एक्समैं)={exp(αβ0)अगर एक्समैं=0,exp(γ+α(β0+β1+β2टी))अगर एक्समैं=1

इसके बाद हम इन्हें एकीकृत करते हैं टी संचयी खतरा समारोह पाने के लिए:

Λ(टी|एक्समैं)={टीexp(αβ0)अगर एक्समैं=0,0टीexp(γ+α(β0+β1+β2τ))τअगर एक्समैं=1={टीexp(αβ0)अगर एक्समैं=0,exp(γ+α(β0+β1))1αβ2(exp(αβ2टी)-1)अगर एक्समैं=1

ये तब हमें जीवित रहने के कार्य देते हैं:

S(t)=exp(Λ(t))={exp(texp(αβ0))if Xi=0,exp(exp(γ+α(β0+β1))1αβ2(exp(αβ2t)1))if Xi=1.

आप फिर नमूना लेकर उत्पन्न करते हैं Xi तथा UUniform(0,1), स्थानापन्न U के लिये S(t) और उपयुक्त सूत्र को पुन: व्यवस्थित करना (आधार पर) Xi) अनुकरण करने के लिए t। यह सीधा बीजगणित होना चाहिए जिसे आप R में कोड कर सकते हैं, लेकिन कृपया मुझे कमेंट द्वारा बताएं यदि आपको और मदद की आवश्यकता हो तो।


1
बीजगणित के लिए बहुत-बहुत धन्यवाद। मैं आर में कोड करूंगा और आगे की मदद के लिए आपसे संपर्क करूंगा।
शेख

क्या एक सही जवाब, @tristan। मेरे पास एक समान प्रश्न था और आपका उत्तर मिला। बस कमाल।
सैम

@tristan मैं पहले समीकरण में अल्फा के अर्थ के बारे में थोड़ा उलझन में हूँ, जहाँ आप शी = 0. देते हैं, क्या आप उस पर थोड़ा विस्तार करना चाहेंगे? धन्यवाद।
स्टैटवॉन्क

1
@Statwonk यह मूल पोस्टर द्वारा प्रदत्त खतरे दर समीकरण से निम्नानुसार है
ट्रिस्टन

क्षमा करें, लेकिन मुझे यकीन नहीं है कि समय का अनुकरण करने के लिए फ़ंक्शन S (t) का उपयोग कैसे करें। मुझे लगता है कि आपको S ^ {- 1} की गणना करनी चाहिए और यह कार्य X_i = 1 के मामले के लिए तुच्छ नहीं है।
पीएमसी
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.