If you want a matrix (perhaps a data.frame so you can add different data
types(gene name, accession, etc) of just the F and Ps from each anova.
Something like this works for me. Though I'm a newbie too, so if there
is a better way someone let me know!
----------------------
#create 2 by # genes matrix, can add row names or col names
#precomputing your matrix will save some computational time as compare
to rbind
#assume genes is some matrix or data.frame
results <- matrix(ncol=2, nrow=length(genes[,1]))
for(i in 1:length(genes[,1]){
an <- anova(lm(formula,data))
results[i,] <- c(an$F[1], an$P[1] )
}
------------------
Don't forget about the multiple testing!
-Mike
Joshi, Nina (NIH/NCI) wrote:
>
>
>I am running the aov function to determine the linear relationship between
>individual genes and a number of covariates in a micro-array data set. My
>aov function is working, but I have not been able to write a code to save
>these results in a matrix. I am using the following code or a slight
>variation of this code:
>
>
>
>for(i in 1(length(m[1,])-20)
>
>
>
>{tmp<- aov(m[,i]~treat, data=m)
>
>
>
>{results [i] <- tmp}}
>
>
>
>Any help would be appreciated.
>
>
>
>
>
>Thanks,
>
>
>
>Nina
>
>
>
>
>
>Nina Joshi, PhD
>
>NIH/NCI/ Genetics Branch
>
>National Naval Medical Center, Bldg. 8, Rm. 5101
>
>8901 Wisconsin Ave.
>
>Bethesda, MD. 20889-5101
>
>(301) 435-5436 - phone
>
>(301) 496-0047 - fax
>
>
>
>
>
>
> [[alternative HTML version deleted]]
>
>______________________________________________
>R-help at stat.math.ethz.ch mailing list
>https://www.stat.math.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
>
>
--
Michael Morley
Bioinformatics Specialist
University of Pennsylvania
Department of Pediatrics
3516 Civic Center Blvd.,
510B Abramson Pediatric Research Center,
Philadelphia, PA 19104-4318.
Phone: (215) 590-7673
FAX: (215) 590-3709