options ls=80 ps=50 nocenter;
data a;
   infile 'oocyte.dat' firstobs=2;
   input  for h89 date dup toocytes gv mi;
data b; set a;
   percenta=gv/toocytes;
   aspera=arsin(sqrt(percenta));
proc sort; by for h89;
proc glm;
   class h89 for;
   model aspera=h89 for h89*for/ ss1 ss2 ss3 solution;
   lsmeans h89 for / stderr;
/*
proc means; by for h89;
var percenta aspera;
*/
proc glm;
   class H89 for date;
   model aspera=date H89 for H89*for;
   lsmeans H89 for H89*for / stderr pdiff;
proc sort; by H89 for;
proc means noprint; by H89 for;
   var aspera;
   output out=plot mean=maspera;
proc plot;
   plot maspera*for=H89;