added interval functionality to LinkPeaks #561
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Hi Tim,
I've got another PR here for some added functionality to the LinkPeaks methods. We're in a situation where we'd like to look at potential very long-range enhancers (> 500kb and < 1mb) across the whole genome, but we've already called linked peaks using the default 500 kb distance. We'd like to avoid re-running the peak-gene linkage on all the peaks within 500 kb of the TSS and only search in our defined interval. This commit allows users to specify a minimum distance for the peak from the TSS along with the already included maximum distance from TSS. I think people will potentially find it useful for enhancer searches by allowing them to exclude promoter regions, or to search a defined interval for a specific gene, along with other use cases.