Search
Question: Pre-filtering DESeqDataSetFromMatrix input
0
gravatar for rbronste
11 weeks ago by
rbronste50
rbronste50 wrote:

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.

ADD COMMENTlink modified 11 weeks ago by James W. MacDonald45k • written 11 weeks ago by rbronste50
1
gravatar for James W. MacDonald
11 weeks ago by
United States
James W. MacDonald45k wrote:

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

ADD COMMENTlink written 11 weeks ago by James W. MacDonald45k

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

ADD REPLYlink written 11 weeks ago by rbronste50
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: 191 users visited in the last hour