Les données sont générées statiquement dans le script via l'étape DATA 'multicenter' et l'instruction DATALINES.
1 Bloc de code
DATA STEP Data
Explication : Création du jeu de données 'multicenter' contenant les résultats de l'étude : le centre, le groupe de traitement, le nombre total de patients (n) et le nombre d'effets secondaires observés.
Copié !
data multicenter;
input center group$ n sideeffect;
datalines;
1 A 32 14
1 B 33 18
2 A 30 4
2 B 28 8
3 A 23 14
3 B 24 9
4 A 22 7
4 B 22 10
5 A 20 6
5 B 21 12
6 A 19 1
6 B 20 3
7 A 17 2
7 B 17 6
8 A 16 7
8 B 15 9
9 A 13 1
9 B 14 5
10 A 13 3
10 B 13 1
11 A 11 1
11 B 12 2
12 A 10 1
12 B 9 0
13 A 9 2
13 B 9 6
14 A 8 1
14 B 8 1
15 A 7 1
15 B 8 0
;
1
DATA multicenter;
2
INPUT center group$ n sideeffect;
3
DATALINES;
4
1 A 3214
5
1 B 3318
6
2 A 304
7
2 B 288
8
3 A 2314
9
3 B 249
10
4 A 227
11
4 B 2210
12
5 A 206
13
5 B 2112
14
6 A 191
15
6 B 203
16
7 A 172
17
7 B 176
18
8 A 167
19
8 B 159
20
9 A 131
21
9 B 145
22
10 A 133
23
10 B 131
24
11 A 111
25
11 B 122
26
12 A 101
27
12 B 9 0
28
13 A 92
29
13 B 96
30
14 A 81
31
14 B 81
32
15 A 71
33
15 B 8 0
34
;
2 Bloc de code
PROC GLIMMIX
Explication : Exécution de la procédure GLIMMIX. Le modèle spécifie une distribution binomiale (syntaxe events/trials) avec 'group' comme effet fixe. Un intercept aléatoire est ajouté pour chaque 'center' afin de capturer la corrélation intra-centre.
Copié !
proc glimmix data=multicenter;
class center group;
model sideeffect/n = group / solution;
random intercept / subject=center;
run;
1
PROC GLIMMIXDATA=multicenter;
2
class center group;
3
model sideeffect/n = group / solution;
4
random intercept / subject=center;
5
RUN;
3 Bloc de code
PROC GLIMMIX
Explication : Seconde exécution identique à la première, mais focalisée sur l'affichage des moyennes des moindres carrés (LS-means). L'option 'ilink' est utilisée pour rapporter les estimations sur l'échelle des données originales (probabilités) plutôt que sur l'échelle du lien (logit).
Copié !
ods select lsmeans;
proc glimmix data=multicenter;
class center group;
model sideeffect/n = group / solution;
random intercept / subject=center;
lsmeans group / cl ilink;
run;
1
ods select lsmeans;
2
PROC GLIMMIXDATA=multicenter;
3
class center group;
4
model sideeffect/n = group / solution;
5
random intercept / subject=center;
6
lsmeans group / cl ilink;
7
RUN;
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: gmxgs1 / PRODUCT: STAT
SAS et tous les autres noms de produits ou de services de SAS Institute Inc. sont des marques déposées ou des marques de commerce de SAS Institute Inc. aux États-Unis et dans d'autres pays. ® indique un enregistrement aux États-Unis. WeAreCAS est un site communautaire indépendant et n'est pas affilié à SAS Institute Inc.
Ce site utilise des cookies techniques et analytiques pour améliorer votre expérience.
En savoir plus.