search for: martinl

...rtin, Have you try to create a new variable for age squared, say agesq? If you fit the model using this new variable you should get the coefficients. So your new model is something like height~age+agesq I hope this helps, Saghir > -----Original Message----- > From: MARTIN Ludovic [SMTP:martinl@mathinfo.ens.univ-reims.fr] > Sent: Wednesday, 07 May, 2003 2:02 PM > To: r-help@stat.math.ethz.ch > Subject: [R] [R ] Query : problems with the arithmetic operator "^" > with function "lme" and "lmList" > > Dear all, > > I've got a proble...

Dear all, I've got a problem in including square variables in lme function. I've tried to work on Dialyzer data of Pinheiro and Bates'book. We fit the heteroscedastic model with: > data(Dialyzer) > fm2Dial.lme<-lme(rate~(pressure+pressure^2+pressure^3+pressure^4)*QB, + Dialyzer,~pressure+pressure^2,weights=varPower(form=~pressure)) We Obtain > fm2Dial.lme Linear

Dear all, I've got a problem in including square variables in lme and lmlist functions. I've tried to work on Oxboys data of Pinheiro and Bates'book, which consist of the heights of 26 boys, each mesured on nine different occasions : > Oxboys Grouped Data: height ~ age | Subject Subject age height Occasion 1 1 -1.0000 140.50 1 2 1 -0.7479 143.40