search for: fm1orth

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 + ag...

There is a bug in 'anova.gls' in the 'nlme' package (3.1-90). The=20 bug is triggered by calling the function with a single 'gls' object=20 and specifying the 'Terms' argument but not the 'L' argument: > library(nlme) > fm1Orth.gls <- gls(distance ~ Sex * I(age - 11), Orthodont, + correlation =3D corSymm(form =3D ~ 1 | Subject), + weights =3D varIdent(form =3D ~ 1 | age)) > anova(fm1Orth.gls) Denom. DF: 104=20 numDF F-value p-value (Intercept) 1 424...

  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

I'am trying to develop some code if R, which would correspond to what I did in SAS. The data look like: Treatment Replicate group1 GSI Control A 1 0.81301 Control B 1 1.06061 Control C 1 1.26350 Control D 1 0.93284 Low A 2 0.79359 Low B

...F-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices 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 i...