अपने आर कोड से ठीक है आप अपने बेसलाइन खतरे के लिए एक घातांक वितरण (निरंतर खतरा) मान रहे हैं। आपके खतरनाक कार्य इसलिए हैं:
एच ( टी ∣)एक्समैं) = {exp( αβ0)exp( γ+ α (β0+β1+β2टी ) )अगर एक्समैं= 0 ,अगर एक्समैं= 1 ।
इसके बाद हम इन्हें एकीकृत करते हैं टी संचयी खतरा समारोह पाने के लिए:
Λ ( टी |एक्समैं)= {टी ऍक्स्प( αβ0)∫टी0exp( γ+ α (β0+β1+β2τ) )घτअगर एक्समैं= 0 ,अगर एक्समैं= 1 ।= {टी ऍक्स्प( αβ0)exp( γ+ α (β0+β1) )1αβ2( ऍक्स्प( αβ2टी ) - 1 )अगर एक्समैं= 0 ,अगर एक्समैं= 1 ।
ये तब हमें जीवित रहने के कार्य देते हैं:
एस( टी )=exp(−Λ(t))={exp(−texp(αβ0))exp(−exp(γ+α(β0+β1))1αβ2(exp(αβ2t)−1))if Xi=0,if Xi=1.
आप फिर नमूना लेकर उत्पन्न करते हैं Xi तथा U∼Uniform(0,1), स्थानापन्न U के लिये S(t) और उपयुक्त सूत्र को पुन: व्यवस्थित करना (आधार पर) Xi) अनुकरण करने के लिए t। यह सीधा बीजगणित होना चाहिए जिसे आप R में कोड कर सकते हैं, लेकिन कृपया मुझे कमेंट द्वारा बताएं यदि आपको और मदद की आवश्यकता हो तो।