from basepair.imports import *
modisco_dir = "/users/amr1/bpnet/basepair/src/cage/modisco/"
mr = ModiscoResult(modisco_dir+"modisco.h5")
mr.open()
len(mr.patterns())
grad_modisco_dir = "/users/amr1/bpnet/basepair/src/cage/grad_modisco/"
grad_mr = ModiscoResult(modisco_dir+"modisco.h5")
grad_mr.open()
len(grad_mr.patterns())
from modisco.visualization import viz_sequence
import matplotlib.pyplot as plt
%matplotlib inline
print("DeepLIFT Modisco Patterns")
for pattern_name in mr.patterns():
viz_sequence.plot_weights(mr.get_pattern(pattern_name).seq, figsize=(20,3))
print("gradxinp Modisco Patterns")
for pattern_name in grad_mr.patterns():
viz_sequence.plot_weights(grad_mr.get_pattern(pattern_name).seq, figsize=(20,3))