Question: DESeq2 with proteomics data
0
gravatar for shyu
13 months ago by
shyu0
shyu0 wrote:

Dear all,

Recently, I am using DESeq2 for differential analysis of proteomics data. But I got an error:

invalid class “DESeqDataSet” object: NA values are not allowed in the count matrix

I checked my input table, it contains no "NA".

Afterwards, I found that somebody had a similar error as mine ( DESeq2 Error NA values ). Based on the response in this post, I also found out my input table contains some values which are higher than 10 billion. Thus, I do a normalization to reduce the number of values. Then, the error is gone, and I got the results.

I would like to confirm that does DESeq2 not accept the values which is extremely high? But this kind of high value is quite common in the intensity of proteomics.

Thank you,

Silas

deseq2 • 592 views
ADD COMMENTlink modified 13 months ago by sven.schenk10 • written 13 months ago by shyu0
Answer: DESeq2 with proteomics data
0
gravatar for Michael Love
13 months ago by
Michael Love25k
United States
Michael Love25k wrote:

We do not have support in DESeq2 for proteomics data.

There have been threads on this before. It's not simple count data as I understand (depending on actually what assay is performed...).

ADD COMMENTlink written 13 months ago by Michael Love25k
Answer: DESeq2 with proteomics data
0
gravatar for sven.schenk
13 months ago by
sven.schenk10
Vienna
sven.schenk10 wrote:

We also tried using DESeq on proteomics data, and this gave extrmely high numbers of rejections (of the null), which didn`t seem to be realistic.

The difference between RNASeq and proteomics data is that protein intensities are continuous data and not simple counts as are RNASeq data. This, I guess poses a problem for DESeq and the normalisation algorithm.

For proteomics data I`d recomment to use either LIMMA or the ROTS package, both of which perform similar and, importantly, can deal with missing data points.

 

ADD COMMENTlink modified 13 months ago • written 13 months ago by sven.schenk10
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 16.09
Traffic: 171 users visited in the last hour