options nocenter ps=64 ls=80;
data a;
   infile 'running.dat' firstobs = 2;
   input sub ord sp sf md $ vo2 tm etot cn ar ecn ahm vhm;
   lvo2 = log10(vo2);
proc glm;
   class sp sf md;
   model lvo2 = sp | sf | md / ss3;
   lsmeans sp*sf*md / stderr;
proc glm;
   class sub sp sf md;
   model lvo2 = sub sp sf md / ss3 solution;
   lsmeans sp sf md / pdiff stderr;
   output out = b predicted = plvo2 residual = rlvo2;
data c; set b;
   code = 0;
   if sp > 3   then code = 4;
   if sf > 140 then code = code + 2;
   if md = 'L' then code = code + 1;
/* use first of each pair to connect points by subject */
proc plot hpercent = 50 vpercent = 50;
   plot ord * lvo2 = sub;
   plot ord * lvo2 = code;
   plot ord * rlvo2 = sub;
   plot ord * rlvo2 = code;