DESeq2 with proteomics data
Entering edit mode
shyu ▴ 10
Last seen 4.7 years ago

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,


deseq2 • 3.3k views
Entering edit mode
Last seen 1 day ago
United States

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...).

Entering edit mode
sven.schenk ▴ 20
Last seen 3.9 years ago

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.



Login before adding your answer.

Traffic: 437 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6