Tests Freeman-Tukey et t avec rééchantillonnage Bootstrap (PROC MULTTEST)

Ce script analyse des données de carcinogénicité animale (présence de tumeurs S1 et S2, temps de décès T) en fonction de la dose administrée. Il utilise la PROC MULTTEST avec l'option BOOTSTRAP pour calculer des valeurs de p ajustées, en appliquant un test de Freeman-Tukey pour les variables binaires et un test t pour la variable continue.

Analyse des données

Type : CREATION_INTERNE

Les données sont générées via un DATA STEP utilisant DATALINES avec l'opérateur '@@' pour lire plusieurs observations par ligne.

1 Bloc de code

DATA STEP Data

Explication :
Création du jeu de données 'a'. Chaque observation représente un animal avec deux indicateurs de tumeur (S1, S2), le temps de décès (T) et la dose reçue. L'opérateur '@@' est utilisé pour lire plusieurs enregistrements sur une seule ligne de données.

Copié !

1	title 'Freeman-Tukey and t Tests with Bootstrap Resampling';
2
3	DATA a;
4	INPUT S1 S2 T Dose @@;
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 Bloc de code

PROC MULTTEST Data

Explication :
Exécution des tests statistiques multiples. L'option 'bootstrap' active le rééchantillonnage (10 000 échantillons, graine fixée). On spécifie un test de Freeman-Tukey (ft) pour les variables binaires S1/S2 et un test de moyenne (mean) pour T. Un contraste linéaire est défini pour tester la tendance à travers les niveaux de dose.

Copié !

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 Bloc de code

PROC PRINT

Explication :
Affichage des 36 premières observations du jeu de données de sortie 'res' contenant les résultats du rééchantillonnage.

Copié !

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

Retour à la liste