Los datos ('length' y 'mass') utilizados para el análisis se crean directamente dentro del script a través de una instrucción DATA STEP con 'datalines', lo que significa que están integrados en el código fuente y no provienen de una fuente externa.
1 Bloque de código
DATA STEP Data
Explicación : Este bloque DATA STEP crea un nuevo conjunto de datos temporal llamado 'lizard'. Define dos variables numéricas, 'length' y 'mass', y las llena con los datos brutos proporcionados directamente en el script a través de la instrucción DATALINES. Este conjunto de datos es la fuente de los datos para los análisis estadísticos y las visualizaciones posteriores.
Explicación : Este procedimiento PROC NLIN ajusta un modelo de regresión no lineal a los datos del conjunto de datos 'lizard'. El modelo especificado es un 'modelo de meseta izquierda', donde la variable dependiente 'mass' es constante ('b0') hasta que 'length' alcanza 'tau', luego sigue una relación lineal. Los parámetros 'b0', 'b1' y 'tau' se inicializan para el procedimiento. La opción 'output out=a p=p' crea un nuevo conjunto de datos 'a' que incluye las variables originales de 'lizard', así como la variable 'p', que representa los valores de masa predichos por el modelo.
¡Copiado!
proc nlin;
parameters b0 = 0.2904 b1 = 0.0189 tau = 23.44;
if length <= tau then do;
model mass = b0;
end;
else do;
model mass = b0 + b1*(length-tau);
end;
output out=a p=p;
run;
1
PROC NLIN;
2
parameters b0 = 0.2904 b1 = 0.0189 tau = 23.44;
3
IFLENGTH <= tau THENDO;
4
model mass = b0;
5
END;
6
ELSEDO;
7
model mass = b0 + b1*(LENGTH-tau);
8
END;
9
OUTPUT out=a p=p;
10
RUN;
3 Bloque de código
PROC GPLOT
Explicación : Este bloque inicializa las opciones gráficas globales (`GOPTIONS`) para definir la apariencia de los gráficos (colores, fuentes, tamaños de texto). También configura un título principal (`title1`), definiciones de símbolos (`symbol1`, `symbol2`) y especificaciones de ejes (`axis1`, `axis2`) que incluyen las etiquetas y los rangos. Luego se utiliza el procedimiento `PROC GPLOT` para generar un diagrama de dispersión simple ('mass' vs 'length') a partir del conjunto de datos 'a', utilizando las definiciones de ejes previamente establecidas.
¡Copiado!
goptions reset=global gunit=pct border cback=white
colors=(black blue green red)
ftitle=swissb ftext=swiss htitle=4 htext=4;
title1 'Mass vs Length with the fit';
symbol1 color=red
interpol=none
value=dot
height=3;
symbol2 color=red
interpol=join;
axis1 label=('Length (mm)')
order=(22 to 28 by 1)
width=3;
axis2 label=('Mass (g)')
order=(0.25 to .45 by 0.05)
width=3;
proc gplot data=a;
plot mass*length/ haxis=axis1 vaxis=axis2;
run;
Explicación : Este último bloque utiliza `PROC GPLOT` para crear un gráfico superpuesto. Traza los datos observados ('mass' vs 'length') y los valores predichos por el modelo ('p' vs 'length') en el mismo gráfico a partir del conjunto de datos 'a'. La opción `overlay` permite visualizar ambas series en el mismo sistema de coordenadas, facilitando la comparación entre las observaciones reales y el ajuste del modelo. Se reutilizan las configuraciones de ejes definidas previamente.
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.
Documentación relacionada
Aucune documentation spécifique pour cette catégorie.
SAS y todos los demás nombres de productos o servicios de SAS Institute Inc. son marcas registradas o marcas comerciales de SAS Institute Inc. en los EE. UU. y otros países. ® indica registro en los EE. UU. WeAreCAS es un sitio comunitario independiente y no está afiliado a SAS Institute Inc.
Este sitio utiliza cookies técnicas y analíticas para mejorar su experiencia.
Saber más.