options ls=75;
data one;
   input plant @;
   do plot=1 to 3;
      input protein @;
      output;
   end;
cards;
1 42.4  41.0  39.6
2 28.6  36.3  42.2
3 43.2  42.1  40.2
4 40.8  41.0  38.9
5 41.0  38.3  41.1
6 39.4  39.5  37.2
7 39.6  40.4  38.9
8 38.1  38.3  37.9
9 35.9  36.1  35.6
10 39.6  39.9  39.7
;
run;
proc glm;
   class plant;
   model protein=plant;
   random plant;
   means plant;
   output out=two p=ymean r=r;
run;

proc gplot;
   plot protein*plant;
   plot r*plant;
   plot r*ymean=plant;
run;

proc univariate plot;
   var r protein;
   histogram r protein;
run;
/* What's going on?  Any of these data badly behaved? */

/*proc mixed;
   class plant;
   model protein=/s cl;
   random plant;
run;*/