[R] structure of mlm objects ?

Prof Brian Ripley ripley at stats.ox.ac.uk
Wed Feb 25 11:34:27 CET 2004

On Wed, 25 Feb 2004, ZABALZA-MEZGHANI Isabelle wrote:

> I am using the function "lm" to fit several responses at the same time (100
> responses). At the end of the fit, I get an object of class "mlm".

Actually of class c("mlm", "lm").

> I would like to know if there is a way to access to each of the 100
> underlying models separately (is it a list, ... ?). Which syntax should I
> use to see and use the 15th model (for instance) just like it is possible
> for classical "lm" objects.

You don't have 100 underlying fits.  All the information is for one set 
of x's and 100 sets of y, and that relating only to x is stored only once.
The following simple function will extract the i'th fit.

mlm2lm <- function(fit, i)
    for(k in c("coefficients", "residuals", "effects", "fitted.values"))
        fit[[k]] <- fit[[k]][, i]
    class(fit) <- "lm"

albeit with the wrong call, so use carefully.

Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

More information about the R-help mailing list