Trying to generate replicate masks for the dba.count function from a .csv where I list 3 replicates under DBA_CONDITION. Want to do dba.counts with minMembers = 3 to include only those that are in all 3 of each replicate group for the binding matrix. Any clues on how to do this? Thanks.
I'm not sure exactly what your sample sheet looks like for this. Do you have three replicates of each condition, where each replicated sample has the same
DBA_CONDITION value, but a different replicate number? And you want to keep peaks for one specific value condition?
Is this the same project as in question #DiffBind contrasts issue And you would to, for example, use
dba.count() to make a binding matrix with all the peaks that are in all three MALE/TREATMENT replicates? If that is what you want to do, this should work:
> masks <- CHIP_adult_BN_peaks$masks > CHIP_adult_BN_mask_MaleTRE <- masks$MALE & masks$TREATMENT > CHIP_adult_BN_peaks_MaleTRE <- dba(CHIP_adult_BN_peaks,mask=CHIP_adult_BN_mask_MaleTRE, minMembers=3) > CHIP_adult_BN_consensus_MaleTRE <- dba(CHIP_adult_BN_peaks_MaleTRE, bRetrieve=TRUE) > CHIP_adult_BN_counts <- dba.count(CHIP_adult_BN_peaks, peaks=CHIP_adult_BN_consensus_MaleTRE)