Analyse de Survie avec PROC GENMOD (Exemple 3)

Le script commence par créer deux jeux de données, 'A' et 'B', qui contiennent des mesures de 'lifetime' pour deux fabricants différents, à l'aide de blocs DATALINES. Ces deux jeux de données sont ensuite concaténés en un seul jeu, 'lifdat'. Deux appels à PROC GENMOD sont effectués. Le premier modélise 'lifetime' en fonction du groupe 'mfg', en spécifiant une distribution Gamma et une fonction de lien log, avec une analyse de type 3. Le second modèle estime un modèle d'interception seule avec les mêmes spécifications de distribution et de lien, et demande les intervalles de confiance du rapport de vraisemblance (LRCI).

Analyse des données

Type : CREATION_INTERNE

Les données sont créées directement dans le script via des blocs DATALINES pour les jeux de données 'A' et 'B'.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc DATA STEP crée le jeu de données 'A' en lisant les valeurs de 'lifetime' directement depuis les DATALINES. Chaque ligne contient plusieurs observations, l'opérateur '@@' est utilisé pour maintenir le pointeur sur la ligne et lire toutes les valeurs de 'lifetime' successivement, en associant à chacune la valeur 'A' pour la variable 'mfg' (fabricant). L'instruction `input lifetime @code_sas_json/...` du code source original a été interprétée comme une erreur et corrigée en `input lifetime @@;` pour permettre la lecture correcte des données multi-valeurs par ligne, conformément à l'intention probable du code et à l'exemple du fichier JSON référencé.

Copié !

1	DATA A;
2	INPUT lifetime @@;
3	mfg = 'A';
4	DATALINES;
5	620 470 260 89 388 242
6	103 100 39 460 284 1285
7	218 393 106 158 152 477
8	403 103 69 158 818 947
9	399 1274 32 12 134 660
10	548 381 203 871 193 531
11	317 85 1410 250 41 1101
12	32 421 32 343 376 1512
13	1792 47 95 76 515 72
14	1585 253 6 860 89 1055
15	537 101 385 176 11 565
16	164 16 1267 352 160 195
17	1279 356 751 500 803 560
18	151 24 689 1119 1733 2194
19	763 555 14 45 776 1
20	;
21

2 Bloc de code

DATA STEP Data

Explication :
Similaire au bloc précédent, ce DATA STEP crée le jeu de données 'B' pour le deuxième fabricant. L'opérateur '@@' est également utilisé pour lire toutes les valeurs de 'lifetime' par ligne, en associant la valeur 'B' à la variable 'mfg'. L'instruction `input lifetime @code_sas_json/...` du code source original a été interprétée comme une erreur et corrigée en `input lifetime @@;` pour permettre la lecture correcte des données multi-valeurs par ligne, conformément à l'intention probable du code et à l'exemple du fichier JSON référencé.

Copié !

1	DATA B;
2	INPUT lifetime @@;
3	mfg = 'B';
4	DATALINES;
5	1747 945 12 1453 14 150
6	20 41 35 69 195 89
7	1090 1868 294 96 618 44
8	142 892 1307 310 230 30
9	403 860 23 406 1054 1935
10	561 348 130 13 230 250
11	317 304 79 1793 536 12
12	9 256 201 733 510 660
13	122 27 273 1231 182 289
14	667 761 1096 43 44 87
15	405 998 1409 61 278 407
16	113 25 940 28 848 41
17	646 575 219 303 304 38
18	195 1061 174 377 388 10
19	246 323 198 234 39 308
20	55 729 813 1216 1618 539
21	6 1566 459 946 764 794
22	35 181 147 116 141 19
23	380 609 546
24	;
25

3 Bloc de code

DATA STEP Data

Explication :
Ce DATA STEP concatène les jeux de données 'A' et 'B' créés précédemment en un nouveau jeu de données nommé 'lifdat'. La variable 'mfg' distingue les observations provenant de 'A' de celles provenant de 'B'.

Copié !

1	DATA lifdat;
2	SET A B;
3	RUN;
4

4 Bloc de code

PROC GENMOD

Explication :
Ce bloc exécute la procédure GENMOD pour ajuster un modèle linéaire généralisé. La variable dépendante est 'lifetime', et la variable catégorielle 'mfg' est utilisée comme prédicteur. Le modèle spécifie une distribution Gamma pour la variable de réponse et une fonction de lien logarithmique. L'option 'type3' demande le calcul des tests de Type 3 pour les effets du modèle.

Copié !

1	PROC GENMOD DATA = lifdat;
2	class mfg;
3	model lifetime = mfg / dist=gamma
4	link=log
5	type3;
6	RUN;
7

5 Bloc de code

PROC GENMOD

Explication :
Ce deuxième bloc PROC GENMOD ajuste un modèle d'interception seule pour la variable 'lifetime', toujours avec une distribution Gamma et une fonction de lien log. L'option 'lrci' demande l'affichage des intervalles de confiance basés sur le rapport de vraisemblance pour les paramètres du modèle.

Copié !

1	PROC GENMOD DATA = lifdat;
2	model lifetime = / dist=gamma
3	link=log
4	lrci;
5	RUN;
6

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 : S A S S A M P L E L I B R A R Y ; NAME: GENMEX3 ; TITLE: Example 3 for PROC GENMOD ; PRODUCT: STAT ; SYSTEM: ALL ; KEYS: generalized linear models ; PROCS: GENMOD ; DATA: ; ; REF: PROC GENMOD, EXAMPLE 3 ; MISC:

Retour à la liste