@@ -39,11 +39,11 @@ dna_rs = @reaction_network begin
39
39
k6, DNAR --> DNA + P
40
40
end
41
41
rates = [:k1 => 0.5 ,
42
- :k2 => (20 * log (2.0 ) / 120.0 ),
43
- :k3 => (log (2.0 ) / 120.0 ),
44
- :k4 => (log (2.0 ) / 600.0 ),
45
- :k5 => 0.025 ,
46
- :k6 => 1.0 ]
42
+ :k2 => (20 * log (2.0 ) / 120.0 ),
43
+ :k3 => (log (2.0 ) / 120.0 ),
44
+ :k4 => (log (2.0 ) / 600.0 ),
45
+ :k5 => 0.025 ,
46
+ :k6 => 1.0 ]
47
47
tf = 1000.0
48
48
u0 = [:DNA => 1 , :mRNA => 0 , :P => 0 , :DNAR => 0 ]
49
49
prob = DiscreteProblem (dna_rs, u0, (0.0 , tf), rates, eval_module = @__MODULE__ )
@@ -103,7 +103,7 @@ oscil_rs = @reaction_network begin
103
103
0.05 , SP2 --> 0
104
104
end
105
105
u0 = [:X => 200.0 , :Y => 60.0 , :Z => 120.0 , :R => 100.0 , :S => 50.0 , :SP => 50.0 ,
106
- :SP2 => 50.0 ] # Hill equations force use of floats!
106
+ :SP2 => 50.0 ] # Hill equations force use of floats!
107
107
tf = 4000.0
108
108
prob = DiscreteProblem (oscil_rs, u0, (0.0 , tf), eval_module = @__MODULE__ )
109
109
"""
@@ -121,7 +121,7 @@ specs_sym_to_name = Dict(:S1 => "R(a,l)",
121
121
:S8 => " A(Y~P,r!1).R(a!1,l)" ,
122
122
:S9 => " A(Y~P,r)" )
123
123
rsi = Dict (:R0 => 1 , :L0 => 2 , :A0 => 3 , :kon => 4 , :koff => 5 ,
124
- :kAon => 6 , :kAoff => 7 , :kAp => 8 , :kAdp => 9 )
124
+ :kAon => 6 , :kAoff => 7 , :kAp => 8 , :kAdp => 9 )
125
125
params = (5360 , 1160 , 5360 , 0.01 , 0.1 , 0.01 , 0.1 , 0.01 , 0.1 )
126
126
rs = @reaction_network begin
127
127
kon, S1 + S2 --> S4
146
146
rates = [:kon , :kAon , :koff , :kAoff , :kAp , :kAdp ] .=>
147
147
params[[rsi[:kon ], rsi[:kAon ], rsi[:koff ], rsi[:kAoff ], rsi[:kAp ], rsi[:kAdp ]]]
148
148
u0 = [:S1 => params[1 ], :S2 => params[2 ], :S3 => params[3 ], :S4 => 0 , :S5 => 0 ,
149
- :S6 => 0 , :S7 => 0 , :S8 => 0 , :S9 => 0 ]
149
+ :S6 => 0 , :S7 => 0 , :S8 => 0 , :S9 => 0 ]
150
150
tf = 100.0
151
151
prob = DiscreteProblem (rs, u0, (0.0 , tf), rates, eval_module = @__MODULE__ )
152
152
"""
@@ -220,7 +220,7 @@ rn = @reaction_network begin
220
220
c8, P2G --> P2 + G
221
221
end
222
222
rnpar = [:c1 => 0.09 , :c2 => 0.05 , :c3 => 0.001 , :c4 => 0.0009 , :c5 => 0.00001 ,
223
- :c6 => 0.0005 , :c7 => 0.005 , :c8 => 0.9 ]
223
+ :c6 => 0.0005 , :c7 => 0.005 , :c8 => 0.9 ]
224
224
varlabels = [" G" , " M" , " P" , " P2" , " P2G" ]
225
225
u0 = [:G => 1000 , :M => 0 , :P => 0 , :P2 => 0 , :P2G => 0 ]
226
226
tf = 4000.0
0 commit comments