/* ------------------------------------------------
	File: plotMeans.genotype.sas

	Purpose: Plot the means for the Genotype data
		set (see genotypes.sas)
    ----------------------------------------------- */

PROC SORT DATA=p7291.genotype OUT=Temp1;
	BY Genotype;
RUN;

PROC FORMAT;
	VALUE dvfmt
		0 = ' '
		1 = 'Y1'
		2 = 'Y2'
		3 = ' ';
RUN;


DATA Temp1;
	SET p7291.genotype;
	Depvar=1; Score = Y1; OUTPUT;
	Depvar=2; Score = Y2; OUTPUT;
	format depvar dvfmt.;
RUN;

PROC SORT DATA=Temp1;
	BY Depvar Genotype;
RUN;

PROC MEANS DATA=Temp1;
	BY Depvar Genotype;
	VAR Score;
	OUTPUT OUT=TempMeans MEAN=Mean STDERR=StdErr;
RUN;

PROC PRINT;
RUN;

DATA TempPlot;
	SET TempMeans;
	Score = Mean; OUTPUT;
	Score = Mean + StdErr; OUTPUT;
	Score = Mean - StdErr; OUTPUT;
	LABEL
		Score = 'Mean'
		Depvar = 'Dependent Variable';
RUN;


TITLE J=l Means for Data Set Genotype;
AXIS1;
AXIS2 order=0 to 3 by 1;
SYMBOL1 c=black f=, v=none l=1 w=4 i=hilotj ci=red;
SYMBOL2 c=black f=, v=none l=2 w=4 i=hilotj ci=green;
SYMBOL3 c=black f=, v=none l=3 w=4 i=hilotj ci=blue;
PROC GPLOT DATA=TempPlot;
	PLOT Score*Depvar=Genotype / VAXIS=axis1 HAXIS=axis2;
RUN;
QUIT;