Publié le :
Statistique CREATION_INTERNE

Contrastes isotoniques avec PROC MIXED

Ce code est également disponible en : Deutsch English Español
En attente de validation
Le script commence par un DATA STEP qui génère le dataset 'FerriteCores'. Ce dataset simule des mesures de 'MagneticForce' pour quatre niveaux de 'Temp' (température), avec cinq répétitions par niveau, les données étant fournies via des DATALINES. Ensuite, PROC MIXED est invoquée pour effectuer une analyse de variance sur 'MagneticForce' en fonction de 'Temp'. La déclaration CLASS spécifie 'Temp' comme variable catégorielle. La déclaration LSMESTIMATE est utilisée pour définir et tester des contrastes isotoniques spécifiques, permettant de comparer les moyennes des niveaux de 'Temp' de manière ordonnée. Les options '/ adjust=simulate(seed=1) cl upper' sont appliquées pour l'ajustement des p-valeurs, le calcul des intervalles de confiance et la définition d'une borne supérieure. Enfin, ODS SELECT LSMestimates est utilisé pour afficher uniquement les résultats pertinents des estimations des moindres carrés.
Analyse des données

Type : CREATION_INTERNE


Les données sont générées directement dans le script via un DATA STEP et des DATALINES pour le dataset 'FerriteCores'. Le champ 'MagneticForce' est lu à partir des DATALINES.

1 Bloc de code
DATA STEP Data
Explication :
Ce bloc DATA STEP crée le dataset 'FerriteCores'. Il initialise une boucle pour générer une variable 'Temp' (température) de 1 à 4, et une boucle interne pour 5 répétitions par 'Temp'. La variable 'MagneticForce' est lue à partir des DATALINES, simulant des mesures pour différentes conditions de température.
Copié !
1DATA FerriteCores;
2 DO Temp = 1 to 4;
3 DO rep = 1 to 5; drop rep;
4 INPUT MagneticForce;
5 OUTPUT;
6 END;
7 END;
8 DATALINES;
910.8 9.9 10.7 10.4 9.7
1010.7 10.6 11.0 10.8 10.9
1111.9 11.2 11.0 11.1 11.3
1211.4 10.7 10.9 11.3 11.7
13;
14 
2 Bloc de code
PROC MIXED
Explication :
Ce bloc utilise PROC MIXED pour analyser le dataset 'FerriteCores'. 'Temp' est définie comme une variable de classification. Un modèle est spécifié avec 'MagneticForce' comme réponse et 'Temp' comme prédicteur. La déclaration LSMESTIMATE est utilisée pour tester des contrastes isotoniques, évaluant des hypothèses ordonnées sur les moyennes des niveaux de 'Temp'. Les options incluent l'ajustement des p-valeurs par simulation, le calcul des intervalles de confiance et la spécification d'une borne supérieure. ODS SELECT LSMestimates limite la sortie aux résultats des estimations des moindres carrés.
Copié !
1PROC MIXED DATA=FerriteCores;
2 class Temp;
3 model MagneticForce = Temp;
4 lsmestimate Temp
5 'avg(1:1)<avg(2:4)' -3 1 1 1 divisor=3,
6 'avg(1:2)<avg(3:4)' -1 -1 1 1 divisor=2,
7 'avg(1:3)<avg(4:4)' -1 -1 -1 3 divisor=3
8 / adjust=simulate(seed=1) cl upper;
9 ods select LSMestimates;
10RUN;
11 
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 : S A S S A M P L E L I B R A R Y