Displaying 1 result from an estimated 1 matches for "myresponsematrix".
2008 Apr 03
1
getting serial anovas from a complex lm object obtained with a matrix of responses
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...