Hi Katie, maybe the easiest solution is to create a new factor that corresponds to the combinations of the three factors A, B and C. A quick and dirty way to create such a factor is: ABC <- factor(paste(A, "x", B, "x", C, sep = "")) ABC and then fit the model using the variable ABC instead of A*B*C. Christian