Exemple 11 pour la PROC GLM : Analyse d'un plan de criblage

Le code définit d'abord un jeu de données `HalfFraction` représentant une demi-fraction d'un plan factoriel 2^4. La `PROC GLM` est utilisée pour modéliser la variable de réponse `rate` en fonction des facteurs `power`, `flow`, `pressure`, et `gap`. Ensuite, pour étudier la structure d'aliasing, les facteurs sont recodés en valeurs -1/1 dans un nouveau jeu de données `Coded`. La `PROC GLM` est à nouveau exécutée sur ces données codées avec l'option `aliasing` pour afficher les relations de confusion entre les effets. Par la suite, les données de l'autre moitié de l'expérience sont définies dans `OtherHalf`. Les deux moitiés sont combinées pour créer le jeu de données complet `FullRep`. Finalement, la `PROC GLM` est exécutée sur ce jeu de données complet pour effectuer une analyse de variance complète, sans la confusion présente dans le plan fractionnaire.

Analyse des données

Type : CREATION_INTERNE

Les données sont créées directement dans le script via des étapes DATA avec des cartes de données (datalines). Les jeux de données `HalfFraction` et `OtherHalf` sont créés de cette manière, puis combinés pour former `FullRep`.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc crée le jeu de données `HalfFraction` en lisant les données directement via l'instruction `datalines`. Il représente la première moitié d'une expérience factorielle.

Copié !

1	DATA HalfFraction;
2	INPUT power flow pressure gap rate;
3	DATALINES;
4	0.8 4.5 125 275 550
5	0.8 4.5 200 325 650
6	0.8 550.0 125 325 642
7	0.8 550.0 200 275 601
8	1.2 4.5 125 325 749
9	1.2 4.5 200 275 1052
10	1.2 550.0 125 275 1075
11	1.2 550.0 200 325 729
12	;
13	RUN;

2 Bloc de code

PROC GLM

Explication :
Ce bloc exécute une analyse de modèle linéaire général sur les données `HalfFraction`. Il déclare les variables de classe et spécifie un modèle factoriel complet (interaction d'ordre 4) pour la variable de réponse `rate`.

Copié !

1	PROC GLM DATA=HalfFraction;
2	class power flow pressure gap;
3	model rate=power\|flow\|pressure\|gap;
4	RUN;

3 Bloc de code

DATA STEP Data

Explication :
Ce bloc crée le jeu de données `Coded` à partir de `HalfFraction`. Il transforme les valeurs des facteurs en un codage -1/+1, ce qui est standard pour analyser la structure d'aliasing dans les plans factoriels.

Copié !

1	DATA Coded; SET HalfFraction;
2	power = -1(power =0.80) + 1(power =1.20);
3	flow = -1(flow =4.50) + 1(flow =550 );
4	pressure = -1(pressure=125 ) + 1(pressure=200 );
5	gap = -1(gap =275 ) + 1(gap =325 );
6	RUN;

4 Bloc de code

PROC GLM

Explication :
Ce bloc ré-exécute l'analyse sur les données codées. L'option `solution` demande l'affichage des estimations des paramètres et l'option `aliasing` est utilisée pour examiner la structure de confusion (aliasing) du plan d'expérience.

Copié !

1
2	PROC GLM
3	DATA=Coded;
4	model rate=power\|flow\|pressure\|gap / solution aliasing;
5	RUN;
6

5 Bloc de code

DATA STEP Data

Explication :
Ce bloc crée le jeu de données `OtherHalf` avec les données de la seconde moitié de l'expérience factorielle, en utilisant également des `datalines`.

Copié !

1	DATA OtherHalf;
2	INPUT power flow pressure gap rate;
3	DATALINES;
4	0.8 4.5 125 325 669
5	0.8 4.5 200 275 604
6	0.8 550.0 125 275 633
7	0.8 550.0 200 325 635
8	1.2 4.5 125 275 1037
9	1.2 4.5 200 325 868
10	1.2 550.0 125 325 860
11	1.2 550.0 200 275 1063
12	;
13	RUN;

6 Bloc de code

DATA STEP Data

Explication :
Ce bloc combine les jeux de données `HalfFraction` et `OtherHalf` pour créer un jeu de données complet `FullRep` qui représente l'expérience factorielle entière.

Copié !

1	DATA FullRep;
2	SET HalfFraction OtherHalf;
3	RUN;

7 Bloc de code

PROC GLM

Explication :
Ce bloc effectue l'analyse de variance finale sur le plan d'expérience complet (`FullRep`), permettant une estimation non biaisée des effets principaux et des interactions.

Copié !

1	PROC GLM DATA=FullRep;
2	class power flow pressure gap;
3	model rate=power\|flow\|pressure\|gap;
4	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 : S A S S A M P L E L I B R A R Y

Retour à la liste

1	DATA Coded; SET HalfFraction;
2	power = -1(power =0.80) + 1(power =1.20);
3	flow = -1(flow =4.50) + 1(flow =550 );
4	pressure = -1(pressure=125 ) + 1(pressure=200 );
5	gap = -1(gap =275 ) + 1(gap =325 );
6	RUN;