getting error with DEseq2
1
1
Entering edit mode
Nemo ▴ 80
@nemo-7332
Last seen 6.3 years ago
India

Hello, 

I am getting an error like this 

Error in DESeqDataSet(se, design = design, ignoreRank) : 
  some values in assay are not integers

Based on previous questions, I tried to solve it but I could not. So let me tell you how to imported the data etc 

countdata<-read.table("data.txt",header=TRUE,row.names=1)

conds<-as.factor(c("Treated","Control"))

coldata <- data.frame(row.names=colnames(countdata), conds)

dds=DESeqDataSetFromMatrix(countData=countdata,colData=coldata,design=~conds)

The following are more information about countdata 

 

deseq2 • 21k views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 1 hour ago
United States

The error says some of your values are not integers, and your data.frame looks like this:

> head(countdata)
         Control    Treated
gene1 0.00842815 0.01090013
gene2 0.25059516 0.98230798
gene3 0.29728284 0.15692372
gene4 0.30918591 0.36842593
gene5 0.31040477 0.47786814
gene6 0.31833954 0.40623309

Which makes sense, because none of those values are integers! If you want to use DESeq2 to analyze your data, you need to get count data which will be strictly integer values.

ADD COMMENT
0
Entering edit mode

@James W. MacDonald how should i get that ? any suggestion 

ADD REPLY
0
Entering edit mode

I don't know. Where did you get these data to begin with? I would imagine you should be talking to whomever you got the data from.

ADD REPLY
0
Entering edit mode

@James W. MacDonald I think, the answer is very easy , I should just get the data normalised , Not a rocket science ! 

ADD REPLY
0
Entering edit mode

No. The data you input has to be integers. Whatever you think normalization is going to do, rest assured it won't. Do you know what an integer is? Maybe that's the problem. If you don't understand that term, then we will never get anywhere.

ADD REPLY
0
Entering edit mode

@James W. MacDonald  of course i do. as I said , if I do make them integer as they are now , i'll lose a lot of info. the best is to normalise them and then use the function "as.integer" in R . 

ADD REPLY
0
Entering edit mode

James is correct.

Do not provide normalized data to DESeq2. 

You need to generate counts using one of the methods described in the documentation:

vignette("DESeq2")

On the very first page of the vignette, after the table of contents, you will find a section called "Why un-normalized counts?" Read this.

ADD REPLY

Login before adding your answer.

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