Visualización Estadística Anotada con PROC GPLOT

Este script primero genera un conjunto de datos de prueba que contiene diversas medidas físicas clasificadas por tipo (TEST). Luego utiliza el procedimiento MEANS para calcular la media, el mínimo y el máximo para cada categoría. Un paso DATA complejo construye un conjunto de datos de anotación (Annotate Facility) para dibujar manualmente barras de error horizontales que representan el intervalo Mín-Máx y marcar la media. Finalmente, el procedimiento GPLOT (SAS^©/GRAPH) muestra el resultado. Nota: Este código utiliza procedimientos gráficos antiguos (SAS^©/GRAPH) en lugar de los modernos procedimientos de ODS Graphics (SGPLOT).

Análisis de datos

Type : CREATION_INTERNE

Los datos se definen directamente en el script a través de la instrucción CARDS en el paso DATA 'a'.

1 Bloque de código

DATA STEP Data

Explicación :
Inicialización de las opciones gráficas globales y creación del conjunto de datos fuente 'a' que contiene las variables TEST (categoría) y BREAKS (valor).

¡Copiado!

1	/* Set graphics options */
2	goptions reset=all cback=white border;
3
4	/* Set system options */
5	DATA a;
6	INPUT TEST $ BREAKS ;
7	CARDS;
8	Cold 5
9	Cold 12
10	Cold 14
11	Cold 22
12	Cold 52
13	Heat 20
14	Heat 25
15	Heat 10
16	Heat 22
17	Heat 47
18	Gases 12
19	Gases 25
20	Gases 33
21	Gases 48
22	Gases 24
23	Pressure 10
24	Pressure 12
25	Pressure 14
26	Pressure 22
27	Pressure 60
28	Xrays 20
29	Xrays 25
30	Xrays 14
31	Xrays 22
32	Xrays 29
33	Humidity 20
34	Humidity 25
35	Humidity 33
36	Humidity 40
37	Humidity 24
38	;
39	RUN;

2 Bloque de código

PROC SORT

Explicación :
Ordenamiento de los datos por la variable de grupo 'TEST' para preparar la agregación.

¡Copiado!

1	/* Sort
2	data by variable TEST */
3	PROC SORT;
4	BY TEST;
5
6	RUN;
7

3 Bloque de código

PROC MEANS Data

Explicación :
Cálculo de las estadísticas descriptivas (media, desviación estándar, min, max) agrupadas por 'TEST'. Los resultados se almacenan en la tabla de salida 'b'.

¡Copiado!

1	/**************************************************/
2	/* Create an output data set, B using PROC MEANS */
3	/* that contain new variables, MEAN, STD, STDERR, */
4	/* MIN, and MAX. */
5	/**************************************************/
6
7	PROC MEANS mean std stderr min max DATA=a;
8	BY TEST;
9	OUTPUT out=b mean=mean min=min max=max;
10	RUN;

4 Bloque de código

DATA STEP Data

Explicación :
Creación de un conjunto de datos especial 'Annotate'. Define comandos gráficos (move, draw) para dibujar líneas azules del min al max, y marcadores rojos para la media, el min y el max. Utiliza el sistema de coordenadas de datos (xsys='2', ysys='2').

¡Copiado!

1	/****************************************************************/
2	/* Create an annotate data set, ANNO to draw the bars at +/- 1, */
3	/* 2, or 3 Standard Deviation or Standard Error of the mean. */
4	/****************************************************************/
5
6	DATA anno;
7	retain xsys ysys '2' when 'a';
8	LENGTH color function $8 ;
9	SET b;
10
11	/* Draw the horizontal line from min to max */
12	function='move'; xsys='2'; ysys='2'; yc=TEST; x=min; color='blue'; OUTPUT;
13	function='draw'; x=max; color='blue'; size=2; OUTPUT;
14
15	/* Draw the MEAN horizontal line making the SIZE bigger */
16	function='move'; xsys='2';ysys='2';yc=TEST;x=mean; color='red'; OUTPUT;
17	function='draw'; x=mean; ysys='9'; y=+2; size=4; OUTPUT;
18	function='draw'; x=mean; y=-4; size=4; OUTPUT;
19
20	/* Draw the line for the MIN value */
21	function='move';xsys='2';ysys='2';yc=TEST;x=min;color='red';OUTPUT;
22	function='draw';x=min;ysys='9';y=+2;size=2;OUTPUT;
23	function='draw';x=min;y=-4;size=2;OUTPUT;
24
25	/* Draw the line for the MAX value */
26	function='move';xsys='2';ysys='2';yc=TEST;x=max;color='red';OUTPUT;
27	function='draw';x=max;ysys='9';y=+2;size=2;OUTPUT;
28	function='draw';x=max;y=-4;size=2;OUTPUT;
29	RUN;

5 Bloque de código

PROC GPLOT

Explicación :
Generación del gráfico final. El procedimiento GPLOT utiliza el conjunto de datos 'anno' para superponer los dibujos personalizados a los ejes definidos. El eje horizontal se fija de 0 a 100.

¡Copiado!

1	axis1 order=(0 to 100 BY 10);
2	symbol1 i=none v=none c=black;
3
4	PROC GPLOT DATA=b ;
5	plot testmean / anno=anno haxis=axis1 href=30 60 90; / The HREF= option draws reference lines */
6	RUN;
7	QUIT;

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.

Retour à la liste