options ls=75 nodate;

data one;
   infile "peanuts3.dat" firstobs=5;
   input junk location supplement $ ppfd y CP;
run;

/*
proc glm ;
   class ppfd location supplement ;
   model CP=ppfd|supplement location|ppfd;
   random location location*ppfd;
   means ppfd*supplement;
   contrast "linear" ppfd -3 -1 1 3;
   contrast "nonlinear" ppfd 1 -1 -1 1;
run;
*/
options ls=95 nodate;
proc mixed method=type3;
   class ppfd location supplement ;
   model CP=ppfd|supplement/ddfm=satterth clparm; 
   random location location*ppfd;
   *lsmeans ppfd|supplement/diff ;
   estimate "theta1" ppfd -1 0  0 1/cl; 
   estimate "theta2" ppfd -1 0  0 1 ppfd*supplement 0 -1 0 0 0 0 0 1/cl;
   estimate "theta3" supplement -1 1 ppfd*supplement 0 0 0 0 0 0 -1 1/cl;
   estimate "linear" ppfd -3 -1 1 3;
   estimate "nonlinear" ppfd 1 -1 -1 1;
   contrast "linear" ppfd -3 -1 1 3;
   contrast "nonlinear" ppfd 1 -1 -1 1;
run;