Thanks both Harold Doran and Prof. Ripley for the suggestion.
Time*Group - 1 or Time*(Group-1) does seem better. However as Prof.
Ripley pointed out, it is a little complicated with the interactions.
For example,
==
> set.seed(1)
> group <- as.factor (sample (c("M","F"), 12, T))
> y <- rnorm(
On Wed, 16 Apr 2008, Doran, Harold wrote:
> R may not be giving you what you want, but it is doing the right thing.
> You can change what the base category is through contrasts but you can't
> get the marginal effects for every level of all factors because this
> creates a linear dependence in the
R may not be giving you what you want, but it is doing the right thing.
You can change what the base category is through contrasts but you can't
get the marginal effects for every level of all factors because this
creates a linear dependence in the model matrix.
> -Original Message-
> Fro
3 matches
Mail list logo