Hey everyone,
I'm trying to use the hg.divide function, and it works when I load different datasets with hg.remote(), but not with hg.cooler(). I've tried various approaches, but I haven't been able to get it to work. Could maybe someone help me?
import higlass as hg
t1 = hg.cooler('../dir/file1.mcool')
t2 = hg.cooler('../dir/file2.mcool')
t3 = hg.divide(t1, t2).opts(
colorRange=["blue", "white"],
valueScaleMin=0.1,
valueScaleMax=10,
)
domain = (2e7, 3e7)
v1 = hg.view(t1, width=4).domain(x=domain)
v2 = hg.view(t2, width=4).domain(x=domain)
v3 = hg.view(t3, width=4).domain(x=domain)
view_lock = hg.lock(v1, v2, v3)
(v1 | v3 | v2).locks(view_lock)
Best,
Henrik