Search
Question: Statistical analysis of Spectral count data : edgeR, limma or others ?
2
gravatar for eleonoregravier
2.4 years ago by
France
eleonoregravier40 wrote:

Hi BioC community,

I have spectral count data (proteomics) to analyze. I have two groups (A ad B) to compare with 8 patients by group. Each sample was processed in duplicate, leading to 16 distinct patients in the study and 32 samples (16 by group).

As for the method to use to detect differentially expressed proteins between the two groups, I read a lot of things in the support site and in publications, in particular :

- Use the classical edgeR with TMM normalization or DEseq (Anders et al, nature protocol, 2013)

- Use edgeR with total spectral count normalisation (like in the msms package)

- Use limma on the log2 (count+1) with quantile normalization and eBayes(trend=TRUE, robust=TRUE) : voom for spectral counts

- Use edgeR and its quasi negative binomial approach : voom for spectral counts

- PLGEM package

- ...etc

As my data are count data, I wanted to use edgeR or DEseq package. But with the post of Gordon Smyth, 5 months ago ( voom for spectral counts ), I wanted to switch to limma with quantile normalization but I have no reference that I can cite.

Did you learn more about comparing the different approach on spectral count data ? Do you have recommendations to give ? Is there an article that compare the different methods and that I can cite please ?

Thanks a lot for your help,

Eléonore

ADD COMMENTlink modified 2.4 years ago by Gordon Smyth32k • written 2.4 years ago by eleonoregravier40
6
gravatar for Gordon Smyth
2.4 years ago by
Gordon Smyth32k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth32k wrote:

I don't have anything to add to my earlier post. As far as a reference, you could cite

  http://www.ncbi.nlm.nih.gov/pubmed/25605792

for limma or

  http://www.statsci.org/smyth/pubs/QLedgeRPreprint.pdf

for quasi-likelihood but we haven't published anything specifically about proteomics spectral data yet. It's too early for comparison papers I think.

ADD COMMENTlink modified 2.4 years ago • written 2.4 years ago by Gordon Smyth32k

Thanks a lot for your answer and the references,

Eleonore

ADD REPLYlink written 2.4 years ago by eleonoregravier40
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: 149 users visited in the last hour