Chcę wygenerować czas przeżycia na podstawie proporcjonalnego modelu zagrożeń Coxa, który zawiera zmienną zależną od czasu. Model jest
gdzie jest generowany z Binomial (1,0.5) i .
Prawdziwe wartości parametrów są stosowane jako
Dla niezależnej od czasu (tj. wygenerowałem w następujący sposób
#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)
Czy ktoś może mi pomóc w generowaniu danych o przeżyciu ze zmienną towarzyszącą w czasie.
Odpowiedzi:
OK z kodu R zakładasz rozkład wykładniczy (stałe zagrożenie) dla twojego podstawowego zagrożenia. Do twoich funkcji zagrożeń należą zatem:
Następnie integrujemy je w odniesieniu do aby uzyskać funkcję skumulowanego zagrożenia:t
To daje nam funkcje przetrwania:
Następnie generowania na podstawie pobranych próbek i , zastępując o i przegrupowanie odpowiednie formuły (na podstawie ) W celu symulacji . Powinna to być prosta algebra, którą możesz następnie napisać w języku R, ale daj mi znać poprzez komentarz, jeśli potrzebujesz dalszej pomocy.Xja U∼Uniform(0,1) U S(t) Xi t
źródło