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 = []
