Question: Differential expression analysis using DESeq2
0
gravatar for prameelaawale
11 weeks ago by
prameelaawale0 wrote:

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 • 177 views
ADD COMMENTlink modified 11 weeks ago by genomax69k • written 11 weeks ago by prameelaawale0
2

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

ADD REPLYlink written 11 weeks ago by swbarnes26.0k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1109 users visited in the last hour