I started trying to use Homer for finding overlaps of my 20 sample data. It finds the overlaps pretty good. However, I cannot get that taking the mean of the overlaps option work. Is there a way to get it worked. I found this github page https://github.com/shraddhapai/mergePeaks_resolve . But that boost compiler gave me lots of errors.
Your method is nice. But I have already got overlapping regions. Homer finds those overlapping regions but I does not calculate the mean of the peak heights. I think your line of code, also finds the overlapping peaks which might generate problems for me?
I meant to use the original bed files, and not the one generated by Homer. You start directly with your 20 or so files by concatenating them in order to merge them. This method will "merge peaks" (union) and does not intersect them. Your question is confusing: Do you want to merge your peaks or to get the overlaps ? These are two different things.
I only wanna get overlaps. But what I am saying is I have already got the overlaps with various programs such as Homer, Bedops/Bedmap. Therefore, I asked my question specifically about is there a way to invoke average option of the Homer.
What I need for my project I have 13 tumour and 7 normal samples. I want to get overlapping regions within these data set. For normal samples, if the overlap region no is 7 then I call it high confidence peaks. where I call less than 7 as low peaks. same for the Tumour but 13 for high confidence and less than 13 for low confidence.
Your method is nice. But I have already got overlapping regions. Homer finds those overlapping regions but I does not calculate the mean of the peak heights. I think your line of code, also finds the overlapping peaks which might generate problems for me?
I meant to use the original bed files, and not the one generated by Homer. You start directly with your 20 or so files by concatenating them in order to merge them. This method will "merge peaks" (union) and does not intersect them. Your question is confusing: Do you want to merge your peaks or to get the overlaps ? These are two different things.
I only wanna get overlaps. But what I am saying is I have already got the overlaps with various programs such as Homer, Bedops/Bedmap. Therefore, I asked my question specifically about is there a way to invoke average option of the Homer.
What I need for my project I have 13 tumour and 7 normal samples. I want to get overlapping regions within these data set. For normal samples, if the overlap region no is 7 then I call it high confidence peaks. where I call less than 7 as low peaks. same for the Tumour but 13 for high confidence and less than 13 for low confidence.