Hi,

I conducted a three-factorial [ pH (2 levels), salt(3 levels), and iron (3 levels) in a completely randomized design conducted with 3 replications. I want to do an orthogonal test on the second and third levels of iron at the first levels of pH and salt in SAS. However, I get the message NOTE: CONTRAST one is not estimable. Could you please tell me where I am making mistake?

Data factorial;

input ph salt iron rep yield;

cards;

1 1 1 1 3.5

1 1 1 2 3.8

1 1 1 3 3.4

1 1 2 1 4.1

1 1 2 2 4.3

1 1 2 3 4.5

1 1 3 1 4.6

1 1 3 2 4.8

1 1 3 3 4.9

1 2 1 1 3.1

1 2 1 2 3.2

1 2 1 3 3

1 2 2 1 3.4

1 2 2 2 3.5

1 2 2 3 3.3

1 2 3 1 3.5

1 2 3 2 3.8

1 2 3 3 4

1 3 1 1 3

1 3 1 2 2.98

1 3 1 3 2.8

1 3 2 1 3.1

1 3 2 2 3.5

1 3 2 3 3.2

1 3 3 1 3.6

1 3 3 2 4.1

1 3 3 3 3.4

2 1 1 1 3.2

2 1 1 2 3.4

2 1 1 3 3.1

2 1 2 1 4.1

2 1 2 2 4

2 1 2 3 4.2

2 1 3 1 4.3

2 1 3 2 4.5

2 1 3 3 4.3

2 2 1 1 2.98

2 2 1 2 2.8

2 2 1 3 2.7

2 2 2 1 3.2

2 2 2 2 3.1

2 2 2 3 3.3

2 2 3 1 3

2 2 3 2 3.4

2 2 3 3 3.6

2 3 1 1 2.65

2 3 1 2 2.75

2 3 1 3 2.5

2 3 2 1 3

2 3 2 2 2.98

2 3 2 3 2.95

2 3 3 1 3.2

2 3 3 2 3.5

2 3 3 3 3.1

;

proc glm;

class ph salt iron rep;

model yield=ph salt iron ph*salt ph*iron salt*iron ph*salt*iron/ss3;

contrast 'one' iron 0 1 -1

ph*salt*iron 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0/e;

run;