Search
Question: Alternative packages to limma for analysis of microarray datasets
0
gravatar for Maximilian
4 months 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 4 months ago by Bernd Klaus440 • written 4 months 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 4 months ago by Ryan C. Thompson5.9k

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 4 months ago • written 4 months ago by Maximilian10
1
gravatar for svlachavas
4 months ago by
svlachavas560
Greece/Athens/National Hellenic Research Foundation
svlachavas560 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 4 months ago • written 4 months ago by svlachavas560
1

Other possibilities include multtest and rankprod.

ADD REPLYlink written 4 months ago by James W. MacDonald44k

Hello James,

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

ADD REPLYlink modified 4 months ago • written 4 months 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 4 months ago • written 4 months ago by Maximilian10
0
gravatar for Bernd Klaus
4 months ago by
Bernd Klaus440
Germany
Bernd Klaus440 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 4 months ago by Bernd Klaus440
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: 175 users visited in the last hour