Entering edit mode
Hi, I was having an issue when using Glimma package to generate plots. My data processing before Glimma was in DESeq2.
dds = DESeq(dds)
glimmaMA(dds)
I tried another plot glimmaMDS(dds), it was operating normally. However, whenever I try the glimmaMA or glimmaVolcano it doesnt work, and the error message is the same:
302 of 12813 genes were filtered out in DESeq2 tests Error in buildXYData(table, status, main, display.columns, anno, counts, : Length of groups must be equal to the number of columns in counts.
The R feedback was
- stop("Length of groups must be equal to the number of columns in counts.\n")
- buildXYData(table, status, main, display.columns, anno, counts, xlab, ylab, status.cols, sample.cols, groups, transform.counts)
- glimmaVolcano.DESeqDataSet(dds)
- glimmaVolcano(dds)
Could you please help with potential cause and solutions regarding this error? Thank you very much!!