Publicado el :
Estadística CREATION_INTERNE

Prueba de Peto sobre Mortalidad y Prevalencia con PROC MULTTEST

Este código también está disponible en: English Français
En espera de validación
Este script crea un conjunto de datos que simula observaciones de animales (tipos de tumores, dosis, tiempo de fallecimiento). Deriva una variable de estratificación (B) basada en el tiempo de fallecimiento (T). Luego, se utiliza el procedimiento MULTTEST para realizar una prueba de Peto unilateral (cola superior) con 20 permutaciones, estratificando por B y analizando el efecto de la dosis a través de un contraste lineal.
Análisis de datos

Type : CREATION_INTERNE


Los datos se definen directamente en el script a través de DATALINES (tabla 'a') y utilizan el formato de lectura continua ( @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json).

1 Bloque de código
DATA STEP Data
Explicación :
Creación de la tabla 'a' a partir de datos internos. La instrucción INPUT utiliza '@@' para leer varias observaciones en la misma línea física. Se calcula una variable de estratificación 'B' condicionalmente en función del valor de 'T'.
¡Copiado!
1DATA a;
2 INPUT S1-S3 T Dose @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json;
3 IF T<=90 THEN B=1; ELSE B=2;
4 DATALINES;
50 0 0 104 0 2 0 1 80 0 0 0 1 104 0
60 0 0 104 0 0 2 0 100 0 1 0 0 104 0
72 0 0 85 1 2 1 0 60 1 0 1 0 89 1
82 0 1 96 1 0 0 0 96 1 2 0 1 99 1
92 1 1 60 2 2 0 0 50 2 2 0 1 80 2
100 0 2 98 2 0 0 1 99 2 2 1 1 50 2
11;
2 Bloque de código
PROC MULTTEST Data
Explicación :
Ejecución del procedimiento MULTTEST para calcular las probabilidades ajustadas. La prueba 'peto' se especifica para las variables S1 a S3, teniendo en cuenta el tiempo T y utilizando permutaciones. La opción 'out=p' guarda los resultados.
¡Copiado!
1PROC MULTTEST DATA=a notables out=p stepsid;
2 test peto(S1-S3 / permutation=20 time=T uppertailed);
3 class Dose;
4 strata B;
5 contrast 'mort-prev' 0 1 2;
6RUN;
3 Bloque de código
PROC PRINT
Explicación :
Visualización del contenido de la tabla de salida 'p' que contiene los resultados estadísticos (valores p).
¡Copiado!
1PROC PRINT DATA=p;
2RUN;
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: MULTEX3