วิธีการสร้างข้อมูลการอยู่รอดด้วยโควาเรียที่ขึ้นอยู่กับเวลาโดยใช้ R

ฉันต้องการสร้างเวลาการเอาชีวิตรอดจากโมเดลอันตรายตามสัดส่วนของค็อกซ์ที่มีเวลาแปรปรวนร่วม รูปแบบคือ

$h(t|X_i) =h_0(t) \exp(\gamma X_i + \alpha m_{i}(t))$

ที่ถูกสร้างขึ้นจากทวินาม (1,0.5) และT$X_i$$m_{i}(t)=\beta_0 + \beta_1 X_{i} + \beta_2 X_{i} t$

ค่าพารามิเตอร์จริงถูกใช้เป็น$\gamma = 1.5, \beta_0 = 0, \beta_1 = -1, \beta_2 = -1.5, h_0(t) = 1$

สำหรับตัวแปรอิสระตามเวลา (เช่นฉันสร้างขึ้นดังนี้$h(t|X_i) =h_0(t) \exp(\gamma X_i)$

#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)

ทุกคนได้โปรดช่วยฉันสร้างข้อมูลการอยู่รอดด้วยตัวแปรแปรผันตามเวลา

ตกลงจากรหัส R ของคุณคุณกำลังสมมติการแจกแจงแบบเอ็กซ์โพเนนเชียล (อันตรายคงที่) สำหรับอันตรายพื้นฐานของคุณ ฟังก์ชั่นอันตรายของคุณคือ:

$$h\left(t \mid X_i\right) = \begin{cases} \exp{\left(\alpha \beta_0\right)} & \text{if $X_i = 0$,} \\ \exp{\left(\gamma + \alpha\left(\beta_0+\beta_1+\beta_2 t\right)\right)} & \text{if $X_i = 1$.} \end{cases}$$

จากนั้นเรารวมสิ่งเหล่านี้ด้วยความเคารพต่อเพื่อให้ได้ฟังก์ชันอันตรายสะสม:$t$

$$\begin{align} \Lambda\left(t\mid X_i\right) &= \begin{cases} t \exp{\left(\alpha \beta_0\right)} & \text{if $X_i=0$,} \\ \int_0^t{\exp{\left(\gamma + \alpha\left(\beta_0+\beta_1+\beta_2 \tau\right)\right)} \,d\tau} & \text{if $X_i=1$.} \end{cases} \\ &= \begin{cases} t \exp{\left(\alpha \beta_0\right)} & \text{if $X_i=0$,} \\ \exp{\left(\gamma + \alpha\left(\beta_0+\beta_1\right)\right)} \frac{1}{\alpha\beta_2} \left(\exp\left(\alpha\beta_2 t\right)-1\right) & \text{if $X_i=1$.} \end{cases} \end{align}$$

สิ่งเหล่านี้ทำให้เรามีฟังก์ชั่นการเอาชีวิตรอด:

$$\begin{align} S\left(t\right) &= \exp{\left(-\Lambda\left(t\right)\right)} \\ &= \begin{cases} \exp{\left(-t \exp{\left(\alpha \beta_0\right)}\right)} & \text{if $X_i=0$,} \\ \exp{\left(-\exp{\left(\gamma + \alpha\left(\beta_0+\beta_1\right)\right)} \frac{1}{\alpha\beta_2} \left(\exp\left(\alpha\beta_2 t\right)-1\right)\right)} & \text{if $X_i=1$.} \end{cases} \end{align}$$

จากนั้นคุณสร้างโดยการสุ่มตัวอย่างและแทนสำหรับและจัดเรียงสูตรที่เหมาะสม (ขึ้นอยู่กับ ) เพื่อจำลองทีนี่ควรเป็นพีชคณิตแบบตรงไปตรงมาจากนั้นคุณสามารถโค้ดใน R แต่โปรดแจ้งให้เราทราบโดยแสดงความคิดเห็นหากคุณต้องการความช่วยเหลือเพิ่มเติม$X_i$$U\sim\mathrm{Uniform\left(0,1\right)}$$U$$S\left(t\right)$$X_i$$t$