Question: Workflows For Microarray Meta-Analysis Projects.
gravatar for Ankur
6.2 years ago by
Ankur90 wrote:

Hi everyone.

I have been contemplating doing a few analyses of publicly available affy array data using R and Bioconductor to define signatures for phenotypes/identify transcriptional biomarkers et cetera. I always work using raw data and the options I have in terms of workflows are

[1] Pool together all the .CEL files, then run it through RMA and limma in one go. [2] Normalise arrays from individual studies with their respective batches, then combine normalised expression values into one expression set for further analysis. [3] Try combining P-values using Stouffer's z, for instance.

Previously, my approach involved looking across differentially expressed genes for each study addressing a question to see which genes were recurrent, but given issues associated with dodgy datasets/small datasets with high adjusted P.values introducing lots of false negatives I am not a fan.

Which workflow would you recommend and why? Also, what other solutions exist to carry out microarray-meta analysis starting from .CEL files and sample group data?

Cheers, Ankur Chakravarthy.

ADD COMMENTlink modified 6.2 years ago by boczniak767630 • written 6.2 years ago by Ankur90
gravatar for boczniak767
6.2 years ago by
boczniak767630 wrote:

I'd recommend you the second approach.
i) Normalize data from different studies separately - each lab and experiment has unique component (noise) added to data.

ii) Next, you could transform the data (again separately for each experiment) with Probability Of Expression (POE) transformation.
Package available here and here

iii) Then you can merge transformed datasets and analyse with statistical test. I've used multtest, but I have only two conditions and t-test sufficed

Other option - you could try packages RankProd or GeneMeta or metahdep
- I haven't tried it by myself (like above procedure) but there are examples of Affy analysis.


ADD COMMENTlink written 6.2 years ago by boczniak767630
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1109 users visited in the last hour