[R] vectorizing ANOVA over a vectorized linear model

Mark Kimpel
Sun, 07 Mar 2010 18:03:03 -0000

Is it possible to vectorize anova over the output of a vectorized lm?  I
have a gene expression matrix with each row being a gene and columns for
samples. There are several factors with interactions. I can get p values by
looping over the matrix with lm and anova, but I would like to make this as
computationally efficient as possible. I am able to vectorize the lm
command, but when I try to use anova on the resultant model object I get
just one anova result.

Is what I want to do possible? And, yes, I am quite conversant with Limma
and other BioC packages, I have my reasons for wanting to use lm and anova.

Thanks,

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) 399-1219 Skype No Voicemail please

	[[alternative HTML version deleted]]