Ejemplo n.º 1
0
import numpy as np
from scipy.integrate import odeint 

from DS_solver import bioDS
from cell_cycle import cell
test=cell(path='SinIR_delay/')
test.load()
test.set("Spo0A_p",0)
test.simulate(tmax=30)
test.update()
test.clear()
for i in range(200):
	test.simulate(tmin=i*0.5,tmax=i*0.5+0.5,tsteps=50)
	print(1)
	test.parameterTurn(sigma=0.25)



test.display(['Ptap','Pmot'],save=True,fname='test')
#v_line=[2,6,8,12,14,18,20,24,26,30,32,36,38]
Ejemplo n.º 2
0
import numpy as np
from scipy.integrate import odeint 
import pickle
from DS_solver import bioDS
from cell_cycle import cell
import matplotlib.pyplot as plt
for ppath in ['SinIR_delayf2_5/','SinIR_delayf2_4/','SinIR_delayf2_d/']:
	print(ppath)	
	test=cell(path=ppath)
	test.load()
	test.set("Spo0A_p",0)
	test.simulate(tmax=30)
	test.update()
	test.clear()
	for i in range(20000):
		test.simulate(tmin=i*0.5,tmax=i*0.5+0.5,tsteps=30)
		test.parameterTurn(sigma=0.25)
		if (i/1000)*1000==i:
			print i
	tt=test.trajectory[0]
	yy=test.getTrajectory('Ptap')
	zz=test.getTrajectory('Pmot')
	pickle.dump([tt,yy,zz],open(ppath+"f.pickle",'wb'))
	for threshold in [60,80,100,120,140,160,180,200,250,300,400,500]:
		chain=[]
		mot=[]
		s=0
		t=0
		if yy[0]>threshold:
			s=1
		for i in range(1,len(tt)):
Ejemplo n.º 3
0
import numpy as np
from scipy.integrate import odeint 
import pickle
from DS_solver import bioDS
from cell_cycle import cell
import matplotlib.pyplot as plt
test=cell(path='SinIR_feedback/')
test.load()
test.set("Spo0A_p",0)
test.simulate(tmax=30)
test.update()
test.clear()
for i in range(20000):
	test.simulate(tmin=i*0.5,tmax=i*0.5+0.5,tsteps=30)
	test.parameterTurn(sigma=0.25)
	if (i/100)*100==i:
		print i
tt=test.trajectory[0]
yy=test.getTrajectory('Ptap')
pickle.dump([tt,yy],open("SinIR_feedback/ttyy-fb_2.pickle",'wb'))
for threshold in [60,80,100,120,140,160,180,200,250,300,400,500]:
	chain=[]
	mot=[]
	s=0
	t=0
	if yy[0]>threshold:
		s=1
	for i in range(1,len(tt)):
		if s==0:
			if yy[i]>threshold:
				s=1
Ejemplo n.º 4
0
import numpy as np
from scipy.integrate import odeint 

from DS_solver import bioDS
from cell_cycle import cell
test=cell(path='SinIR_delayf2_d/')
test.load()
test.set("Spo0A_p",0)
test.simulate(tmax=30)
test.update()
test.clear()
for i in range(200):
	test.simulate(tmin=i*0.5,tmax=i*0.5+0.5,tsteps=50)
	print(1)
	test.parameterTurn(sigma=0.25)



test.display(['Ptap','Pmot'],save=True,fname='test')
#v_line=[2,6,8,12,14,18,20,24,26,30,32,36,38]