This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DESeq2 error message on R3.3

Hi, I ran the following script on R version 3.2 and everything was fine:

conds<- c("kd1","kd1","kd2","kd2","GFP","GFP")

Design<-data.frame(condition=conds,
                   row.names=paste(conds,rep(c(1,2),3),
                   sep="_"))
dds <- DESeqDataSetFromMatrix(countData = assay(se),
                             colData = Design,
                             design = ~ condition)

however, when I updated R into 3.3 I start getting the following error message:

Error in FUN(X[[i]], ...) : 
  assay colnames() must be NULL or equal colData rownames()

I tried to follow many solutions but non of them worked .

can anyone help

thanks much

rna-seq deseq2
• 3,288 views • link
written by ta_awwad ta_awwad

1 answer

So what's colnames(se)?

• 1 views • link
written by Devon Ryan Devon Ryan 
   >  colnames(se)
[1] "kd1_1.bam" "kd1_2.bam" "kd2_1.bam"  "kd2_2.bam"  "GFP_1.bam"  
[6] "GFP_2.bam"

whereas, 

> Design
                  condition
kd1_1          kd1
kd1_2          kd1
kd2_1          kd2
kd2_2          kd2
GFP_1        GFP
GFP_2        GFP
• 1 views • link
written by ta_awwad ta_awwad

There's your problem then.

• 1 views • link
written by Devon Ryan Devon Ryan

got it .. thanks much

• 1 views • link
written by ta_awwad ta_awwad

thanks I solved the problem

• 1 views • link
written by ta_awwad ta_awwad

Log in to answer this question.

More posts like this