Hi All,

I am getting an error msgs while running DESeq2 Bioconductor in R studio. When I tried without using the gene_ID column and it does work but then I am unable to see which genes are DEGs.

My Input file:

gene_ID t1.control t2.control t3.control t4.sample t5.sample t6.sample

ENSCAFG1 363 358 409 549 243 437

ENSCAFG2 294 246 242 429 211 300

ENSCAFG3 60 96 26 90 16 17

ENSCAFG4 93 81 89 126 64 79

ENSCAFG5 208 108 174 259 159 177

ENSCAFG6 467 394 412 570 311 373

My R code is:

library(DESeq2)

setwd("/file/path/")

countdata=read.table("file5", header = TRUE)

condition <- factor(c(rep("control",3), rep("sample",3)))

condition

timepoints <- factor(c(rep("t1",1), rep("t2",1), rep("t3",1), rep("t4",1), rep("t5",1), rep("t6",1)))

timepoints

sampleTable <- data.frame(condition = as.factor(condition),
                          timepoints = as.factor(timepoints))

sampleTable

rownames(sampleTable) <- colnames(countdata)

#### Error in `.rowNamesDF<-`(x, value = value) : invalid 'row.names' length

deseq <- DESeqDataSetFromMatrix(countData = countdata,
                                    colData = sampleTable,
                                    design = ~condition+timepoints)

deseq
deseq <- deseq[ rowSums(counts(deseq)) > 1, ] ### remove rows that have only 0 and 1
d.deseq <- DESeq(deseq)                    
d.deseq

####  Error in DESeqDataSetFromMatrix(countData = countdata, colData = sampleTable,  : 
  ncol(countData) == nrow(colData) is not TRUE

Much appreciate for any help.