Load the results of the global run

In [1]:
from basepair.modisco.results import ModiscoResult
from basepair.imports import *

from basepair.config import get_data_dir
ddir = get_data_dir()

mdir = Path(f"{ddir}/processed//chipnexus/exp/models/oct-sox-nanog-klf/models/n_dil_layers=9/modisco/valid/merged")

/users/avsec/bin/anaconda3/envs/chipnexus/lib/python3.6/site-packages/h5py/__init__.py:36: FutureWarning: Conversion of the second argument of issubdtype from `float` to `np.floating` is deprecated. In future, it will be treated as `np.float64 == np.dtype(float).type`.
  from ._conv import register_converters as _register_converters
Using TensorFlow backend.
In [2]:
mr = ModiscoResult(mdir / "modisco.h5")
mr.open()
In [3]:
mr.metaclusters()
Out[3]:
['metacluster_0', 'metacluster_1']
In [4]:
mr.plot_pattern("metacluster_0", "pattern_1", trim_frac=0.08, kind=['sequence', 'merged_contrib_scores']);
In [5]:
mr.plot_all_patterns(trim_frac=0.08, kind=['sequence', 'merged_contrib_scores'])

# seqlets assigned to patterns: 6532 / 224651 (3%)