Publicado el :
Estadística CREATION_INTERNE

Medidas repetidas multipoblación con PROC CATMOD

Este código también está disponible en: Deutsch English Français
En espera de validación
Este script ilustra el uso del procedimiento CATMOD para analizar respuestas binarias (0 o 1) registradas durante cuatro ensayos sucesivos en sujetos de tres grupos diferentes. Realiza un análisis de medidas repetidas ajustando primero un modelo saturado sobre las probabilidades marginales y luego un modelo con efectos anidados.
Análisis de datos

Type : CREATION_INTERNE


Los datos se generan directamente en el script a través de la instrucción DATALINES. La instrucción INPUT utiliza el operador ' @@' para leer varias observaciones por línea.

1 Bloque de código
DATA STEP Data
Explicación :
Creación del conjunto de datos 'group'. Las variables a, b, c, d representan las respuestas a los 4 ensayos. 'Group' es el identificador del grupo y 'wt' es la ponderación (frecuencia). El operador ' @@' permite leer varias observaciones en una sola línea física.
¡Copiado!
1DATA group;
2 INPUT a b c d Group wt @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json;
3 DATALINES;
41 1 1 1 2 2 0 0 0 0 2 2 0 0 1 0 1 2 0 0 1 0 2 2
50 0 0 1 1 4 0 0 0 1 2 1 0 0 0 1 3 3 1 0 0 1 2 1
60 0 1 1 1 1 0 0 1 1 2 2 0 0 1 1 3 5 0 1 0 0 1 4
70 1 0 0 2 1 0 1 0 1 2 1 0 1 0 1 3 2 0 1 1 0 3 1
81 0 0 0 1 3 1 0 0 0 2 1 0 1 1 1 2 1 0 1 1 1 3 2
91 0 1 0 1 1 1 0 1 1 2 1 1 0 1 1 3 2
10;
2 Bloque de código
PROC CATMOD
Explicación :
Ejecución de PROC CATMOD. La primera parte ajusta un modelo saturado para medidas repetidas (definidas por 'repeated Trial 4') sobre las probabilidades marginales. La segunda instrucción MODEL (modo interactivo) prueba una estructura específica donde el efecto de los ensayos está anidado dentro del grupo 3.
¡Copiado!
1title 'Multiple-Population Repeated Measures';
2PROC CATMOD DATA=group;
3 weight wt;
4 response marginals;
5 model a*b*c*d=Group _response_ Group*_response_
6 / freq;
7 repeated Trial 4;
8 title2 'Saturated Model';
9RUN;
10 
11 model a*b*c*d=Group _response_(Group=3)
12 / noprofile noparm design;
13 title2 'Trial Nested within Group 3';
14QUIT;
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, SAS/STAT User's Guide