This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error in my DESeq script

Could someone please help me troubleshoot this line?

dds <- DESeqDataSetFromMatrix( countData, DataFrame(condition), ~ condition)

I get the error:

Error in if (all(var == var[1])) { : missing value where TRUE/FALSE needed

rna-seq r deseq
• 2,080 views • link
written by sandKings • sandKings

Try removing transcripts with 0 expression value from your dataset.

• 0 views • link
written by Arup Ghosh Arup Ghosh

You probably get NAs somewhere in the building process so var==var[1] is neither TRUE nor FALSE. Do you have NAs in the countData matrix?

• 0 views • link
written by Asaf Asaf

While typing the response to your question, I figured out what I did wrong. The csv file which had my samples and treatments outlined, was for a different cell lineage. I accidentally copied the same file in 2 different folders. So, I didn't have the matching names in the countData. Thanks!!

• 0 views • link
written by sandKings • sandKings

It's called rubber-duck debugging. You take a rubber duck, put it next to your screen and explain to it what you did. Most of the time you'll catch the bug. Happy to be a rubber duck. For the next time, try to set up a start to end pipeline without human intervention, these types of errors can have a major effect.

• 0 views • link
written by Asaf Asaf

hahaha! rubber-duck debugging sounds very promising and discreet! Thanks!

• 0 views • link
written by sandKings • sandKings

0 answers

No answers yet.

Log in to answer this question.

More posts like this