Natural Recovery Time - Exemple PROC LIFEREG

Ce script illustre l'utilisation de la procédure LIFEREG pour analyser des données de survie avec censure arbitraire (intervalles). Il modélise la relation entre le temps de récupération, l'âge et le sexe, incluant une interaction. Deux analyses sont effectuées pour générer des graphiques de probabilité (probplots) basés sur des données de référence spécifiques (xrow1 et xrow2).

Analyse des données

Type : CREATION_INTERNE

Toutes les données (tables 'mice', 'xrow1', 'xrow2') sont définies directement dans le script via des blocs DATALINES.

1 Bloc de code

DATA STEP Data

Explication :
Création du jeu de données principal 'mice' contenant les variables explicatives (sexe, age) et les bornes de l'intervalle de temps (time1, time2).

Copié !

1	title 'Natural Recovery Time';
2	DATA mice;
3	INPUT sex age time1 time2;
4	DATALINES;
5	1 57 631 631
6	1 45 . 170
7	1 54 227 227
8	1 43 143 143
9	1 64 916 .
10	1 67 691 705
11	1 44 100 100
12	1 59 730 .
13	1 47 365 365
14	1 74 1916 1916
15	2 79 1326 .
16	2 75 837 837
17	2 84 1200 1235
18	2 54 . 365
19	2 74 1255 1255
20	2 71 1823 .
21	2 65 537 637
22	2 33 583 683
23	2 77 955 .
24	2 46 577 577
25	;

2 Bloc de code

DATA STEP Data

Explication :
Création d'un jeu de données de référence 'xrow1' pour la prédiction ou le traçage (Sexe 1, Age 50).

Copié !

1	DATA xrow1;
2	INPUT sex age time1 time2;
3	DATALINES;
4	1 50 . .
5	;

3 Bloc de code

DATA STEP Data

Explication :
Création d'un second jeu de données de référence 'xrow2' (Sexe 2, Age 60.6).

Copié !

1	DATA xrow2;
2	INPUT sex age time1 time2;
3	DATALINES;
4	2 60.6 . .
5	;

4 Bloc de code

PROC LIFEREG

Explication :
Exécution du modèle de survie Weibull avec les données 'mice'. L'option 'xdata=xrow1' est utilisée pour personnaliser le graphique de probabilité. L'instruction 'model' définit les bornes de censure (time1, time2) et les covariables.

Copié !

1	ods graphics on;
2	PROC LIFEREG DATA=mice xdata=xrow1;
3	class sex;
4	model (time1, time2) = age sex age*sex / dist=Weibull;
5
6	probplot / nodata
7	plower=.5
8	vref(intersect) = 75
9	vreflab = '75 Percent';
10	inset;
11	RUN;

5 Bloc de code

PROC LIFEREG

Explication :
Seconde exécution du modèle identique au premier, mais en utilisant 'xrow2' comme données de référence pour le graphique.

Copié !

1	PROC LIFEREG DATA=mice xdata=xrow2;
2	class sex;
3	model (time1, time2) = age sex age*sex / dist=Weibull;
4
5	probplot / nodata
6	plower=.5
7	vref(intersect) = 75
8	vreflab = '75 Percent';
9	inset;
10	RUN;
11	title;

Ce matériel est fourni "tel quel" par We Are Cas. Il n'y a aucune garantie, expresse ou implicite, quant à la qualité marchande ou à l'adéquation à un usage particulier concernant le matériel ou le code contenu dans les présentes. We Are Cas n'est pas responsable des erreurs dans ce matériel tel qu'il existe maintenant ou existera, et We Are Cas ne fournit pas de support technique pour celui-ci.

Informations de Copyright : SAS SAMPLE LIBRARY - NAME: LIFEREX4

Retour à la liste