Question: DESeq2/tximport 'duplicate rownames not allowed' error
0
gravatar for ctl
14 months ago by
ctl0
London/United Kingdom
ctl0 wrote:

Hi, 

I have used Salmon to quantify some 126 samples using the mouse ensembl reference cDNA transcriptome. Then I tried to obtained 

abundance values using tximport as follows:

txi.salmon <- tximport(files, type = "salmon", tx2gene = tx2gene, countsFromAbundance = "lengthScaledTPM", ignoreTxVersion = F)

names(txi.salmon)
head(txi.salmon$counts)

###Create a matrix of abundance (a.k.a lengthScaledTPM) from the tximport object and import it as a text dataframe

cts <- txi.salmon$counts
abn <- txi.salmon$abundance
samples_rc <- data.frame(cts)
samples_tpm <- data.frame(abn)

This has worked without errors many times. 

Then I wanted to used these values with DESeq2, but I've gotten the following error: 

ddsTxi <- DESeqDataSetFromTximport(txi.salmon, colData = metadata, design = ~condition + day + strain + centre)
Error in `rownames<-`(`*tmp*`, value = colnames(countData)) : 
  duplicate rownames not allowed

This is the first time this occurs and I have no idea what may be wrong, since I have used the same step many times successfully. 

Now, I understand that this is complaining about repeated gene names in the rows, but this has never happened before. 

 Can you give me an idea of what may have gone wrong?

 

Thanks 

 

 

 

deseq2 tximport • 436 views
ADD COMMENTlink written 14 months ago by ctl0
Answer: C: DESeq2/tximport 'duplicate rownames not allowed' error
0
gravatar for ctl
14 months ago by
ctl0
London/United Kingdom
ctl0 wrote:

Please ignore this. I found the mistake. 

My metadata table had some repeated samples, after removing them everything worked just fine. 

Thanks! 

ADD COMMENTlink written 14 months ago by ctl0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 303 users visited in the last hour