Pre-filtering DESeqDataSetFromMatrix input
1
0
Entering edit mode
rbronste ▴ 60
@rbronste-12189
Last seen 5.0 years ago

A little confused about this, wanted too filter out lines in the matrix with read counts lower than 10 to speed up DESeq2 processing using the following:

dds <- dds[ rowSums(counts(dds)) > 10, ]

However just as before the filtering, after the filtering countData looks like this:

head(countData)
  MBV1 MBV2 MBV3 FBV1 FBV2 FBV3 MBE7 MBE8 MBE9 FBE1 FBE2 FBE3
1    1    1    1    1    1    1    1   66    1    1   50   34
2   11    1    1    1    1    1    6   98    1   11  100    1
3    1    1    1    1    1    1    1    1    1  116  108    1
4    1    1   22    2   84    1    1    4    1   64    1   40
5    1    1   18   74   74    1  102    1  126   22    1    1
6    1    1    1    1   44    1    1    1  122    1    1    1

Wondering where I am going wrong? Thank you.

deseq2 differential binding analysis • 890 views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 1 day ago
United States

Check your parentheses. You are filtering on the rowSums of the matrix, not the number of samples with counts > 10.

ADD COMMENT
0
Entering edit mode

Thank you! Stupid mistake on my part for sure...

ADD REPLY

Login before adding your answer.

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