# is monostable at each FP layer.
import DSGRN
import Data_added_GetParameterIndex as DGetParaIndex
import MorseList_to_ParaList
import Pitogether
getparameterindex_object = DGetParaIndex.Database.GetParameterIndex
database = DSGRN.Query.Database.Database("EMT_Zeb1notE.db") ###****** Choose
# database file
# properly ******###
#Get the list of ParameterIndices with monostable.
monostable_query_object = DSGRN.Query.MonostableQuery.MonostableQuery(database)
monostable_pilists = MorseList_to_ParaList.MorseList_to_ParaList(
monostable_query_object.matches(), getparameterindex_object, database)
monostable_pitogether = Pitogether.pitogether(monostable_pilists)
#Get the list of ParameterIndices with each point layer.
import PointLayers
Pl_pitogether = PointLayers.Pl_pitogether
#Get the list of ParameterIndices of each hexcode layer.
import Zeb1HexcodeLayer
hexlayers = Zeb1HexcodeLayer.zeb1_layers ###****** Choose this according
# to the databese file being used ******###
#Get the number of ParameterIndices with each point layer in each
# hexcode layer.
Plmono_in_hexl = []
with open("players.txt", "a") as f:
import DSGRN
import Data_added_GetParameterIndex as DGetParaIndex
import MorseList_to_ParaList
import Pitogether
getparameterindex_object = DGetParaIndex.Database.GetParameterIndex
database = DSGRN.Query.Database.Database("EMT_Ovol2notE.db") ###****** Change
# this database
# file
# as needed ******###
bounds002 = {"Zeb1": 0, "Snail1": 0, "Ovol2": 2}
matches002 = DSGRN.Query.SingleFixedPointQuery.SingleFixedPointQuery(
database, bounds002).matches()
pilists_002 = MorseList_to_ParaList.MorseList_to_ParaList(
matches002, getparameterindex_object, database)
P002_pitogether = Pitogether.pitogether(pilists_002)
bounds330 = {"Zeb1": 3, "Snail1": 3, "Ovol2": 0}
matches330 = DSGRN.Query.SingleFixedPointQuery.SingleFixedPointQuery(
database, bounds330).matches()
pilists_330 = MorseList_to_ParaList.MorseList_to_ParaList(
matches330, getparameterindex_object, database)
P330_pitogether = Pitogether.pitogether(pilists_330)
import Nstable_in_Hexlayers
print('The number of parameter nodes in each hexcode layer for each nstable '
'cases:')
nstable_pitogether = Nstable_in_Hexlayers.nstable_pitogether
hexlayers = Nstable_in_Hexlayers.hexlayers
#In each nstable and each hexcode layer, number of parameter nodes with E
import DSGRN
import Data_added_GetParameterIndex as DGetParaIndex
import MorseList_to_ParaList
import Pitogether
getparameterindex_object = DGetParaIndex.Database.GetParameterIndex
database = DSGRN.Query.Database.Database("EMT_Zeb1notE.db") ###****** Change
# this database file
# as needed ******###
#Get the list of ParameterIndices with monostable.
monostable_query_object = DSGRN.Query.MonostableQuery.MonostableQuery(database)
monostable_pilists = MorseList_to_ParaList.MorseList_to_ParaList(
monostable_query_object.matches(),getparameterindex_object,database)
monostable_pitogether = Pitogether.pitogether(monostable_pilists)
#Get the list of ParameterIndices with bistable.
bistable_query_object = DSGRN.Query.BistableQuery.BistableQuery(database)
bistable_pilists = MorseList_to_ParaList.MorseList_to_ParaList(
bistable_query_object.matches(),getparameterindex_object,database)
bistable_pitogether = Pitogether.pitogether(bistable_pilists)
nstable_query = DSGRN.Query.NstableQuery.NstableQuery
tristable_query_object = nstable_query(database, 3)
tri_matches_pre = tristable_query_object.matches()
fourstable_query_object = nstable_query(database, 4)
import DSGRN
import Data_added_GetParameterIndex as DGetParaIndex
import MorseList_to_ParaList
import Pitogether
getparameterindex_object = DGetParaIndex.Database.GetParameterIndex
database = DSGRN.Query.Database.Database("EMT_Zeb1notE.db") ###****** Choose
# the database
# file
# properly ******###
#Get the list of ParameterIndices with monostable.
monostable_query_object = DSGRN.Query.MonostableQuery.MonostableQuery(database)
monostable_pilists = MorseList_to_ParaList.MorseList_to_ParaList(
monostable_query_object.matches(), getparameterindex_object, database)
monostable_pitogether = Pitogether.pitogether(monostable_pilists)
#Get the list of ParameterIndices with FP(0,0,*,*,*,2).
bounds002 = {"Zeb1": 0, "Snail1": 0, "Ovol2": 2}
matches002 = DSGRN.Query.SingleFixedPointQuery.SingleFixedPointQuery(
database, bounds002).matches()
print('Number of MorseGraphIndices with FP(0,0,*,*,*,2):')
print(len(matches002))
pilists_002 = MorseList_to_ParaList.MorseList_to_ParaList(
matches002, getparameterindex_object, database)
#Put the ParameterIndices together
P002_pitogether = Pitogether.pitogether(pilists_002)
print('Number of ParameterIndices with FP(0,0,*,*,*,2):')
print(len(P002_pitogether))
#Get the list of ParameterIndices with FP(3,3,*,*,*,0).
import DSGRN
import Data_added_GetParameterIndex as DGetParaIndex
import MorseList_to_ParaList
import Pitogether
getparameterindex_object = DGetParaIndex.Database.GetParameterIndex
database = DSGRN.Query.Database.Database("EMT_Ovol2notE.db") ### For different
# databases,
# just need to replace the database file here.
#Get the list of ParameterIndices with monostable.
monostable_query_object = DSGRN.Query.MonostableQuery.MonostableQuery(database)
monostable_pilists = MorseList_to_ParaList.MorseList_to_ParaList(
monostable_query_object.matches(), getparameterindex_object, database)
monostable_pitogether = Pitogether.pitogether(monostable_pilists)
#Get the list of ParameterIndices with each FP(i,j,*,*,*,k).
FP_pitogether = []
for i in range(0, 4):
for j in range(0, 4):
for k in range(0, 3):
bounds = {"Zeb1": i, "Snail1": j, "Ovol2": k}
matches = DSGRN.Query.SingleFixedPointQuery\
.SingleFixedPointQuery(database,
bounds).matches()
FPpilists = MorseList_to_ParaList.MorseList_to_ParaList(
matches, getparameterindex_object, database)
FP_pitogether.append(Pitogether.pitogether(FPpilists))
mono_at_FP = []
bounds322 = {"Zeb1":3,"Snail1":2,"Ovol2":2}
matches322 = DSGRN.Query.SingleFixedPointQuery.SingleFixedPointQuery(database, bounds322).matches()
bounds330 = {"Zeb1":3,"Snail1":3,"Ovol2":0}
matches330 = DSGRN.Query.SingleFixedPointQuery.SingleFixedPointQuery(database, bounds330).matches()
bounds331 = {"Zeb1":3,"Snail1":3,"Ovol2":1}
matches331 = DSGRN.Query.SingleFixedPointQuery.SingleFixedPointQuery(database, bounds331).matches()
bounds332 = {"Zeb1":3,"Snail1":3,"Ovol2":2}
matches332 = DSGRN.Query.SingleFixedPointQuery.SingleFixedPointQuery(database, bounds332).matches()
#Collect in layers.
matches_Pl1 = matches002
pilist_Pl1 = MorseList_to_ParaList.MorseList_to_ParaList(matches_Pl1,
getparameterindex_object,database)
Pl1_pis = Pitogether.pitogether(pilist_Pl1)
matches_Pl2 = set(list(matches012)+list(matches102)+list(matches001))
pilist_Pl2 = MorseList_to_ParaList.MorseList_to_ParaList(matches_Pl2,
getparameterindex_object,database)
Pl2_pis = Pitogether.pitogether(pilist_Pl2)
matches_Pl3 = set(list(matches202) + list(matches112) + list(matches022) +
list(matches101) + list(matches011) + list(matches000))
pilist_Pl3 = MorseList_to_ParaList.MorseList_to_ParaList(matches_Pl3,
getparameterindex_object,database)
Pl3_pis = Pitogether.pitogether(pilist_Pl3)
matches_Pl4 = set(list(matches302) + list(matches212) + list(matches122) +
list(matches032) + list(matches201) + list(matches111) +
list(matches021) + list(matches100) + list(matches010))