options formdlim = '-';

data trigly;
  infile 'trigly.txt' firstobs = 2; 
  input method day y;
       
proc sort;
  by method day;

proc means noprint;
  by method day;
  var y;
  output out = means mean = y;
  
proc plot;
  plot y*day = method;
  title 'Plots of means on each day';

proc mixed data = trigly method = type3;
  class method day;
  model y = method;
  random day method*day;
  
  lsmeans method;
  
  estimate 'method 1 - method 2' method 1 -1;
  title 'Broad sense inference';
  
proc mixed data = trigly method = type3;
  class method day;
  model y = day method;
  random  method*day;
  
  lsmeans method;
  
  estimate 'method 1 - method 2' method 1 -1;
  title 'Intermediate sense inference: fixed days and random method*day';  

proc mixed data = trigly method = type3;
  class method day;
  model y = day method  method*day;
  
  lsmeans method;
  
  estimate 'method 1 - method 2' method 1 -1;
  title 'Narrow sense inference: fixed days and method*day';  
run;