search for: fm2orth

Hi all, I'm trying to figure out the example about mixed models in the Pinheiro and Bates book (Mixed-Effects Models in S and S-Plus, 2000, pp. 135-137). One thing I don't understand is: When I run the command fm1Orth.lm <- lm( distance ~ age, Orthodont ) followed by fm2Orth.lm <- update( fm1Orth.lm, formula = distance ~ Sex*age ) and then do summary(fm2Orth.lm) in the output, there's "SexFemale" instead of "Sex" (see below). Why? Call: lm(formula = distance ~ Sex + age + Sex:age, data = Orthodont) Residuals: Min 1Q Median...

  Hola,   Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?   http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5   Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.   Muchas gracias.   Eva [[alternative HTML

Dear R-users: On page 155 of "Mixed-effects Models in S and S-Plus", the degree of freedoms of the anova comparison of lme and lm are 8 and 5. But when I use the following SAS code: proc glm data=ortho2; class gender; model distance = age|gender / solution ; run; The df is 3. Could you please explain this to me? Thanks Joe

...Devices utils datasets methods base other attached packages: [1] MEMSS_0.3-6 lme4_0.999375-32 Matrix_0.999375-31 lattice_0.17-26 loaded via a namespace (and not attached): [1] grid_2.10.0 tools_2.10.0 > fm1Orth.lis<-lmList(distance~age|Subject, data=Orthodont) > fm2Orth.lis<-lmList(distance~I(age-11)|Subject, data=Orthodont) Error in unique(c("AsIs", oldClass(x))) : object 'age' not found Error in unique(c("AsIs", oldClass(x))) : object 'age' not found Error in unique(c("AsIs", oldClass(x))) : object 'age' n...