Quantification of pseudo-bulked signal at cluster specific peaks from snap.tl.merge_peaks function

[skpeterson]

Is anyone aware of a way to quantify signal in the merged peak set that is obtained using the function snap.tl.merge_peaks? The output of the merge_peaks fn is a boolean dataframe. I was hoping for a pseudo-bulk signal based on cluster on the merged peak set, ie if instead of boolean 'true'/'false' it would be a signal value. The desired output would have the column names 'Peaks' 'Cluster0' 'Cluster1' and the rows would be the genomic regions 'chr#:xx-xx' and then the signal_value would be shown in the dataframe.