data one;
  input container temp chamber humidity $ variety $ y;
  cards;
1 50 1 low A 54
1 50 1 low B 53
1 50 1 low C 51
1 50 1 low D 51
1 50 1 low E 50
1 50 2 high A 49
1 50 2 high B 50
1 50 2 high C 52
1 50 2 high D 50
1 50 2 high E 50
2 50 3 low A 50
2 50 3 low B 49
2 50 3 low C 51
2 50 3 low D 51
2 50 3 low E 48
2 50 4 high A 46
2 50 4 high B 45
2 50 4 high C 44
2 50 4 high D 45
2 50 4 high E 45
3 50 5 low A 48
3 50 5 low B 47
3 50 5 low C 46
3 50 5 low D 46
3 50 5 low E 47
3 50 6 high A 46
3 50 6 high B 47
3 50 6 high C 42
3 50 6 high D 45
3 50 6 high E 46
4 70 7 low A 47
4 70 7 low B 47
4 70 7 low C 44
4 70 7 low D 43
4 70 7 low E 44
4 70 8 high A 41
4 70 8 high B 40
4 70 8 high C 37
4 70 8 high D 39
4 70 8 high E 40
5 70 9 low A 43
5 70 9 low B 43
5 70 9 low C 43
5 70 9 low D 41
5 70 9 low E 41
5 70 10 high A 34
5 70 10 high B 39
5 70 10 high C 35
5 70 10 high D 34
5 70 10 high E 34
6 70 11 low A 49
6 70 11 low B 48
6 70 11 low C 45
6 70 11 low D 47
6 70 11 low E 44
6 70 12 high A 30
6 70 12 high B 32
6 70 12 high C 32
6 70 12 high D 32
6 70 12 high E 31
;

proc mixed data=one method=type3;
  class container temp humidity variety;
  model y=temp
          humidity temp*humidity
		  variety temp*variety humidity*variety temp*humidity*variety;
  random container(temp) humidity*container(temp);
  lsmeans temp / cl diff;
  lsmeans humidity / cl diff;
  lsmeans variety / adjust=tukey cl diff;
run;