Hi I had a question about the following example in the diffbind tutorial:
Plot effect of a range of filter values and then apply filter
data(tamoxifen_counts) rate.max <- dba.count(tamoxifen, peaks=NULL, filter=0:250) rate.sum <- dba.count(tamoxifen, peaks=NULL, filter=0:250,filterFun=sum) plot(0:250,rate.max/rate.max,type='l',xlab="Filter Value",ylab="Proportion Retained Sites") lines(0:250,rate.sum/rate.sum,col=2) tamoxifen <- dba.count(tamoxifen,peaks=NULL,filter=125,filterFun=sum) tamoxifen
I understand everything done here except for why a filter of 125 was applied as the final step, how is this filter value obtained as the optimal from the proportion of retained sites? Thank you.