Hi I can do this: formula = as.factor(outcome) ~ . in glm and other model building functions. I think there is a way to get the product of the determinants (that is d1 * d2, d1 * d3, etc) and also another way to get all the polynomials (that is like poly(d1,2) would produce for a single determinant). Can anyone tell me how you write them? Stephen [[alternative HTML version deleted]]