options nocenter ps=50 ls=76;
data hardy;
   infile 'hardyraw.dat' firstobs=2;
   input potato regime temp photo leak;
   /* 0...7 plotting code */
   code = 4 * (potato - 1) + 2 * regime + (temp - 1);
   leak = sqrt(leak+1.11); /* stabilize the variance */

proc print;
   var potato regime temp code photo leak;
run;