search for: myresponsematrix

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...