Filtering of RNA-seq data
Hi! I'm going to filter my RNA-seq count data, and I'm wondering how I remove rows where only one sample has a count and the rest of the sample has 0 counts. I'm using DESeq and have only found how to remove sums of counts.
Thanks!
• 1,231 views
•
link
updated
by
Papyrus
•
written
by
BioQueen •
•
written
by
BioQueen •
1 answer
For example, this selects rows with all 0s, or all 0s except 1 sample. You can change the number to be more or less stringent:
filter <- apply(counts(deseq.object),1,function(x){
sum(x==0) >= (ncol(deseq.object) - 1)
})
# Now filter these rows
deseq.object2 <- deseq.object[!filter,]
Log in to answer this question.
More posts like this
-
Limma with already normalised data
written by BioQueen •
-
Batch effect removal
written by BioQueen •
-
Merging of RNA-seq datasets from different studies
written by BioQueen •
-
Batch effects
written by BioQueen •
-
Difference between dataset analysis
written by BioQueen •
-
Omics data collection and merging of datasets from different studies
written by BioQueen •
-
remove Batch Effect between several **RNA-seq** studies (no replication)
written by snstab •
-
Millions of Reads Needed for Differential Gene Expression
written by contact •
-
ssGSEA on RNA-Seq data from TCGA
written by oriolebaltimore
-
Fitting Distribution Of The Gene Expression Levels Of Rna Seq Data
written by k.nirmalraman
