what could be the source of this erro
0
0
Entering edit mode
Abdul • 0
@9ec788b1
Last seen 15 months ago
Sweden

Enter the body of text here

Error in DESeqDataSetFromMatrix(countData = as.matrix(count_Table1), colData = sample_table, : ncol(countData) == nrow(colData) is not TRUE

timepoint<-factor(rep(c("8h","24h","48h","72h","96h"),each=37),levels = c("8h","24h","48h","72h","96h"))
sample_table<-data.frame("samples"=colnames(count_Table1),"timepoint"=timepoint,"user.ID"= a)
library(edgeR)
meanlog2CPM<-rowMeans(log2(cpm(count_Table1) + 1))
hist(meanlog2CPM)
dim(count_Table1)
count_Table1<-count_Table1[meanlog2CPM >1,]
dim(count_Table1)
library(DESeq2)
dds<-DESeqDataSetFromMatrix(countData=count_Table1,colData=sample_table,design = ~timepoint)
normcount<-rlog(dds,blind =FALSE)
DrugVsDiseasedata • 715 views
ADD COMMENT
0
Entering edit mode

It looks like the explanation is right there. Can you explain and show us how you are sure the number of columns in the count data is the same number of rows in the coldata?

Also, why are you putting log CPM into DESeq2? You must know that's not valid.

ADD REPLY
0
Entering edit mode

Thank you for the your suggestions and reply. Yes the number of rows and columns are different and i was wondring how to make them same since I'm new in R . Thank you

ADD REPLY
0
Entering edit mode

Again, why are you putting logCPM into DESeq? If you put the wrong input in, your output will be garbage. If you want to know why your row and column names don't match, you are going to have to look at them.

ADD REPLY

Login before adding your answer.

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