Los datos se generan estáticamente en el script a través del paso DATA 'multicenter' y la instrucción DATALINES.
1 Bloque de código
DATA STEP Data
Explicación : Creación del conjunto de datos 'multicenter' que contiene los resultados del estudio: el centro, el grupo de tratamiento, el número total de pacientes (n) y el número de efectos secundarios observados.
¡Copiado!
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 Bloque de código
PROC GLIMMIX
Explicación : Ejecución del procedimiento GLIMMIX. El modelo especifica una distribución binomial (sintaxis events/trials) con 'group' como efecto fijo. Se agrega un intercepto aleatorio para cada 'center' para capturar la correlación intracentro.
¡Copiado!
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 Bloque de código
PROC GLIMMIX
Explicación : Segunda ejecución idéntica a la primera, pero enfocada en mostrar las medias de mínimos cuadrados (LS-means). La opción 'ilink' se utiliza para reportar las estimaciones en la escala de los datos originales (probabilidades) en lugar de en la escala del enlace (logit).
¡Copiado!
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;
Este material se proporciona "tal cual" por We Are Cas. No hay garantías, expresas o implícitas, en cuanto a la comerciabilidad o idoneidad para un propósito particular con respecto a los materiales o el código contenidos en este documento. We Are Cas no es responsable de los errores en este material tal como existe ahora o existirá, ni We Are Cas proporciona soporte técnico para el mismo.
Información de copyright : SAS SAMPLE LIBRARY / NAME: gmxgs1 / PRODUCT: STAT
SAS y todos los demás nombres de productos o servicios de SAS Institute Inc. son marcas registradas o marcas comerciales de SAS Institute Inc. en los EE. UU. y otros países. ® indica registro en los EE. UU. WeAreCAS es un sitio comunitario independiente y no está afiliado a SAS Institute Inc.
Este sitio utiliza cookies técnicas y analíticas para mejorar su experiencia.
Saber más.