This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DeSeq2 "some values in assay are negative"

Hello everyone, I'm back again with RNA-seq data analysis and I encountered an error that seems to be common for the work that I'm doing

dds <- DESeqDataSetFromMatrix(countdata, coldata, design=~condition)

Error in DESeqDataSet(se, design = design, ignoreRank) : some values in assay are negative

I tried this solve in this way:

mat <- countdata[,-1]

rownames(mat) <- countdata[,1]

message: non-unique values when setting 'row.names': ‘0’, ‘1’, ‘2’, ‘3’, ‘4’, ‘5’, ‘6’, ‘7’, ‘8’, ‘9’, ‘10’, ‘11’, ‘12’, ‘13’, ‘14’, ‘15’, ‘16’, ‘17’, ‘18’, ‘19’, ‘20’, ‘21’, ‘22’, ‘23’, ‘24’, ‘25’, ‘26’, ‘27’, ‘28’, ‘29’, ‘30’, ‘31’, ‘32’, ‘33’, ‘34’, ‘35’, ‘36’, ‘37’, ‘38’, ‘39’, ‘40’, ‘41’, ‘42’, ‘43’, ‘44’, ‘45’, ‘46’, ‘47’, ‘48’, ‘49’, ‘50’, ‘51’, ‘52’, ‘53’, ‘54’, ‘55’, ‘56’, ‘57’, ‘58’, ‘59’, ‘60’, ‘61’, ‘62’, ‘63’, ‘64’, ‘65’, ‘66’, ‘67’, ‘68’, ‘69’, ‘70’, ‘71’, ‘72’, ‘73’, ‘74’, ‘75’, ‘76’, ‘77’, ‘78’, ‘79’, ‘80’, ‘81’, ‘82’, ‘83’, ‘84’, ‘85’, ‘86’, ‘87’, ‘88’, ‘89’, ‘90’, ‘91’, ‘92’, ‘93’, ‘94’, ‘95’, ‘96 [... truncated]duplicate 'row.names' are not allowed In addition: WarningError in .rowNamesDF<-(x, value = value) :

but I can't find duplicated rows...

I accept any suggestion and thank you for help

rnaseq
• 4,543 views • link
written by Morris_Chair Morris_Chair

head(countdata)

head(coldata)

• 0 views • link
written by ATpoint ATpoint

Hi ATpoint,

here is the result

 > head(countdata)
                   SRR5088824_starReadsPerGene.out.tab.x SRR5088825_starReadsPerGene.out.tab.x
ENSG00000000003.10                                    70                                    70
ENSG00000000005.5                                      0                                     0
ENSG00000000419.8                                     83                                    83
ENSG00000000457.9                                     62                                    62
ENSG00000000460.12                                    39                                    39
                   SRR5088834_starReadsPerGene.out.tab.x SRR5088835_starReadsPerGene.out.tab.x
ENSG00000000003.10                                    70                                    70
ENSG00000000005.5                                      0                                     0
ENSG00000000419.8                                     83                                    83
ENSG00000000457.9                                     62                                    62
ENSG00000000460.12                                    39                                    39
                   SRR5088842_starReadsPerGene.out.tab.x SRR5088843_starReadsPerGene.out.tab.x
ENSG00000000003.10                                    70                                    70
ENSG00000000005.5                                      0                                     0
ENSG00000000419.8                                     83                                    83
ENSG00000000457.9                                     62                                    62
ENSG00000000460.12                                    39                                    39
                   SRR5088844_starReadsPerGene.out.tab.x SRR5088859_starReadsPerGene.out.tab.x
ENSG00000000003.10                                    70                                    70
ENSG00000000005.5                                      0                                     0
ENSG00000000419.8                                     83                                    83
ENSG00000000457.9                                     62                                    62
ENSG00000000460.12                                    39                                    39
                   SRR5088860_starReadsPerGene.out.tab.x SRR5088861_starReadsPerGene.out.tab.x
ENSG00000000003.10                                    70                                    70
ENSG00000000005.5                                      0                                     0
ENSG00000000419.8                                     83                                    83
ENSG00000000457.9                                     62                                    62
ENSG00000000460.12                                    39                                    39


 getOption("max.print") -- omitted 1 rows ]


> > head(coldata)

   condition 
SRR5088824_starReadsPerGene.out.tab.x responsive
SRR5088825_starReadsPerGene.out.tab.x responsive
SRR5088834_starReadsPerGene.out.tab.x responsive
SRR5088835_starReadsPerGene.out.tab.x responsive
SRR5088842_starReadsPerGene.out.tab.x responsive
SRR5088843_starReadsPerGene.out.tab.x responsive

I shorted the countdata list

Thank you

• 0 views • link
written by Morris_Chair Morris_Chair

Please format data and code with the code button, not as a quote, it's messy and unreadable that way. The counts must be a numeric matrix. Rows being genes and columns being samples.

• 0 views • link
written by ATpoint ATpoint

even if I convert the data.frame as matrix doesn't help

countdata<-as.matrix(countdata)
dds <- DESeqDataSetFromMatrix(countdata, coldata, design=~condition)
>Error in DESeqDataSet(se, design = design, ignoreRank) : 
>some values in assay are negative
• 0 views • link
written by Morris_Chair Morris_Chair

Did you make sure that:

all(rownames(coldata) == colnames(count_matrix))

?

• 0 views • link
written by Marco Pannone • Marco Pannone

Hi Marco, It comes as "TRUE "

all(rownames(coldata) == colnames(countdata))
>[1] TRUE

thanks

• 0 views • link
written by Morris_Chair Morris_Chair

yes I did it but without success :(

• 0 views • link
written by Morris_Chair Morris_Chair

When you want to reassign rownames, try this way:

mat <- countdata[,-1]
names_countdata <- countdata[,1]
rownames(mat) <- make.names(names_countdata, unique=TRUE)
• 0 views • link
written by Marco Pannone • Marco Pannone

Hi Marco I think I found the error, the problem is on my merged files which duplicates values in the same boxes and does some more errors..

Thank you for your answers :)

• 0 views • link
written by Morris_Chair Morris_Chair

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Nebula 30x read length distribution
    written by berndmann • berndmann
    I used samtools view {path}/{sample_name}.bam -F 4 | cut -f 10 | perl -ne \'chomp;print length($_) . "\\n"\' | sort -n | uniq -c > …
  • How to change the number of seeds/predictions run by Alphafold
    written by umjpyw • umjpyw
    Hi everyone, I'm running alphafold 2, but I'm wanting the output to give me 100 predictions (and hence 100 structure I can put into PyMol) …
  • LOWESS Regression analysis in R
    written by georgians • georgians
    Hello Everyone, Below is the table each id column contains 7 visits and their value is in AC107c column. I want to use lowess function …
  • Removing adapters for assembly
    written by Ric Ric
    Hi, I ran fastp with this parameters `fastp --detect_adapter_for_pe --correction ... ` and it removed the adapters and trimmed the paired-end reads. Which read length …
  • DESeq2 and integers values
    written by Morris_Chair Morris_Chair
    Dear all, I have an experiment with two different conditions, two treated and two untreated samples from where I performed RNA sequencing. here is the …
  • BiomaRt: archive version of ensembl not working
    written by salamandra salamandra
    Hi, I want to access an archived version of biomart and select the human dataset. mart=useMart("ensembl_mart_92",host="http://apr2018.archive.ensembl.org") works but if I do: mart=useMart("ensembl_mart_92",host="http://apr2018.archive.ensembl.org", dataset = "hsapiens_gene_ensembl") …
  • DEseq matrix preparation
    written by Learner Learner
    I recently had a question related to having data as integer for DEseq . One data looks like this will have no problem > head(df) …
  • Mauvealigner Error
    written by HG HG
    <p>Dear All, I used MIRA for genome assembly and i got contigfile " Ecoli-1_out.unpadded.fasta". Now i want to do alignment with this contig with another …
  • Extracting Set Of Numbers From A File
    written by Cynthia • Cynthia
    <p>In a text file there are numbers. for eg:</p> <blockquote> <p>human 0 1 2 3 4 5 6 7 8 9 10 11 12 13 …
  • Synonmous And Non Synonmous Snps
    written by Payal Payal
    <p>Dear friends,</p> <p>I've compared the two sequences and i ve got the locations of SNPs between the two sequences. for eg:</p> <p><strong>Sequence 1</strong>: ---------GCCTGCTGCTGGCCGGGCGGGGGACGGGG------------CGGGACCGGAGCCGGAGCTGCGGGGCGCACCGGCTAGA</p> <p><strong>sequence …