Analyse de survie de base R Flashcards

1
Q

Une fonction clé pour l’analyse des données de survie dans R est la fonction

A

Surv()

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
2
Q

role de Surv()

A

spécifier le type de données de survie que nous avons, à savoir,
censuré à droite, censuré à gauche, censuré par intervalle.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
3
Q

la fonction Surv()accepte comme arguments

A

premier argument les temps de survie observés,

et comme second l’indicateur d’événement.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
4
Q

La fonction utilisée pour calculer l’estimation de Kaplan-Meier est

A

survfit()

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
5
Q

exemple de estimation de kaplan meir

A

fit = survfit(Surv(time, status) ~ 1, data = )

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
6
Q

La sortie de base de la fonction survfit()donne

A

le nombre de sujets,
le nombre d’événements,
le temps de survie médian
et son intervalle de confiance à 95%.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
7
Q

le code pour calculer les probabilités de survie à 1000 et 2000 jours

A

summary(KM_fit, times = c(1000, 2000))

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
8
Q

voulons trouver à combien de jours la probabilité de survie est égale à 0,7 et à combien de jours elle est égale à 0,6 - le code est

A

quantile(KM_fit, probs = 1 - c(0.7, 0.6))

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
9
Q

Pour calculer l’estimateur de Breslow, nous utilisons

A

la fonction survfit() avec la même syntaxe que pour l’estimateur de Kaplan-Meier

la seule différence est que nous devons maintenant définir l’argument typesur “fleming-harrington”.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
10
Q

Nous avons deux tests statistiques standard pour tester une hypothèse simple avec une variable de regroupement

A

logrank

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
11
Q

quel test de comparaison choisir

A

Le test du log-rank est le test le plus puissant lorsque l’hypothèse des risques proportionnels (PH) est satisfaite.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
12
Q

Pour vérifier l’hypothèse l’hypothèse de risque proportionnel

A

nous pouvons tracer les fonctions de risque cumulées pour les deux groupes;
lorsque PH est satisfait, les deux courbes seront proportionnelles l’une à l’autre (c’est-à-dire qu’elles s’éloignent progressivement l’une de l’autre).

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
13
Q

formule pour tracer la fonction de risque cumulé

A

plot(KM_fit, fun = “cumhaz”)

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
14
Q

les differente fonctions daln l’analue de survie

A
La fonction de survie S(t)
La fonction de probabilité de densité f (t)
La fonction de répartition F(t)
La fonction de risque h(t)
La fonction de risque cumul´ee H(t)
How well did you know this?
1
Not at all
2
3
4
5
Perfectly
15
Q

La fonction de probabilit´e de densit´e f (t)

A

repr´esente la limite de

la probabilit´e que l’´ev´enement se produise au temps t

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
16
Q

La fonction de survie S(t)

A

probabilité que l’evenement se produise apres le temps t

S(t) = Pr(T > t)

17
Q

La fonction de répartition F(t)

A

probabilité que l’évènement se produise après le temps t.

F(t) = Pr(T ≤ t) = 1 − S(t)

18
Q

La fonction de risque h(t)

A

représente la limite de la probabilité que l’´évènement se produise au temps t, sachant qu’aucun évènement ne s’est produit jusqu’au temps t.

19
Q

La fonction de risque cumulée H(t)

A

Ht = -log(S(t))

20
Q

Formulation générale d’un modelé `a risques proportionnels (PH)

A

λ(t|Xi) = λ0(t)exp(βi*Xi)

21
Q

λ0(t)

A

est la fonction de risque de base (groupe de référence :x1 = x2 = ::: = 0).

22
Q

estimation de λ0(t)

A

λ0(t) n’est pas estimée.

23
Q

Pourquoi le Modelé et dit semi-paramétrique

A

ˆ Partie λ0(t) est non-paramétrique (pas d’estimation)
ˆ Partie exp(βi*Xi) et paramétrique (fonction Exponentielle avec prédicteur
linéaire)

24
Q

La proportionnalité des risques

A

βk (k = 1; :::; p) ne dependent pas du temps

25
Q

pour vérifier graphiquement PH

A

Tracer log(−log(S(t)))

26
Q

La proportionnalité des risques est respectée

A

si les courbes sont parallèles

27
Q

La proportionnalité des risques

Avantages de l’approche visuelle

A

Interprétation des résultats.
Évaluation possible de la correction pour gérer la
non-proportionnalité des risques

28
Q

La proportionnalité des risques

Limites de l’approche visuelle

A

ˆ Interprétations subjectives
ˆ Seulement une covariable
ˆ Covariables catégorielles uniquement

29
Q

modéliser la non-proportionnalité du risque lié à la

covariable X1

A

λ(t| X1; :::; Xi) = λ0(t)exp( β1(t) X1 + β2X2 + ::: + βiXi)

β1(t) X1 = La relation entre la Covariable change avec le temps.

on dit que la valeur est dépendante du temp

30
Q

modeliser un effet non-proportionnel de la covariable
X1 ?
du model λ(t| X1; :::; Xi) = λ0(t)exp( β1(t) X1 + β2X2 + … βiXi)

A

λ(t| X1; :::; Xi) = λ0(t)exp( β1(t) X1 + γ1X1 + β2X2 + …+ βiXi)
Si t ≤ u alors HR = exp(β1).
Si t > u alors HR = exp(β1 + γ1).
Si γ1 = 0, alors il n’y a aucun effet dépendant du temps.
on dit que l’effet et dépendant du temps

31
Q

la démarché dans la prise en compte de variable temporelle

A

modifier la structure de la base de données

32
Q

exemple d’une indication de stratification

A

une variable catégorielle qui ne satisfait pas l’hypothèse des risques proportionnels,

33
Q

Pour ajuster le modèle de Cox stratifié,

A

inclure le facteur de stratification dans la formule du modèle utilisant la fonction strata()

34
Q

commentaire du code suivant:
coxph(Surv(years, status2) ~ age + drug + log(serBilir) + strata(sex), data = pbc2.id)

et caractéristique de la stratification

A

correspond à un modèle de Cox stratifié pour l’ensemble de données PBC avec différentes fonctions de risque de base pour chacun sex

Notez qu’une caractéristique de la stratification est que nous corrigeons (de la manière la plus générale) l’analyse pour sex mais nous n’obtenons aucun coefficient pour sex.

35
Q

Si nous voulons permettre que l’effet de age sur le logarithme diffère entre les hommes et les femmes,

coxph(Surv(years, status2) ~ age + drug + log(serBilir),
data = pbc2.id)

A

nous incluons le terme d’interaction entre age et strata(sex)

coxph(Surv(years, status2) ~ age * strata(sex) + drug + log(serBilir),
data = pbc2.id)

36
Q

Lorsque nous calculons les probabilités de survie à partir d’un modèle de Cox stratifié, nous obtenons

A

une courbe différente par strate
Nous observons que même si nous avons précisé que le patient pour lequel nous voulons calculer les probabilités de survie est un homme, nous obtenons deux courbes à la fois pour les hommes et les femmes. C’est parce qu’il sexs’agit d’un facteur de stratification et non d’une covariable.

37
Q

Il est possible de stratifier selon plusieurs variables, comme le sexe et la présence d’œdèmes. Le code à renseigner est
Dans ce cas,

A

alors strata(sex, edema).

Dans ce cas, les différents groupes sont constitués par toutes les combinaisons possibles des
différentes variables., les différents groupes sont constitués par toutes les combinaisons possibles des
différentes variables.