Ce script reproduit l'exemple 1 de la documentation SAS©/STAT pour l'instruction EFFECTPLOT. Il analyse des données expérimentales (Frankel, 1961) visant à maximiser le rendement du mercaptobenzothiazole (MBT) en fonction du temps et de la température. Le code utilise PROC ORTHOREG pour ajuster un modèle de surface de réponse et génère plusieurs types de graphiques (ajustement global, vues en tranches, et observations avec bruit 'jitter') pour visualiser les relations entre les variables.
Analyse des données
Type : CREATION_INTERNE
Les données sont définies directement dans le script via l'étape DATA 'd' et l'instruction DATALINES. Elles contiennent 3 variables numériques : Time, Temp, et MBT.
1 Bloc de code
DATA STEP Data
Explication :
Création du jeu de données 'd' contenant les résultats de l'expérience. L'utilisation de '@@' dans l'instruction INPUT permet de lire plusieurs observations sur une même ligne de données.
Création du jeu de données 'd' contenant les résultats de l'expérience. L'utilisation de '@@' dans l'instruction INPUT permet de lire plusieurs observations sur une même ligne de données.
Copié !
| 1 | DATA d; |
| 2 | INPUT Time Temp MBT @@; |
| 3 | label Time = "Reaction Time (Hours)" |
| 4 | Temp = "Temperature (Degrees Centigrade)" |
| 5 | MBT = "Percent Yield Mercaptobenzothiazole"; |
| 6 | DATALINES; |
| 7 | 4.0 250 83.8 20.0 250 81.7 12.0 250 82.4 |
| 8 | 12.0 250 82.9 12.0 220 84.7 12.0 280 57.9 |
| 9 | 12.0 250 81.2 6.3 229 81.3 6.3 271 83.1 |
| 10 | 17.7 229 85.3 17.7 271 72.7 4.0 250 82.0 |
| 11 | ; |
2 Bloc de code
PROC ORTHOREG
Explication :
Ajustement du modèle de régression avec PROC ORTHOREG. L'instruction EFFECTPLOT FIT génère un graphique de l'ajustement du modèle avec 'Time' en axe X, séparé par 'Temp'.
Ajustement du modèle de régression avec PROC ORTHOREG. L'instruction EFFECTPLOT FIT génère un graphique de l'ajustement du modèle avec 'Time' en axe X, séparé par 'Temp'.
Copié !
| 1 | ods graphics on; |
| 2 | PROC ORTHOREG DATA=d; |
| 3 | model MBT=Time|Time|Temp|Temp; |
| 4 | effectplot fit(x=time plotby=temp); |
| 5 | RUN; |
3 Bloc de code
PROC ORTHOREG
Explication :
Second ajustement identique, mais utilisant EFFECTPLOT SLICEFIT pour afficher des courbes d'ajustement pour des valeurs spécifiques (tranches) de la variable 'Temp' (229, 250, 271, 280).
Second ajustement identique, mais utilisant EFFECTPLOT SLICEFIT pour afficher des courbes d'ajustement pour des valeurs spécifiques (tranches) de la variable 'Temp' (229, 250, 271, 280).
Copié !
| 1 | PROC ORTHOREG DATA=d; |
| 2 | model MBT=Time|Time|Temp|Temp; |
| 3 | effectplot slicefit(x=time sliceby=temp=229 250 271 280); |
| 4 | RUN; |
4 Bloc de code
PROC ORTHOREG
Explication :
Troisième ajustement utilisant EFFECTPLOT par défaut mais avec l'option OBS(JITTER) pour afficher les observations réelles superposées avec un léger bruit aléatoire (jitter) pour éviter le chevauchement visuel.
Troisième ajustement utilisant EFFECTPLOT par défaut mais avec l'option OBS(JITTER) pour afficher les observations réelles superposées avec un léger bruit aléatoire (jitter) pour éviter le chevauchement visuel.
Copié !
| 1 | PROC ORTHOREG DATA=d; |
| 2 | model MBT=Time|Time|Temp|Temp; |
| 3 | effectplot / obs(jitter(seed=39393)); |
| 4 | RUN; |
| 5 | ods graphics off; |
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 : Frankel (1961) ... From Myers, Response Surface Methodology 1976.
Le Conseil de l'Expert
Stéphanie
Spécialiste Machine Learning et IA.
« L'instruction EFFECTPLOT est universelle : elle fonctionne de la même manière avec PROC GLM, LOGISTIC ou GENMOD. Pour ce cas précis sur le MBT, surveillez particulièrement la courbure à 280°C ; la chute brutale du rendement indique souvent une décomposition thermique que seule une visualisation par tranches permet de bien quantifier pour la production. »