Saya ingin menghasilkan waktu bertahan hidup dari model bahaya proporsional Cox yang mengandung kovariat tergantung waktu. Modelnya adalah
di mana dihasilkan dari Binomial (1,0.5) dan .m i ( t ) = β 0 + β 1 X i + β 2 X i t
Nilai parameter sebenarnya digunakan sebagai
Untuk kovariat independen waktu (yaitu saya membuat sebagai berikut
#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)
Adakah yang bisa membantu saya menghasilkan data survival dengan kovariat yang bervariasi waktu.
Jawaban:
OK dari kode R Anda mengasumsikan distribusi eksponensial (bahaya konstan) untuk bahaya baseline Anda. Karena itu fungsi bahaya Anda:
Kami kemudian mengintegrasikan ini sehubungan dengan untuk mendapatkan fungsi bahaya kumulatif:t
Ini kemudian memberi kita fungsi bertahan hidup:
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 tXsaya U∼Uniform(0,1) U S(t) Xi t
sumber