I'm using DiffBind for ATAC-seq data, and everything runs smoothly except one thing -> I cant make a dba.report for one of my contrasts. All steps from loading -> counting -> analyzing seems to work, I get differently called peaks etc. I have 4 different conditions and 6 contrasts (output after dba.analyze):
6 Contrasts: Group1 Members1 Group2 Members2 DB.DESeq2 1 Olig2_Non_injured 2 Olig2_Injured 2 1 2 Olig2_Non_injured 2 WT_Non_injured 2 0 3 Olig2_Non_injured 2 WT_Injured_24h 2 50 4 Olig2_Injured 2 WT_Non_injured 2 2509 5 Olig2_Injured 2 WT_Injured_24h 2 8 6 WT_Non_injured 2 WT_Injured_24h 2 952
What really puzzles me is that I have done the exact same workflow for all my samples, but it's just one single contrast (contrast = 1) that produce this error:
Error in `colnames<-`(`*tmp*`, value = c("Chr", "Start", "End", "Conc", : 'names' attribute  must be the same length as the vector 
Furthermore, the same samples are (as you see in my contrast table) included in other comparisons, which works fine, hence, it's only for contrast=1 this happens.
Do anyone have any idea what could be the reason for this?