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 PRINT;

proc mixed data=influent covtest cl;
   class influent;
   model y = /solution;
   random influent / solution;

   estimate 'influent 1' intercept 1 | influent 1 0 0 0 0 0;
   estimate 'influent 2' intercept 1 | influent 0 1 0 0 0 0;
   estimate 'influent 3' intercept 1 | influent 0 0 1 0 0 0;
   estimate 'influent 4' intercept 1 | influent 0 0 0 1 0 0;
   estimate 'influent 5' intercept 1 | influent 0 0 0 0 1 0;
   estimate 'influent 6' intercept 1 | influent 0 0 0 0 0 1;

   estimate 'influent 1U' | influent 1 0 0 0 0 0;
   estimate 'influent 2U' | influent 0 1 0 0 0 0;
   estimate 'influent 3U' | influent 0 0 1 0 0 0;
   estimate 'influent 4U' | influent 0 0 0 1 0 0;
   estimate 'influent 5U' | influent 0 0 0 0 1 0;
   estimate 'influent 6U' | influent 0 0 0 0 0 1;
run;