[BioC] converting probe to gene level
Naomi Altman
naomi at stat.psu.edu
Thu Jan 13 05:28:29 CET 2011
I think in general the point of algorithms like rma is that you can
do better than a simple average. For example, you might use median
polish or at least a robust method such as a trimmed mean. Of
course, if you only have 2 or 3 probes per gene, averaging is about
the best you can do.
Naomi
At 05:55 PM 1/12/2011, Wei Shi wrote:
>Hi Viritha:
>
> You do not need a function, but the code below should do this:
>
>For array1:
>tapply(x$array1, factor(x$genes), mean) # x is a data frame
>containing the data in your example.
>
>Cheers,
>Wei
>
>On Jan 13, 2011, at 8:38 AM, viritha kaza wrote:
>
> > Hi Group,
> > Is there any function which would help in converting data from probe level
> > to gene level by averaging the expression of all the probes
> corresponding to
> > a gene.
> > For eg:
> > array1 array2
> > probe 1 Gene 1 2 4
> > probe 2 Gene 1 4 3
> > probe 3 Gene 2 5 4
> >
> > result:
> > array1 array2
> > Gene1 3 3.5
> > Gene 2 5 4
> > Could any one give some ideas to perform this?
> > Thanks,
> > Viritha
> >
> > [[alternative HTML version deleted]]
> >
> > _______________________________________________
> > Bioconductor mailing list
> > Bioconductor at r-project.org
> > https://stat.ethz.ch/mailman/listinfo/bioconductor
> > Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>
>______________________________________________________________________
>The information in this email is confidential and inten...{{dropped:11}}
More information about the Bioconductor
mailing list