/* classicalgrapefruit.sas */
title 'Grapefruit Data';

data fruit;                                   /* Skip the header */
     infile '/folders/myfolders/unigrapefruit.data.txt' firstobs=2;
     input store price sales; 

proc glm;
     title2 'Classical mixed model repeated measures';
     class price store;
     model sales = price store;
     means price;
     random store  / test;   

/* Follow up with matched t-tests, but they require a multivariate format.
   Read the data in again, with multiple lines of data per case. */

data mvfruit;
     infile '/folders/myfolders/unigrapefruit.data.txt' firstobs=2; 
     input store1 price1 sales1
           store2 price2 sales2
           store3 price3 sales3;
     One_vs_Two   = sales1-sales2;
     One_vs_Three = sales1-sales3;
     Two_vs_Three = sales2-sales3;

proc means n mean stddev t probt;
     title2 'Matched t-tests: 0.05/3 = 0.0167';
     var One_vs_Two One_vs_Three Two_vs_Three;