options nocenter ls=79 ps=100;

data soil;
   infile 'maxwell.dat';
   input soil$ rep acc$ score1 score2 score3 ht1 score4 ht2 score5 ht3 score6;
   if acc = 'NOTAUTO';

data times; set soil;
   score = score1; ht = .; time = 1; output;
   score = score2; ht = .; time = 2; output;
   score = score3; ht = ht1; time = 3; output;
   score = score4; ht = ht2; time = 4; output;
   score = score5; ht = ht3; time = 5; output;
   score = score6; ht = .; time = 6; output;
   drop score1--score6;

proc sort; by soil rep;
proc means noprint; by soil rep;
   var score ht;
   output out=mean mean=mscore mht;
proc print;
   var soil rep mscore mht;
run;