options ps=54 ls=76;

data swamp;
infile 'C:\data\swamp1.dat';
input size swampid watlev;
run;

proc means mean;
class size;
var watlev;
run;

proc means mean std;
class swampid;
var watlev;
run;

proc glm;
class size swampid;
model watlev = size swampid(size);
test h=size e=swampid(size);
run;

quit;