Question: "Duplicate rownames not allowed" when the rownames are not duplicated
gravatar for Kristin Muench
5 months ago by
United States
Kristin Muench380 wrote:


I'm getting an error from DESeq2 that insists I have duplicates somewhere in my data, but I don't think I do.

I loaded transcript data from kallisto using tximport(). I wanted HUGO gene-wise counts, not transcript counts, so I aggregated these and fed the result into a DESeq object using DESeqDataSetFromMatrix. I double-checked to make sure that neither the rownames or column names in this count matrix were duplicated, but DESeq is still throwing an error.

Any ideas as to what might be the source of the issue, or what exactly DESeq is pointing to when it says "names(x)"?

Thanks for your help!


# pull count data out of tximport-created data frame for manipulation
hugoCrunch <- data.frame(assay(dds_kal))

# use my own function to add HUGO names using bioMart
hugoCrunch <- addHUGOnames(hugoCrunch, 'hsapiens_gene_ensembl', 'ensembl_transcript_id')

## add back into original data frame
hugoCrunch_data <- hugoCrunch[,-c(52:54)]
hugoCrunch_names <- hugoCrunch$Gene # these are the HUGO genes produced with addHUGOnames. There are duplicates here

# add genes together
hugoCrunch_crunched <- aggregate(hugoCrunch_data[], by = list(hugoCrunch_names), FUN = sum) 
row.names(hugoCrunch_crunched) <- hugoCrunch_crunched[,1]
 # row names are now HUGO genes, which are not duplicated

# remove row names and the NA gene counts
hugoCrunch_crunched <- hugoCrunch_crunched[-1,-1] 

# ensure everything is unique
rownames(hugoCrunch_crunched) = make.names(rownames(hugoCrunch_crunched), unique=TRUE)
colnames(hugoCrunch_crunched) = make.names(colnames(hugoCrunch_crunched), unique=TRUE)

# put data back into a DESeq object
dds_kal_agg <- DESeqDataSetFromMatrix(countData = hugoCrunch_crunched,
                              colData = sampleTable_plus,
                              design= ~ Genotype)

# do sampleTable_plus IDs have duplicates in the first row (DESeq ID row)?


# DESeq
dds_deseq <- DESeq(dds_kal_agg)

Error in rownames<-(*tmp*, value = names(x)) : duplicate rownames not allowed


R version 3.5.1 (2018-07-02)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS High Sierra 10.13.6

Matrix products: default
BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib

[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] parallel  stats4    stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] gdtools_0.1.7               ggrepel_0.8.0               ggplot2_3.0.0               limma_3.36.5                pheatmap_1.0.10            
 [6] biomaRt_2.36.1              DESeqAid_0.2                DESeq2_1.20.0               SummarizedExperiment_1.10.1 DelayedArray_0.6.6         
[11] BiocParallel_1.14.2         matrixStats_0.54.0          Biobase_2.40.0              GenomicRanges_1.32.7        GenomeInfoDb_1.16.0        
[16] IRanges_2.14.12             S4Vectors_0.18.3            BiocGenerics_0.26.0        

loaded via a namespace (and not attached):
 [1] httr_1.3.1             bit64_0.9-7            splines_3.5.1          Formula_1.2-3          assertthat_0.2.0       latticeExtra_0.6-28   
 [7] blob_1.1.1             GenomeInfoDbData_1.1.0 progress_1.2.0         yaml_2.2.0             pillar_1.3.0           RSQLite_2.1.1         
[13] backports_1.1.2        lattice_0.20-35        glue_1.3.0             digest_0.6.17          RColorBrewer_1.1-2     XVector_0.20.0        
[19] checkmate_1.8.5        colorspace_1.3-2       htmltools_0.3.6        Matrix_1.2-14          plyr_1.8.4             XML_3.98-1.16         
[25] pkgconfig_2.0.2        genefilter_1.62.0      zlibbioc_1.26.0        purrr_0.2.5            xtable_1.8-3           scales_1.0.0          
[31] svglite_1.2.1          htmlTable_1.12         tibble_1.4.2           annotate_1.58.0        withr_2.1.2            nnet_7.3-12           
[37] lazyeval_0.2.1         survival_2.42-6        magrittr_1.5           crayon_1.3.4           memoise_1.1.0          foreign_0.8-71        
[43] prettyunits_1.0.2      tools_3.5.1            data.table_1.11.4      hms_0.4.2              stringr_1.3.1          locfit_1.5-9.1        
[49] munsell_0.5.0          cluster_2.0.7-1        AnnotationDbi_1.42.1   bindrcpp_0.2.2         compiler_3.5.1         rlang_0.2.2           
[55] grid_3.5.1             RCurl_1.95-4.11        rstudioapi_0.7         htmlwidgets_1.3        labeling_0.3           bitops_1.0-6          
[61] base64enc_0.1-3        gtable_0.2.0           curl_3.2               DBI_1.0.0              R6_2.2.2               gridExtra_2.3         
[67] knitr_1.20             dplyr_0.7.6            bit_1.1-14             bindr_0.1.1            Hmisc_4.1-1            stringi_1.2.4         
[73] Rcpp_0.12.19           geneplotter_1.58.0     rpart_4.1-13           acepack_1.4.1          tidyselect_0.2.4


Ah, of course I found a relevant post as soon as I post my own:

My issue was solved once I changed the DESeq line to:

# DESeq
dds_deseq <- DESeq(dds_kal_agg, minRep=Inf)

I'm not entirely sure how this worked. When I follow the code on this website, it does look like my data have a few outliers:

outliers in data

deseq2 • 503 views
ADD COMMENTlink modified 5 months ago by Michael Love1.8k • written 5 months ago by Kristin Muench380
gravatar for Michael Love
5 months ago by
Michael Love1.8k
United States
Michael Love1.8k wrote:

To propagate the answer from Bioconductor support site back here:

The problem is that you have duplicate columns of colData(dds), which breaks some code where replaceOutliers adds a column to colData(dds) and adds some metadata about that column.


So a solution is to only have unique column names for colData(dds), which is probably a good idea anyway.

ADD COMMENTlink modified 5 months ago • written 5 months ago by Michael Love1.8k

Thank you for taking a look at this!

ADD REPLYlink written 5 months ago by Kristin Muench380

Thanks for checking in and following up, Michael!

ADD REPLYlink written 5 months ago by Kevin Blighe39k
Please log in to add an answer.


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