Square heatmaps#
CLI mode#
conf-hicmap_5.yml
trackc:
- ax: t1
height: 3
track_type: hicmap
track_para:
mat:
method: extractCisContact
clr: ./github/seqyuan/trackc_data/examples/HPDE6C7.chr18.mcool::/resolutions/25000
mapC:
map_type: squ
label: upper
- ax: t2
height: 3
hspace: 0.1
track_type: hicmap
track_para:
mat2:
method: extractCisContact
clr: ./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000
mapC:
map_type: squ
label: lower
- ax: t3
height: 3
hspace: 0.1
track_type: hicmap
track_para:
mat:
method: extractCisContact
clr: ./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000
mapC:
map_type: squ
symmetric: True
- ax: t4
height: 3
hspace: 0.1
track_type: hicmap
track_para:
mat:
method: extractCisContact
clr: ./github/seqyuan/trackc_data/examples/HPDE6C7.chr18.mcool::/resolutions/25000
mat2:
method: extractCisContact
clr: ./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000
mapC:
cmap:
- PuBu
- Purples
map_type: squ
label:
- Normal
- Tumor
trackc cli conf-hicmap_1.yml -s 3,1 -r 18:64000000-71000000
API mode#
[10]:
import trackc as tc
ten = tc.tenon(figsize=(4, 1))
ten.add(pos="bottom", height=4)
region = "18:64000000-71000000"
clr1 = "./github/seqyuan/trackc_data/examples/HPDE6C7.chr18.mcool::/resolutions/25000"
clr2 = "./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000"
mat = tc.tl.extractCisContact(clr=clr1, region=region)
mat2 = tc.tl.extractCisContact(clr=clr2, region=region)
tc.pl.mapC(
ax=ten.axs(0),
mat=mat,
mat2=mat2,
cmap=[tc.pa.cmaps.fruitpunch, "PuBu"],
map_type="squ",
label=["Normal", "Tumor"],
)
no max min range
maxrange: 160.0 minrange: 1.0
no max min range
maxrange: 174.0 minrange: 1.0