"David Barron" <[EMAIL PROTECTED]> writes:

> I'm not sure if this will help, but it's worth a try.  Do the
> regression as I suggested before, extract the model matrix and remove
> the "offending" column.  I'm assuming you don't know in advance how
> many levels there are in the factor.  Then use this to perform the
> regression.  Something like this:
>
> m1 <- lm(x ~ fd:y + fd)
> mm <- model.matrix(m1)
> nl <- length(levels(fd))
> newdat <- mm[,-c(1,nl)]
> lm(x ~ newdat)

Thanks a lot! With a small change that solved my problem!

-- 
Andreas Gösele                   Omnis enim res, quae dando non deficit,
Inst. f. Gesellschaftspolitik    dum habetur et non datur,
Kaulbachstr. 31a, 80539 München  nondum habetur, quomodo habenda est.
E-mail: [EMAIL PROTECTED]      (Augustinus)

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to