Dieses Skript reproduziert Beispiel 1 der SAS©/STAT-Dokumentation für die EFFECTPLOT-Anweisung. Es analysiert experimentelle Daten (Frankel, 1961) mit dem Ziel, die Ausbeute an Mercaptobenzothiazol (MBT) in Abhängigkeit von Zeit und Temperatur zu maximieren. Der Code verwendet PROC ORTHOREG, um ein Antwortflächenmodell anzupassen und generiert verschiedene Arten von Grafiken (globale Anpassung, Schnittansichten und Beobachtungen mit 'Jitter'-Rauschen), um die Beziehungen zwischen den Variablen zu visualisieren.
Datenanalyse
Type : INTERNE_ERSTELLUNG
Die Daten werden direkt im Skript über den DATA-Schritt 'd' und die DATALINES-Anweisung definiert. Sie enthalten 3 numerische Variablen: Time, Temp und MBT.
1 Codeblock
DATA STEP Data
Erklärung :
Erstellung des Datensatzes 'd', der die Ergebnisse des Experiments enthält. Die Verwendung von ' @@' in der INPUT-Anweisung ermöglicht das Lesen mehrerer Beobachtungen auf derselben Datenzeile.
Erstellung des Datensatzes 'd', der die Ergebnisse des Experiments enthält. Die Verwendung von ' @@' in der INPUT-Anweisung ermöglicht das Lesen mehrerer Beobachtungen auf derselben Datenzeile.
Kopiert!
| 1 | DATA d; |
| 2 | INPUT Time Temp MBT @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json; |
| 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 Codeblock
PROC ORTHOREG
Erklärung :
Anpassung des Regressionsmodells mit PROC ORTHOREG. Die EFFECTPLOT FIT-Anweisung generiert eine Grafik der Modellanpassung mit 'Time' auf der X-Achse, getrennt nach 'Temp'.
Anpassung des Regressionsmodells mit PROC ORTHOREG. Die EFFECTPLOT FIT-Anweisung generiert eine Grafik der Modellanpassung mit 'Time' auf der X-Achse, getrennt nach 'Temp'.
Kopiert!
| 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 Codeblock
PROC ORTHOREG
Erklärung :
Zweite identische Anpassung, jedoch unter Verwendung von EFFECTPLOT SLICEFIT, um Anpassungskurven für spezifische Werte (Slices) der Variable 'Temp' (229, 250, 271, 280) anzuzeigen.
Zweite identische Anpassung, jedoch unter Verwendung von EFFECTPLOT SLICEFIT, um Anpassungskurven für spezifische Werte (Slices) der Variable 'Temp' (229, 250, 271, 280) anzuzeigen.
Kopiert!
| 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 Codeblock
PROC ORTHOREG
Erklärung :
Dritte Anpassung unter Verwendung von EFFECTPLOT als Standard, aber mit der Option OBS(JITTER), um die tatsächlichen Beobachtungen überlagert mit einem leichten Zufallsrauschen (Jitter) anzuzeigen, um visuelle Überlappungen zu vermeiden.
Dritte Anpassung unter Verwendung von EFFECTPLOT als Standard, aber mit der Option OBS(JITTER), um die tatsächlichen Beobachtungen überlagert mit einem leichten Zufallsrauschen (Jitter) anzuzeigen, um visuelle Überlappungen zu vermeiden.
Kopiert!
| 1 | PROC ORTHOREG DATA=d; |
| 2 | model MBT=Time|Time|Temp|Temp; |
| 3 | effectplot / obs(jitter(seed=39393)); |
| 4 | RUN; |
| 5 | ods graphics off; |
Dieses Material wird von We Are Cas "wie besehen" zur Verfügung gestellt. Es gibt keine ausdrücklichen oder stillschweigenden Garantien hinsichtlich der Marktgängigkeit oder Eignung für einen bestimmten Zweck in Bezug auf die hierin enthaltenen Materialien oder Codes. We Are Cas ist nicht verantwortlich für Fehler in diesem Material, wie es jetzt existiert oder existieren wird, noch bietet We Are Cas technischen Support dafür an.
Urheberrechtsinformationen : Frankel (1961) ... From Myers, Response Surface Methodology 1976.
