Missing values - metaArray, GeneMeta, RankProd, metahdep
2
0
Entering edit mode
boczniak767 ▴ 720
@maciej-jonczyk-3945
Last seen 1 day ago
Poland
Dear List Members, I'd like to perform a metaanalysis for few (minimum three) datasets. Some of them are one channel (Affymetrix), some two-colour data. Obviously, when I join, say three studies I will have big dataset in which some data will be missing. Just because some probes will be present in two experiments and absent in third experiment. I'd prefer not to use *only* probes which are present on all arrays, but also this which are present on two arrays. My *question* is: which package (metaArray, GeneMeta, RankProd, metahdep) could handle metaanalysis with missing data - in tutorials/manuals it is not stated if missing data is allowed. Although metaArray cites MergeMaid package as a example software for merging data with non-overlapping genes. Second *question* which of the above package will best handle dataset from merged two-colour and affymetrix studies? I'd also be grateful for directions to additional tutorials, materials, etc. Best Regards, Maciej Jończyk Maciej Jończyk, MSc Department of Plant Molecular Ecophysiology Institute of Plant Experimental Biology Faculty of Biology, University of Warsaw 02-096 Warszawa, Miecznikowa 1
RankProd GeneMeta metaArray RankProd GeneMeta metaArray • 1.4k views
ADD COMMENT
0
Entering edit mode
John Stevens ▴ 30
@john-stevens-3739
Last seen 9.6 years ago
Hi Maciej, My metahdep package doesn't require all datasets to come from the same array version -- so while the package's tutorial vignette doesn't explicitly say anything about missing values, the example there does allow for missing values in this sense (some genes not represented in all datasets). As to your second question, if you're combining studies with fundamentally different platforms (Affy for some, two-color for others, for example), you'll need to be careful to define an effect size (or some measure of differential expression) that means the same thing on all platforms, and for which you can calculate a meaningful standard error. I think trans-platform cases like this lead to the appeal of alternative scale methods such as those in the metaArray package (which uses the 'probability of expression' or POE scale). Regards, John Stevens ________________________________________ From: bioconductor-bounces@r-project.org [bioconductor- bounces@r-project.org] on behalf of mjonczyk@biol.uw.edu.pl [mjonczyk@biol.uw.edu.pl] Sent: Tuesday, August 09, 2011 1:45 PM To: bioconductor at r-project.org Subject: [BioC] Missing values - metaArray, GeneMeta, RankProd, metahdep Dear List Members, I'd like to perform a metaanalysis for few (minimum three) datasets. Some of them are one channel (Affymetrix), some two-colour data. Obviously, when I join, say three studies I will have big dataset in which some data will be missing. Just because some probes will be present in two experiments and absent in third experiment. I'd prefer not to use *only* probes which are present on all arrays, but also this which are present on two arrays. My *question* is: which package (metaArray, GeneMeta, RankProd, metahdep) could handle metaanalysis with missing data - in tutorials/manuals it is not stated if missing data is allowed. Although metaArray cites MergeMaid package as a example software for merging data with non-overlapping genes. Second *question* which of the above package will best handle dataset from merged two-colour and affymetrix studies? I'd also be grateful for directions to additional tutorials, materials, etc. Best Regards, Maciej Jończyk Maciej Jończyk, MSc Department of Plant Molecular Ecophysiology Institute of Plant Experimental Biology Faculty of Biology, University of Warsaw 02-096 Warszawa, Miecznikowa 1 _______________________________________________ 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
ADD COMMENT
0
Entering edit mode
boczniak767 ▴ 720
@maciej-jonczyk-3945
Last seen 1 day ago
Poland
Hi John, thanks for reply. I've performed a part of example analysis from metaArray, I haven't succeeded to perform analysis on a dataset from three partially overlapping experiments. Zscore function gave results only for the common part. If I found a way to transform various experiments to effect-size scale (probably with inverse-variance method as recommended in Ramasamy et al. 2008 - very useful article) or something similar I will try metahdep. Best Regards, Maciej Jończyk > Hi Maciej, > > My metahdep package doesn't require all datasets to come from the same array version -- so while the package's tutorial vignette doesn't explicitly say anything about missing values, the example there does allow for missing values in this sense (some genes not represented in all datasets). > > As to your second question, if you're combining studies with fundamentally different platforms (Affy for some, two-color for others, for example), you'll need to be careful to define an effect size (or some measure of differential expression) that means the same thing on all platforms, and for which you can calculate a meaningful standard error. I think trans-platform cases like this lead to the appeal of alternative scale methods such as those in the metaArray package (which uses the 'probability of expression' or POE scale). > > Regards, > > John Stevens > > > > ________________________________________ > From: bioconductor-bounces at r-project.org [bioconductor-bounces at r-project.org] on behalf of mjonczyk at biol.uw.edu.pl [mjonczyk at biol.uw.edu.pl] > Sent: Tuesday, August 09, 2011 1:45 PM > To: bioconductor at r-project.org > Subject: [BioC] Missing values - metaArray, GeneMeta, RankProd, metahdep > > Dear List Members, > > I'd like to perform a metaanalysis for few (minimum three) datasets. > Some of them are one channel (Affymetrix), some two-colour data. > Obviously, when I join, say three studies I will have big dataset in > which some data will be missing. Just because some probes will be present > in two experiments and absent in third experiment. I'd prefer not to > use *only* probes which are present on all arrays, but also this which > are present on two arrays. > > My *question* is: which package (metaArray, GeneMeta, RankProd, metahdep) > could handle metaanalysis with missing data - in tutorials/manuals it is > not stated if missing data is allowed. Although metaArray cites MergeMaid > package as a example software for merging data with non-overlapping genes. > > Second *question* which of the above package will best handle dataset > from merged two-colour and affymetrix studies? > > I'd also be grateful for directions to additional tutorials, materials, etc. > > Best Regards, > Maciej Jończyk > > Maciej Jończyk, MSc > Department of Plant Molecular Ecophysiology > Institute of Plant Experimental Biology > Faculty of Biology, University of Warsaw > 02-096 Warszawa, Miecznikowa 1
ADD COMMENT

Login before adding your answer.

Traffic: 702 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6