options nocenter ls = 80 ps = 64;

data drying;
   infile 'drying.dat' firstobs = 2;
   input obs cow period sample drying dmrem cprem;
proc mixed data = drying;
   class cow period sample drying;
   model dmrem = cow | period | drying;
   random sample(period) drying*sample(period);
   lsmeans period*drying;
proc mixed data = drying;
   class cow period sample drying;
   model dmrem = cprem cow | period | drying;
   random sample(period) drying*sample(period);
   lsmeans period*drying;

proc sort; by drying period;
proc plot hprecent = 33 vpercent = 33; by drying period;
    plot dmrem * cprem;