Este script reproduce el ejemplo 1 de la documentación de SAS©/STAT para la instrucción EFFECTPLOT. Analiza datos experimentales (Frankel, 1961) con el objetivo de maximizar el rendimiento del mercaptobenzotiazol (MBT) en función del tiempo y la temperatura. El código utiliza PROC ORTHOREG para ajustar un modelo de superficie de respuesta y genera varios tipos de gráficos (ajuste global, vistas en secciones y observaciones con ruido 'jitter') para visualizar las relaciones entre las variables.
Análisis de datos
Type : CREATION_INTERNE
Los datos se definen directamente en el script a través del paso DATA 'd' y la instrucción DATALINES. Contienen 3 variables numéricas: Time, Temp y MBT.
1 Bloque de código
DATA STEP Data
Explicación :
Creación del conjunto de datos 'd' que contiene los resultados del experimento. El uso de ' @@' en la instrucción INPUT permite leer múltiples observaciones en una misma línea de datos.
Creación del conjunto de datos 'd' que contiene los resultados del experimento. El uso de ' @@' en la instrucción INPUT permite leer múltiples observaciones en una misma línea de datos.
¡Copiado!
| 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 Bloque de código
PROC ORTHOREG
Explicación :
Ajuste del modelo de regresión con PROC ORTHOREG. La instrucción EFFECTPLOT FIT genera un gráfico del ajuste del modelo con 'Time' en el eje X, separado por 'Temp'.
Ajuste del modelo de regresión con PROC ORTHOREG. La instrucción EFFECTPLOT FIT genera un gráfico del ajuste del modelo con 'Time' en el eje X, separado por 'Temp'.
¡Copiado!
| 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 Bloque de código
PROC ORTHOREG
Explicación :
Segundo ajuste idéntico, pero utilizando EFFECTPLOT SLICEFIT para mostrar curvas de ajuste para valores específicos (secciones) de la variable 'Temp' (229, 250, 271, 280).
Segundo ajuste idéntico, pero utilizando EFFECTPLOT SLICEFIT para mostrar curvas de ajuste para valores específicos (secciones) de la variable 'Temp' (229, 250, 271, 280).
¡Copiado!
| 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 Bloque de código
PROC ORTHOREG
Explicación :
Tercer ajuste utilizando EFFECTPLOT por defecto pero con la opción OBS(JITTER) para mostrar las observaciones reales superpuestas con un ligero ruido aleatorio (jitter) para evitar la superposición visual.
Tercer ajuste utilizando EFFECTPLOT por defecto pero con la opción OBS(JITTER) para mostrar las observaciones reales superpuestas con un ligero ruido aleatorio (jitter) para evitar la superposición visual.
¡Copiado!
| 1 | PROC ORTHOREG DATA=d; |
| 2 | model MBT=Time|Time|Temp|Temp; |
| 3 | effectplot / obs(jitter(seed=39393)); |
| 4 | RUN; |
| 5 | ods graphics off; |
Este material se proporciona "tal cual" por We Are Cas. No hay garantías, expresas o implícitas, en cuanto a la comerciabilidad o idoneidad para un propósito particular con respecto a los materiales o el código contenidos en este documento. We Are Cas no es responsable de los errores en este material tal como existe ahora o existirá, ni We Are Cas proporciona soporte técnico para el mismo.
Información de copyright : Frankel (1961) ... From Myers, Response Surface Methodology 1976.
