Alternative packages to limma for analysis of microarray datasets
2
0
Entering edit mode
Maximilian ▴ 10
@maximilian-11742
Last seen 4.4 years ago

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

limma microarray deg • 2.2k views
ADD COMMENT
1
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY
1
Entering edit mode
svlachavas ▴ 780
@svlachavas-7225
Last seen 8 hours ago
Germany/Heidelberg/German Cancer Resear…

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 COMMENT
1
Entering edit mode

Other possibilities include multtest and rankprod.

ADD REPLY
0
Entering edit mode

Hello James,

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

ADD REPLY
1
Entering edit mode

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 REPLY
0
Entering edit mode
Bernd Klaus ▴ 600
@bernd-klaus-6281
Last seen 3.0 years ago
Germany

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 COMMENT

Login before adding your answer.

Traffic: 470 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