Re: [R] How to prevent inclusion of intercept in lme with interaction

2009-04-01 Thread Dieter Menne
ONKELINX, Thierry inbo.be> writes: > With t*treat the model allows for a different slope AND a different > intercept for each treatment. If you only want different slopes and all > intercepts equal to 0, then t:treat - 1 or t + t:treat - 1 is the model > you are looking for. Thanks Thierry and C

Re: [R] How to prevent inclusion of intercept in lme with interaction

2009-04-01 Thread Christian Ritz
Hi Dieter, the following model assumes a linear relationship between the response "newbone" and the independent variable "t" with a common intercept equal to 0 and treatment-dependent slopes: grd.lme0 <- lme(newbone~t:treat-1, data=grd, random=~1|subject) summary(grd.lme0) Christian

Re: [R] How to prevent inclusion of intercept in lme with interaction

2009-04-01 Thread ONKELINX, Thierry
Verzonden: woensdag 1 april 2009 14:11 Aan: r-h...@stat.math.ethz.ch Onderwerp: [R] How to prevent inclusion of intercept in lme with interaction Dear friends of lme, After so many year with lme, I feel ashamed that I cannot get this to work. Maybe it's a syntax problem, but possibly a la

[R] How to prevent inclusion of intercept in lme with interaction

2009-04-01 Thread Dieter Menne
Dear friends of lme, After so many year with lme, I feel ashamed that I cannot get this to work. Maybe it's a syntax problem, but possibly a lack of understanding. We have growth curves of new dental bone that can well be modeled by a linear growth curve, for two different treatments and several