This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Differential expression analysis using DESeq2

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

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

0 answers

No answers yet.

Log in to answer this question.