In [1]:
import numpy as np
motif_name = 'E2F6'
aitacdir = "/mnt/lab_data2/msharmin/oc-atlas/DanSkinData/fold_a_alex/{}".format(motif_name)

ave_filt_infl = np.load(aitacdir+"/ave_filt_infl.npy")
In [8]:
from matlas.matches import DenovoAitac
motif_name = 'E2F6'
aitacdir = "/mnt/lab_data2/msharmin/oc-atlas/DanSkinData/fold_a_alex/{}".format(motif_name)

ob = DenovoAitac(aitacdir, influence=ave_filt_infl)
# ob.fetch_tomtom_matches(
#             meme_db="/mnt/lab_data/kundaje/users/msharmin/annotations/HOCOMOCOv11_core_pwms_HUMAN_mono.renamed.nonredundant.annotated.meme",
#             database_name="HOCOMOCO.nonredundant.annotated",
#             save_report=True, tomtom_dir= "{0}/{1}_tomtomout".format(aitacdir, "HOCOMOCO.nonredundant.annotated"))
ob.load_matched_motifs(database_name="HOCOMOCO.nonredundant.annotated")
ob.get_motif_per_celltype(match_threshold=0.05, database_name="HOCOMOCO.nonredundant.annotated")
pattern_tab, pattern_dict = ob.visualize_pattern_table()
tf_tab, tf_dict = ob.visualize_tf_table("Aitac")
In [9]:
from vdom.helpers import (b, summary, details)
from IPython.display import display

display(details(summary('Click here for ', b('Denovo Patterns'), ' by ', b('{}'.format('Aitac')),
                        ' in ', b(motif_name),
                        ": #{}".format(len(pattern_dict)),
                       ), pattern_tab))
Click here for Denovo Patterns by Aitac in E2F6: #64
Pattern NameTF Name(s)AitacInfluence
filter456.919789252985563e-33
filter406.9099284916703e-33
filter316.887741778710959e-33
filter626.884043993217736e-33
filter176.877881017395697e-33
filter36.874183231902473e-33
filter216.866787660916027e-33
filter276.855694304436356e-33
filter196.854461709271948e-33
filter186.851996518943132e-33
filter466.849531328614316e-33
filter76.832274996312607e-33
filter366.822414234997344e-33
filter126.819949044668529e-33
filter16.79652973654478e-33
filter166.792831951051557e-33
filter296.787901570393925e-33
filter56.787901570393925e-33
filter246.787901570393925e-33
filter46.787901570393925e-33
filter116.787901570393925e-33
filter106.787901570393925e-33
filter386.787901570393925e-33
filter06.787901570393925e-33
filter426.787901570393925e-33
filter566.787901570393925e-33
filter556.787901570393925e-33
filter546.787901570393925e-33
filter576.787901570393925e-33
filter536.787901570393925e-33
filter15HCLUST-100_MYB.UNK.0.A6.785436380065109e-33
filter486.785436380065109e-33
filter96.782971189736294e-33
filter496.782971189736294e-33
filter206.778040809078663e-33
filter476.776808213914254e-33
filter256.765714857434585e-33
filter446.759551881612545e-33
filter136.755854096119322e-33
filter266.754621500954914e-33
filter236.745993334804059e-33
filter306.739830358982019e-33
filter286.732434787995573e-33
filter436.732434787995573e-33
filter86.72873700250235e-33
filter356.725039217009125e-33
filter416.718876241187087e-33
filter396.711480670200639e-33
filter526.710248075036232e-33
filter326.710248075036232e-33
filter6HCLUST-104_ARID5B.UNK.0.A6.707782884707416e-33
filter376.701619908885376e-33
filter346.699154718556561e-33
filter606.69422433789893e-33
filter226.685596171748075e-33
filter616.68313098141926e-33
filter636.675735410432812e-33
filter146.65971167329551e-33
filter26.653548697473471e-33
filter336.623966413527683e-33
filter506.60301229573275e-33
filter586.577127797280185e-33
filter516.41565783074276e-33
filter596.363888833837631e-33
In [10]:
display(details(summary('Click here for ', b('Motifs'), ' by ', b('{}'.format('Aitac')),
                        ' in ', b(motif_name),
                        ": #{}".format(len(tf_dict)),
                       ), tf_tab))
Click here for Motifs by Aitac in E2F6: #2
TF NamePattern(s)
HCLUST-100_MYB.UNK.0.A
Pattern NameAitacSignificance
filter150.0317837
HCLUST-104_ARID5B.UNK.0.A
Pattern NameAitacSignificance
filter60.033925199999999996
In [ ]: