I have fit a model to ~20k different genes and would now like to extract the
p-val for one of the effects, again for each individual gene.

My model code is:
mod <- lm(myResponseMatrix~ Time)

Were I to do this with a response vector rather than matrix, the following
would work:
anovaResult<- anova(mod)

With a more complex lm object, such as that generated above, it seems anova
is trying to compare all the models to each other, a Herculean task for ~20k
models and NOT what I want.

Without resorting to a slow loop/apply approach, is there a way to get
lightening fast anova results in a fashion similar to that obtained from lm?

Mark

-- 
Mark W. Kimpel MD ** Neuroinformatics ** Dept. of Psychiatry
Indiana University School of Medicine

15032 Hunter Court, Westfield, IN 46074

(317) 490-5129 Work, & Mobile & VoiceMail
(317) 663-0513 Home (no voice mail please)

******************************************************************

	[[alternative HTML version deleted]]

