How to use filterbyExpr with DESeq2?
Hello I am trying to use edgeR's filterbyExpr with DESeq2. I am running into an issue where I cannot use filterByExpr with DESeq2, unless I do not specify keep.lib.sizes = FALSE (which is recommended by edgeR) how do I use filterByExpr on my DESeqDataSet object.
dds <- DESeqDataSetFromTximport(txi=txi.kallisto.tsv, colData=meta, design= ~1)
# Filter by expression (Get rid of low variance )
keep.exprs <- filterByExpr(dds)
dds_pre_filt<- dds[keep.exprs,]
dim(dds_pre_filt)
dds_deseq <- DESeq(dds_pre_filt) #
norm_dds<-vst(dds_deseq) # vst normalization
What I want is to be able to to is:
dds_pre_filt<- dds[keep.exprs,, keep.lib.sizes = FALSE]
and then precede with normalization
• 2,313 views
•
link
updated
by
Ram
•
written
by
biotrekker •
•
written
by
biotrekker •
1 answer
You cannot run it directly on a DESeqDataSet (dds) but must run it on counts(dds). The libsize option is irrelevant here, it only applies if you run it on a DGEList from edgeR.
dds <- DESeqDataSetFromTximport(txi=txi.kallisto.tsv, colData=meta, design= ~your_design)
keep <- filterByExpr(counts(dds))
dds <- DESeq(dds[keep,])
Log in to answer this question.
More posts like this
-
What is the right choice of normalization and preprocessing before WGCNA using RNAseq data in a coh…
written by sukeshinik5 •
-
How to get normalized counts corrected for variables in the design matrix of DESeq2 as an input for…
written by Sara •
-
error while plot data with deSeq2
written by bioinformaticssrm2011
-
How to subset the output of tximport in DE analysis using DEseq2
written by tianshenbio
-
tximport output should be integer values, IS NOT.
written by c_u
-
adding batch variable in colData(dds)
written by Mozart
-
How to normalize GTEX gene counts with DESeq2?
written by kakukeshi
-
Error in running DESeqDataSetFromTximport
written by ankur •
-
DESeq2 object problem with padj
written by pennakiza
-
Accessing data in within a DESeqDataSet class from DESeq2 library
written by irritable_phd_syndrome •
