Hello,

I am new to the field of bioinformatics and currently using your R package (DiffBind_2.6.6) for the analysis of my ATAC-seq data to identify differences in chromatin accessibility between groups. Due to the nature of the project I want to construct a specific consensus set to use for dba.count, but I am struggling to achieve this. Your help/advice would be much appreciated.

In the consensus peakset I would like to include:

       - the peaks that are present in at least 6 out of 7 control samples (regardless of the experimental group)

       - the peaks in at least 4 out of 8 experimental samples (regardless of the control group)

Is there an easy way to construct a consensus peakset meeting the requirements stated above?

As I couldn't figure out how to construct a consensus peakset as such, I constructed two separate consensus peaksets (one for each of the requirement using the code below), but I have not yet succeeded in merging these to use them in dba.count.

control_consensus <- dba.peakset(control, minOverlap = 6)
consensus_peaks_control <- dba.peakset(control_consensus, control_consensus$masks$Consensus, bRetrieve = TRUE)

exp_consensus <- dba.peakset(exp, minOverlap = 4)
consensus_peaks_exp <- dba.peakset(exp_consensus, exp_consensus$masks$Consensus, bRetrieve = TRUE)

If a more direct method is not available, how would you advice to combine these two separate peaksets for further use?

Thank you very much,

Best,

Rob