Mercurial > repos > rliterman > csp2
diff CSP2/CSP2_env/env-d9b9114564458d9d-741b3de822f2aaca6c6caa4325c4afce/lib/python3.8/site-packages/pybedtools/stats.py @ 68:5028fdace37b
planemo upload commit 2e9511a184a1ca667c7be0c6321a36dc4e3d116d
author | jpayne |
---|---|
date | Tue, 18 Mar 2025 16:23:26 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/CSP2/CSP2_env/env-d9b9114564458d9d-741b3de822f2aaca6c6caa4325c4afce/lib/python3.8/site-packages/pybedtools/stats.py Tue Mar 18 16:23:26 2025 -0400 @@ -0,0 +1,26 @@ +import multiprocessing +from . import helpers +import pybedtools + + +def random_jaccard(x, y, genome_fn, shuffle_kwargs, jaccard_kwargs): + z = x.shuffle(g=genome_fn, **shuffle_kwargs).sort() + result = z.jaccard(y, **jaccard_kwargs) + helpers.close_or_delete(z) + return result + + +def random_intersection(x, y, genome_fn, shuffle_kwargs, intersect_kwargs): + z = x.shuffle(g=genome_fn, **shuffle_kwargs) + zz = z.intersect(y, stream=True, **intersect_kwargs) + result = len(zz) + helpers.close_or_delete(z, zz) + return result + + +def random_intersection_bp(x, y, genome_fn, shuffle_kwargs, intersect_kwargs): + z = x.shuffle(g=genome_fn, **shuffle_kwargs) + zz = z.intersect(y, stream=True, **intersect_kwargs) + result = sum(len(i) for i in zz) + helpers.close_or_delete(z, zz) + return result