Search
0
13 months ago by
rbronste60
rbronste60 wrote:

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.

modified 13 months ago by Rory Stark2.5k • written 13 months ago by rbronste60
1
13 months ago by
Rory Stark2.5k
CRUK, Cambridge, UK
Rory Stark2.5k wrote:

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
minMembers=3)
peaks=CHIP_adult_BN_consensus_MaleTRE)