I am trying the following code and getting an error I am unsure how to fix:
countData<-read.csv("p14_binding.csv") rownames(countData) <- countData[,1] countData[,1] <- NULL sampleNames<-c("MB1", "MB2", "MB3", "FB1", "FB2", "FB3", "MM1", "MM2", "MM3", "FM1", "FM2", "FM3") sampleSex<-c("male", "male", "male", "female", "female", "female", "male", "male", "male", "female", "female", "female") sampleAge<-c("p14", "p14", "p14", "p14", "p14", "p14", "p14", "p14", "p14", "p14", "p14", "p14") colData<-data.frame(sampleName=sampleNames, sex=sampleSex, age=sampleAge) dds<-DESeqDataSetFromMatrix(countData = countData, colData = colData, design = ~ sex + age + age:sex)
Receiving the following error, I think having something to do with the way my input matrix is designed - as it was imported from diffbind so is not a list of genes but rather a list of binding counts from which I cut the chromosome start-end sites:
Error in DESeqDataSet(se, design = design, ignoreRank) :
design contains one or more variables with all samples having the same value,
remove these variables from the design
Thanks for your help!
LANGUAGE_DETECTION = []