DiffBind: Model must include count data for contrasts. I am trying to use DiffBind for my ATACseq experiment going through the vignette for DiffBind. I was able to import the data and make a heatmap with plot(data). But somehow calling data &lt;- dba.contrast(data, categories=DBA_TISSUE) function is not working for me giving a warning message:


  Warning message:
  Model must include count data for contrasts. 


I tried running beforehand: data &lt;- dba.count(data) but it is giving an error:


  Error: Error processing one or more read files. Check warnings().
  In addition: Warning messages:
  1:
  2:
  3:
  4:
  5:
  6:
  7:  


My sample_sheet.csvlooks like that:



I just need to compare Dura vs Microglia which is Tissue column. How could this error be fixed? Any suggestions would be greatly appreciated.
The problem was that I was using wrong files. I ran a pipeline from here and it generates several folders with many files. The .bam files that should be used are located in the aligned_mm10 folder and are with the extension: .pe.q10.sort.rmdup.bam.