Hi,

 I wanna keep the genes equal or lower than 10 counts in at least one sample within Deseq2. is the following correct? Doesn't "sum" in "rowSums" indicate adding function?

keep <- rowSums(counts(ddsHTSeq) >= 10) >= 1

Thanks in advance,

N. hosseinzadeh

Yes. You read that from the inside out.

counts(ddsHTSeq) >= 10

returns a boolean matrix with TRUE for any gene with counts >= 10

rowSums() >= 1

computes the number of TRUE values you have, per gene (TRUE/FALSE are coerced to 1/0) and tests that you have >= 1.