Differential expression analysis using DESeq2
0
0
Entering edit mode
5.0 years ago

Hello everyone: I am new in using DEseq2. I am trying to analyze differential expression data using a publicly available data for my final project. It consists of 3 replications of control and three replications for treatment. I am trying to compare differential gene expression between treatment and control but i am getting error message which i do not know who to fix. Here is my code. Please help me fix the problem. Codes:

countdata<-read.csv("final project.csv",header = T,row.names=1)
groups<-c("control","control","control","treatment","treatment","treatment")
coldata<-data.frame(row.names=colnames(countdata),groups)
dds = DESeqDataSetFromMatrix(countData=countdata,
                             colData=colData,
                             design=~groups)
dds = DESeq(dds)
nrow(dds)

Error mesage:> dds = DESeqDataSetFromMatrix(countData=countdata,
+                              colData=colData,
+                              design=~groups)
Error in `rownames<-`(`*tmp*`, value = colnames(countData)) : 
  attempt to set 'rownames' on an object with no dimensions

Thank you in advance.

rna-seq • 1.4k views
ADD COMMENT
2
Entering edit mode

Come on. Give people something they can work with. Tell us what head countdata looks like, at least.

ADD REPLY

Login before adding your answer.

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