Pruebas de Freeman-Tukey y t con remuestreo Bootstrap (PROC MULTTEST)

Este script analiza datos de carcinogenicidad animal (presencia de tumores S1 y S2, tiempo de fallecimiento T) en función de la dosis administrada. Utiliza PROC MULTTEST con la opción BOOTSTRAP para calcular valores p ajustados, aplicando una prueba de Freeman-Tukey para variables binarias y una prueba t para la variable continua.

Análisis de datos

Type : CREATION_INTERNE

Los datos se generan mediante un DATA STEP que utiliza DATALINES con 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 'a'. Cada observación representa un animal con dos indicadores de tumor (S1, S2), el tiempo de fallecimiento (T) y la dosis recibida. El operador ' @@' se utiliza para leer múltiples registros en una sola línea de datos.

¡Copiado!

1	title 'Freeman-Tukey and t Tests with Bootstrap Resampling';
2
3	DATA a;
4	INPUT S1 S2 T Dose @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json;
5	DATALINES;
6	0 1 104 1 1 0 80 1 0 1 104 1
7	0 1 104 1 0 1 100 1 1 0 104 1
8	1 0 85 2 1 0 60 2 0 1 89 2
9	1 0 96 2 0 1 96 2 1 0 99 2
10	1 0 60 3 1 0 50 3 1 0 80 3
11	0 1 98 3 0 1 99 3 1 0 50 3
12	;

2 Bloque de código

PROC MULTTEST Data

Explicación :
Ejecución de las pruebas estadísticas múltiples. La opción 'bootstrap' activa el remuestreo (10 000 muestras, semilla fija). Se especifica una prueba de Freeman-Tukey (ft) para las variables binarias S1/S2 y una prueba de media (mean) para T. Se define un contraste lineal para probar la tendencia a través de los niveles de dosis.

¡Copiado!

1	PROC MULTTEST DATA=a bootstrap nsample=10000 seed=37081 outsamp=res;
2	test ft(S1 S2 / lowertailed) mean(T / lowertailed);
3	class Dose;
4	contrast 'Linear Trend' 0 1 2;
5	RUN;

3 Bloque de código

PROC PRINT

Explicación :
Visualización de las primeras 36 observaciones del conjunto de datos de salida 'res' que contiene los resultados del remuestreo.

¡Copiado!

1	PROC PRINT DATA=res(obs=36);
2	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 - Example provided by Keith Soper, Merck

Retour à la liste