Ajustement de p-valeurs brutes avec PROC MULTTEST

Le script commence par créer une table de données nommée 'a' qui contient les p-valeurs brutes. Chaque observation correspond à un test statistique identifié par une variable 'Test'. Ensuite, la procédure PROC MULTTEST est appelée avec cette table en entrée (option `inpvalues=a`). Elle calcule et affiche les p-valeurs ajustées selon trois méthodes différentes demandées : Holm, Hochberg, et FDR. C'est un exemple typique de gestion du problème des comparaisons multiples en statistique.

Analyse des données

Type : CREATION_INTERNE

Les données (p-valeurs brutes) sont directement intégrées dans le code via une instruction DATALINES, créant la table 'a'.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc DATA STEP crée la table de travail 'a'. Il lit les données intégrées dans le code via `datalines`. Chaque ligne contient un identifiant de test (Test) et la p-valeur brute associée (Raw_P).

Copié !

1	title 'Inputting Raw p-Values';
2
3	DATA a;
4	INPUT Test$ Raw_P;
5	DATALINES;
6	test01 0.28282 test02 0.30688 test03 0.71022
7	test04 0.77175 test05 0.78180 test06 0.88581
8	test07 0.54685 test08 0.84978 test09 0.24228
9	test10 0.58977 test11 0.03498 test12 0.41607
10	test13 0.31631 test14 0.05254 test15 0.45061
11	test16 0.75758 test17 0.12496 test18 0.49485
12	test19 0.21572 test20 0.50505 test21 0.94372
13	test22 0.81260 test23 0.77596 test24 0.36889
14	;
15

2 Bloc de code

PROC MULTTEST

Explication :
Cette procédure statistique prend en entrée la table 'a' contenant les p-valeurs. Elle applique plusieurs méthodes de correction pour les tests multiples : Holm (holm), Hochberg (hoc), et le contrôle du taux de fausses découvertes (fdr). Les résultats avec les p-valeurs ajustées sont affichés dans la sortie standard.

Copié !

1	PROC MULTTEST inpvalues=a holm hoc fdr;
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. REF: SAS/STAT User's Guide, PROC MULTTEST chapter

Retour à la liste