Search
Question: Alternative packages to limma for analysis of microarray datasets
0
gravatar for Maximilian
10 weeks ago by
Maximilian10
Maximilian10 wrote:

Hello Bioconductor forum,

I used Bioconductor package limma for analysis of differentially expressed genes in this Agilent microarray dataset available under https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=+GSE59408.

Now my question is: are there any alternative packages to limma for analysis of differentially expressed genes in microarray datasets? I would like to compare results for the differentially expressed genes I found with limma. I did some research but only found alternative packages to limma for RNA-seq data and not for microarray data.

Kind regards,

Max

ADD COMMENTlink modified 9 weeks ago by Bernd Klaus430 • written 10 weeks ago by Maximilian10
1

By my understanding, the alternatives to limma for RNA-seq, such as edgeR and DESeq, were created because limma's statistical model was not (at that time) a good fit for count data. Do you have concerns about limma's statistical model for your data, or is there a specific statistical method you're looking to apply and you're hoping to find a package to do it?

ADD REPLYlink written 10 weeks ago by Ryan C. Thompson5.8k

Hello Ryan,

thank you for your answer. No, I don't have concerns about limma' statistical model for my data, the results for the DEGs with limma are in accord with the gene expression described in the original publication Prediction of antibiotic resistance by gene expression profiles available under https://www.nature.com/articles/ncomms6792.

However, the analysis with limma revealed some additional, possibly interesting genes in this dataset and I thought about using a different package/software/method to see if these results are reproducible.

ADD REPLYlink modified 10 weeks ago • written 10 weeks ago by Maximilian10
1
gravatar for svlachavas
10 weeks ago by
svlachavas540
Greece/Athens/National Hellenic Research Foundation
svlachavas540 wrote:

Dear Max,

personally i would not consider anything better than limma for the analysis of microarray data (not only-)-this has been proven in numerous scientific published papers and pipelines-rather than your approach of using alternative R packages, you could focus on the interpretation of your DE list. Does functional enrichment of your DE genes provides "sensible results" regarding your biological question ? Or for instance, have you performed any literature mining to search for any similar transcriptomic analyses that have answered "similar experimental designs ? Also in parallel, you could create some additional plots, like a heatmap to inspect the expression pattern of your DE genes (and lots of more...)

Nevertheless, i could for instance mention the samr R package and cyberT test as alternatives, but i clearly think that limma provides the best choise.

Best,

Efstathios

ADD COMMENTlink modified 10 weeks ago • written 10 weeks ago by svlachavas540
1

Other possibilities include multtest and rankprod.

ADD REPLYlink written 10 weeks ago by James W. MacDonald43k

Hello James,

thank you for your suggestions. I will try them out.

ADD REPLYlink modified 10 weeks ago • written 10 weeks ago by Maximilian10
1

Hello Efstathios,

thank you for your answer. As I said in my answer to Ryan, the results for the DEGs with limma are in accord with the gene expression described in the original publication. However, there were some additional, possibly interesting genes in the results and I wanted to see if these results are reprducible with other packages/software/methods. Thank you for your suggestions samr and cyberT.

ADD REPLYlink modified 10 weeks ago • written 10 weeks ago by Maximilian10
0
gravatar for Bernd Klaus
9 weeks ago by
Bernd Klaus430
Germany
Bernd Klaus430 wrote:

if you are interested in two group comparisons only, the st CRAN package implements various alternative variance shrinkage estimators, the CAT score even takes correlations between the genes into account:

 

http://strimmerlab.org/software/st/index.html

 

ADD COMMENTlink written 9 weeks ago by Bernd Klaus430
Please log in to add an answer.

Help
Access

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