#!/usr/bin/env python
# coding: utf-8

# In[1]:


import stablemotifs
import biolqm
from colomoto_jupyter import tabulate


# In[2]:


model = biolqm.load("MyeloidRules.txt", "booleannet")


# In[3]:


get_ipython().run_line_magic('time', 'stm = stablemotifs.load(model)')


# In[4]:


tabulate(stm.attractors)


# In[5]:


stm.stable_motifs


# In[6]:


s = stm.reprogramming_to_attractor(EgrNab=1, CEBPA=1)
s


# In[7]:


s.as_graph()


# In[8]:


s.as_table()


# In[ ]: