OK dari kode R Anda mengasumsikan distribusi eksponensial (bahaya konstan) untuk bahaya baseline Anda. Karena itu fungsi bahaya Anda:
h ( t ∣ Xsaya) = { exp( α β0)exp( γ+ α ( β0+ β1+ β2t ) )jika Xsaya= 0 ,jika Xsaya= 1 .
Kami kemudian mengintegrasikan ini sehubungan dengan untuk mendapatkan fungsi bahaya kumulatif:t
Λ ( t ∣ Xsaya)= { t exp( α β0)∫t0exp( γ+ α ( β0+ β1+ β2τ) )dτjika Xsaya= 0 ,jika Xsaya= 1 .= { t exp( α β0)exp( γ+ α ( β0+ β1) ) 1α β2( exp( α β2t ) - 1 )jika Xsaya= 0 ,jika Xsaya= 1 .
Ini kemudian memberi kita fungsi bertahan hidup:
S( t )= exp( - Λ ( t ) )= { exp( - t exp(α β0) )exp( - exp( γ+ α (β0+β1) )1αβ2(exp( αβ2t ) - 1 ) )jika Xsaya= 0 ,jika Xsaya= 1 .
Anda kemudian menghasilkan dengan mengambil sampel dan , menggantikan untuk dan menyusun ulang rumus yang sesuai (berdasarkan ) untuk mensimulasikan . Aljabar ini harus langsung Anda dapat kode di R tetapi tolong beri tahu saya dengan komentar jika Anda memerlukan bantuan lebih lanjut. U ∼ U n i f o r m ( 0 , 1 ) U S ( t ) X i tXsayaU∼Uniform(0,1)US(t)Xit