Hi! Am a complete beginner to using DESeq2. I have read the vignette and am encountering an error when I am trying to create a DESeqDataSet from the count data. I am getting an error that some of the values in the assay are negative. However, I have scanned the data set and there are no negative integers. Also i removed the genes with zero values from my data set and started again but still getting the same error. Would be grateful if some one can help me here.
Code should be placed in three backticks as shown below
``dds <-DESeqDataSetFromMatrix(countData=countData,colData=colData, design= genotype+BDQ_treatment)enter code here
Error in DESeqDataSet(se, design = design, ignoreRank) : some values in assay are negative
sessionInfo( )
```
In order to back up your assertion that there are no negative values, you should provide evidence. For example
The first test is what is giving the error, and I would be quite surprised if it returns
FALSE
on a direct test and thenTRUE
once placed within aSummarizedExperiment
.