I have quantified the RNA-seq samples by Salmon. 2 groups are wild-type and 4 groups are Dnmt2 knocked-out. I've put all in one dataset for DE analysis. box plot of their normalized counts shows the median of knocked-out samples are zero, and maximum of 20 reads are assigned to each transcript. now to perform DE by DESeq2 I have 2 questions: 1- whether zero values should be deleted before? 2- to do DE what minimum of counts has to be selected?
thanks in advance
Sorry, I think, I made a mistake cause I've calculated log2(1+counts), and made a box plot. it's y axis is between 0 to 20. what does it mean?
log2 of 20 is typical. No problem with this data, or anything you've described above.
Perhaps you can look at the vignette and workflow so you get an idea of what typical RNA-seq count datasets look like.
Excuse me, to make plotMA with res following DESeq2 vignette, it makes a plot based on expression log ratio and log expression, not LFC and mean of normalized counts, plotMA(ddsTxi, alpha= 0.1, main = "", xlab = "mean of normalized counts", ylim, mle = TRUE) I face with an error: Error in as.vector(x) : no method for coercing this S4 class to a vector. whether it should be converted to data.frame?
What is
class(ddsTxi)
?If it is a
DESeqDataSet
orDESeqResults
object it should work.yes, it's a DESeqDataSet. but it doesn't work!
Can you try
DESeq2::plotMA()
. Maybe you are using another package that masks our plotMA method.exactly, it works. thanks alot
Excuse me, I would appreciate if help me to making a heatmap of the count matrix which I performed according to DESeq2 vignette,
select <- order(rowMeans(counts(ddstxi,normalized=TRUE)), decreasing=TRUE)[1:20]
df <- as.data.frame(colData(ddstxi)[,c("group")])
pheatmap(assay(ntd)[select,], clusterrows=FALSE, showrownames=FALSE, clustercols=FALSE, annotationcol=df)
but it face with this error,
Error in check.length("fill") : 'gpar' element 'fill' must not be length 0
Not sure of that error. It’s coming from pheatmap not DESeq2 so check what you are inputting and check the help files from that package.