XCHART Statement: ANOM Procedure

See ANMXEX6 in the SAS/QC Sample LibraryThis example illustrates the use of analysis of means in an experiment with two factors where an interaction effect is present. The following data set cleaninginteract is a modified version of the data set Cleaning, which includes an interaction effect for position and depth.

Consider the following data set cleaninginteract:

data cleaninginteract;  
   do position = 1 to 5;
      do depth = 1 to 3;
         do rep = 1 to 2;
            input concentration @@;
            output;
         end; 
      end;
   end;
   datalines;
   15 16 15 14 19 5
   15 16 14 14 0  1
   19 15 16 16 11 8
   18 16 24 23 8  14
   15 12 23 24 8  11
   ;  
;

The following statements use PROC GLM to test for an interaction:

ods graphics off;
proc glm data=cleaninginteract;
   class position depth;
   model concentration = position depth position*depth;
run;

The analysis of variance results in Output 4.7.1 indicate a significant interaction between position and depth.

Since an interaction effect is present, an appropriate way to analyze the data is to create an ANOM chart for the cell means.

In order to create the chart you first need to compute the cell means and a new group variable which designates the cells. The following statements use PROC MEANS for this purpose.

proc means data=cleaninginteract n mean std;
   class position depth; 
   var concentration;
   types position*depth;
   output out=cellmeans mean=concentrationX std=concentrationS;
run;
data cellmeans; set cellmeans;
   rename _FREQ_ = concentrationN;
   pos = put(position, z1.);
   dep = put(depth, z1.);
   cell  = cat('P',pos, 'D', dep);
   drop _TYPE_ pos dep;
run;

The cell means are stored in the data set cellmeans shown in Output 4.7.2:

The data set cellmeans has the structure of a SUMMARY= input data set for the ANOM procedure. For details concerning a SUMMARY= data set, see the section Creating ANOM Charts for Means from Group Summary Data.

The following statements use cellmeans to create the ANOM chart for the cell means using SUMMARY= option:

title "ANOM for Cell Means of Position and Depth";
proc ANOM summary = cellmeans;
   xchart concentration * cell / turnhlabels;      
   label concentrationX = 'Mean of Concentration';
   label cell           = 'Cell';
run;

The chart is shown in Output 4.7.3:

Output 4.7.3: ANOM for Cell Means of Position and Depth

The chart shows that the cell means for P2D3, P4D2, and P5D2 are significantly different from the average concentration level.