Error with DEseq2: DESeqDataSetFromMatrix. error indicating row names of the colData are not in the same order as the colnames of the countData
1
0
Entering edit mode
lara • 0
@lara-18927
Last seen 6.9 years ago

I am receiving an error that "row names of the colData are not the same as in the countData.  But I can't find any discrepancy. 

thank you

 

 

> count_tab <- read.table("ASVs_counts.txt", header=T, row.names=1, check.names=F)
> tax_tab <- as.matrix(read.table("ASVs_taxonomy.txt", header=T, row.names=1, check.names=F, na.strings="", sep="\t"))
> sample_info_tab <- read.table("sample_info_milos.txt", header=T, row.names=1, check.names=F)

 

> deseq_counts <- DESeqDataSetFromMatrix(count_tab, colData = sample_info_tab, design = ~type) 
Error in DESeqDataSetFromMatrix(count_tab, colData = sample_info_tab,  : 
  rownames of the colData:
   0.3.6.O.DNA,0.3.6.C,0.3.6.C.DNA,0.3.6.O,0.9.6.C.2.DNA,0.9.6.C.2,100.3.24.C.DNA,100.3.24.C,100.3.24.O.DNA,100.3.24.O,100.3.6.C.DNA,100.3.6.C,100.3.6.O.2.DNA,100.3.6.O.2,100.9.6.C.DNA,100.9.6.C,100.9.6.O.DNA,100.9.6.O,50.3.24.O.DNA,50.3.24.O,50.3.6.C.DNA,50.3.6.C,50.3.6.O.DNA,50.3.6.O,50.9.6.O.DNA,50.9.6.O 
  are not in the same order as the colnames of the countData:
   0.3.6.C.DNA,0.3.6.C,0.3.6.O.DNA,0.3.6.O,0.9.6.C.2.DNA,0.9.6.C.2,100.3.24.C.DNA,100.3.24.C,100.3.24.O.DNA,100.3.24.O,100.3.6.C.DNA,100.3.6.C,100.3.6.O.2.DNA,100.3.6.O.2,100.9.6.C.DNA,100.9.6.C,100.9.6.O.DNA,100.9.6.O,50.3.24.O.DNA,50.3.24.O,50.3.6.C.DNA,50.3.6.C,50.3.6.O.DNA,50.3.6.O,50.9.6.O.DNA,50.9.6.O

 

deseq2 • 2.5k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 2 days ago
United States

I see a difference.

You can use == to compare two vectors in R.

ADD COMMENT

Login before adding your answer.

Traffic: 792 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6