option ls=76 nocenter;
data influent;
  input influent y type @@;
  datalines;
1   21  2 1   27  2 1   29  2 1   17  2 1   19  2
1   12  2 1   29  2 1   20  2 1   20  2 2   21  2
2   11  2 2   18  2 2   9   2 2   13  2 2   23  2 
2   2   2 3   20  1 3   19  1 3   20  1 3   11  1
3   14  1 4   14  2 4   24  2 4   30  2 4   21  2
4   31  2 4   27  2 5   7   1 5   15  1 5   18  1
5   4   1 5   28  1 6   41  3 6   42  3 6   35  3
6   34  3 6   30   3
;


proc mixed data=influent covtest cl; 
   class type influent;
   model y=type/solution;
   random influent(type)/solution;
   estimate 'influent 1' intercept 1 type 0 1 0 | 
                         influent(type) 1 0 0 0 0 0;
   estimate 'influent 2' intercept 1 type 0 1 0 | 
                         influent(type) 0 1 0 0 0 0;
   estimate 'influent 3' intercept 1 type 1 0 0 | 
                         influent(type) 0 0 1 0 0 0;
   estimate 'influent 4' intercept 1 type 0 1 0 | 
                         influent(type) 0 0 0 1 0 0;
   estimate 'influent 5' intercept 1 type 1 0 0 | 
                         influent(type) 0 0 0 0 1 0;
   estimate 'influent 6' intercept 1 type 0 0 1 | 
                         influent(type) 0 0 0 0 0 1;
   lsmeans type / diff;
run;