Hello, Is it possible to cast the output of lm.fit to an lm object? I've 10,000 linear models for a gene expression experiment, all of which have the same model matrix. Maybe calling lm.fit on a model matrix and a data vector is faster than lm. I'd like to use each fit for an anova as well as comparing different models via anova. kind regards, Arne