Hi,

I need help for using one way anova for microarray data.
In 20 files, I have 5 treatment and 4 replicates for each treatment.

I have the following code:

.
.
.
treatment <- gl(5,4,20, label=c("T1","T2","T3","T4","T5"))
# define ANOVA function
aof <- function(x) {
  m<-data.frame(treatment, x);
  anova(aov(x ~ treatment, m))
}

anovaresults <- apply(sdata, 1, aof)


Can you suggest how can I do one way ANOVA for 5 treatment and to make
all the comparations between treatments?

	[[alternative HTML version deleted]]

