Dear R-help,
I posted this on 4 Oct but got no response (I wasn't even told to go away and
do some more background reading ;) ). I am reposting it in the, perhaps, vain
hope that someone with knowledge of the subject will reply, if only to point me
in a different direction to which I am now
I f I understand you correctly, I don't think that your model is doing
what you think it is.
Look at the model.matrix. Consider a toy example:
x - 1:10
y - factor(letters[1:2])
dd - expand.grid(x, y)
dd$resp - rnorm(20)
model.matrix(~Var2/Var1+1, dd)
(Intercept) Var2b Var2a:Var1