DESEQ object filter
1
0
Entering edit mode
Bine ▴ 40
@bine-23912
Last seen 3 months ago
UK

Hey guys,

Is this correct? I want to reduce my countData (cts):

The aim is to get at least 8 reads in 20% of the samples (which is 50 samples):

keep<- rowSums(cts>8)>=50

cts <-cts[keep,]

Thank you, Bine

DESeq2 • 864 views
ADD COMMENT
2
Entering edit mode
Kevin Blighe ★ 3.9k
@kevin
Last seen 1 day ago
Republic of Ireland

Hi, this is not quite a DESeq2 question, but, nevertheless, yes, this code will return a boolean vector that represents those genes (rows) that have at least 8 read counts in at least 50 samples:

rowSums(cts>=8) >= 50

Kevin

ADD COMMENT
1
Entering edit mode

Thank you very much for confirming!

ADD REPLY

Login before adding your answer.

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