ok I created a matrix C with
A B CA1 B1C1
..
.
the columns contain the gene expression values..
I ran the following t.test:
apply(C, 1, function(x) t.test( x[1:3], x[4,6] )$p.value )
which outputs out 16063 pvalues
Mark,
Please do not post the same question to both R-help and BioC-help
mailing lists because 1) there are many people who are on both lists and
2) people's replies will be archived in two different places making it
harder to others to search in future.
Please see the responses on BioC-help