On 9/29/07, Dieter Menne <[EMAIL PROTECTED]> wrote: > On 9/28/07, Anouk Simard <Anouk.Simard <at> bio.ulaval.ca> wrote: > > > > I would like to extract the fitted values from a model using LMER but > > > only for the fix portion of the model and not for the fix and random > > > portion (e.g it is the procedure outpm or outp in SAS). > .. > Quoting Douglas Bates <bates <at> stat.wisc.edu> > > > > The calculation of the level = 0 fitted values in the new > > > representation of the fitted model is quite easy. It is > > > > > > [EMAIL PROTECTED] %*% fixef(fm1)" > > > > > Douglas Bates replied > > > Yes. > > Mmm.. If I consider a few hundreds of messages in this list on the subject, > it's considered dirty to use the [EMAIL PROTECTED] or fm1$X construct. So > should we > expect that this is final (no longer fitted(), augPred etc), or is it work > in prooogresssssssss?
This is to indicate that I am rather slow in getting the lme4 package finished? :-) I agree that I am and no one is more impatient with the progress than I. You are correct that it is considered bad form to reach into a fitted model structure and pull out individual components or slots. My brief answer agreeing that this was still the best way of obtaining these particular values should have emphasized the "still". I haven't thought of a good way of specifying this particular form in, say, the "fitted" method. There is a clean way of specifying the levels of random effects incorporated into the fitted values for models with nested random effects only. However, for crossed or partially crossed random effects this is less clearly defined. Perhaps a better alternative would be to define a method for the "model.matrix" generic for lmer models. Then the calculation could be written as model.matrix(fm1, type = "fixed") %*% fixef(fm1) using only extractor functions. It may be surprising that sometimes the most difficult part of the design of the code in a package is deciding on the names and arguments of functions or methods. In my experience rushing such decisions often leads to a lot of maintenance headaches. ______________________________________________ R-help@r-project.org 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.
