/* Data on Wonka(R) Bottle Caps candies.
   Halloween 2011.
   Rick Wicklin, SAS, http://blogs.sas.com/content/iml */
data Counts;
input Cola RootBeer Grape Cherry Orange;
Package = _N_;
datalines;
0 0 3 0 0
0 1 1 0 1
1 1 0 0 1
0 1 2 0 0
1 0 0 1 1
1 0 2 0 0
1 1 0 1 0
2 1 0 0 0
1 0 1 0 1
1 0 1 1 0
0 0 0 2 1
0 1 0 2 0
0 0 2 0 1
0 1 2 0 0
1 1 0 1 0
1 0 0 1 0
0 1 1 0 1
0 1 0 0 2
0 1 0 2 0
0 0 0 1 2
0 1 0 2 0
0 1 1 0 1
1 0 1 0 1
0 1 1 0 1
0 0 0 2 1
2 0 0 0 1
1 1 1 0 0
0 1 0 2 0
0 0 0 3 0
0 0 1 0 2
1 0 1 1 0
0 1 0 0 2
0 0 0 1 2
0 0 0 0 3
;
run;

/* Transpose from wide (1 subject var, 5 count vars)
   to long format (Flavor, Count)
   See http://blogs.sas.com/content/iml/2011/01/31/reshaping-data-from-wide-to-long-format/ */
proc transpose data=Counts 
   out=BottleCaps(rename=(Col1=Count)) name=Flavor;                    
   by Package;                   /* for each subject */
   var Cola RootBeer Grape Cherry Orange; /* make a row for these variables */
run;

proc freq data=BottleCaps;
label Flavor = "Flavor";
weight Count;
tables Flavor / nocum chisq plots=deviationplot;
run;