Error in DESeqDataSet(se, design = design, ignoreRank) : all variables in design formula must be columns in colData
1
0
Entering edit mode
yueli7 ▴ 20
@yueli7-8401
Last seen 3.5 years ago
China

Hello,

How can I figure out this?

Thanks in advance!

> head(huang_group_list)
                  x
1    549_repair.bam
2      6_repair.bam
3  6plus_repair.bam
4     14_repair.bam
5 14plus_repair.bam

 

> colData
                                   x
X549_repair.bam       549_repair.bam
X6_repair.bam           6_repair.bam
X6plus_repair.bam   6plus_repair.bam
X14_repair.bam         14_repair.bam
X14plus_repair.bam 14plus_repair.bam

 

> dds <- DESeqDataSetFromMatrix(countData = huang_exprSet, colData = colData, design = ~ huang_group_list)
Error in DESeqDataSet(se, design = design, ignoreRank) :
  all variables in design formula must be columns in colData


 


 

 

Actually, I try to paste

 

> head(huang_exprSet)
                  X549_repair.bam X6_repair.bam X6plus_repair.bam
ENSG00000223972.5               2             0                 0
 

 X14_repair.bam X14plus_repair.bam
ENSG00000223972.5              1                  0


 

          


 

It always said:  Language "af" is not one of the supported languages ['en']!


 

deseq2 • 7.3k views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 13 hours ago
United States

The error says what the problem is:

Error in DESeqDataSet(se, design = design, ignoreRank) :
  all variables in design formula must be columns in colData

So you have to have a column in your colData object called 'huang_group_list'. And given that the underscores won't pass R's make.names filter, you might now have huang.group.list as the column name.

ADD COMMENT
1
Entering edit mode

Thanks,

> colData
             group_list
treated1fb      treated
treated2fb      treated
treated3fb      treated
untreated1fb  untreated
untreated2fb  untreated
> dds <- DESeqDataSetFromMatrix(countData = exprSet, colData = colData, design = ~ group_list)


 

ADD REPLY

Login before adding your answer.

Traffic: 691 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