Discussion :

[yo-stat]

Bonjour,

Je travaille actuellement sur un modèle de durée de vie paramétrique à partir de la procédure LIFEREG avec pour un 1er test la fonction Weibull.

J'aimerai recalculer la fonction de survie à partir des paramètres sortis par la procédure LIFEREG. Mais je n'arrive pas à retomber sur la fonction de survie. Je sais que la durée est transformée en entrée mais je ne m'y retrouve pas.

Savez-vous comment interpréter les paramètres de la fonction Weibull (exemple ci-dessous) pour reconstituer la fonction de survie :

Parameter Estimate
Intercept 6.9251
Scale 0.7404
Weibull Scale 1017.541
Weibull Shape 1.3507

[fafabzh6]

Bonjour,
Peux-tu nous montrer le code SAS que tu as utilisé pour reconstituer ta survie à l'aide du modèle de Weibull ?
Merci

[yo-stat]

Voici mon code,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
proc lifereg data=BASE ;
  model duree*censure(1) =  / dist=weibull;
  output out=out_test cdf=Proba predicted=Predtime quantiles=0.5 ;
run ;

Je n'ai toujours pas débloqué mon problème.
En gros, en faisant 1-Proba de la table out_test, j'obtiens la survie mais uniquement pour les durées observées. Pour extrapoler les résultats, je veux exploiter les paramètres de la fonction weibull pour recalculer la fonction de survie avec la formule S(t)=exp(-(B.t)**a)
sachant que la procédure lifereg travaille sur le log.

Merci,

[yo-stat]

J'ai réussi à recalculer la fonction de survie avec les paramètres sortis par la porc lifereg :

survie_weibull= exp(-exp(-intercept/scale)*(duree)** Weibull_shape);

A présent, il faut que je sois capable de calculer les variables d'intérêt, à savoir la moyenne de durée de vie pour ma population ainsi que l'espérance de vie résiduelle qui correspond à un calcul d'intégral.

ça me semble compliqué à implémenter sous SAS...

[fafabzh6]

Avec SAS IML il me semble que tu peux calculer une intégrale !
Peut-être que ce message pourra aussi t'aider !
Bon courage