@@ -15,9 +15,9 @@ export prob_jump_dnarepressor, prob_jump_constproduct, prob_jump_nonlinrxs,
1515 prob_jump_diffnetwork
1616
1717"""
18- General structure to hold JumpProblem info. Needed since
19- the JumpProblem constructor requires the algorithm, so we
20- don't create the JumpProblem here.
18+ General structure to hold JumpProblem info. Needed since
19+ the JumpProblem constructor requires the algorithm, so we
20+ don't create the JumpProblem here.
2121"""
2222struct JumpProblemNetwork
2323 network:: Any # Catalyst network
@@ -49,7 +49,7 @@ Nsims = 8000
4949expected_avg = 5.926553750000000e+02
5050prob_data = Dict (" num_sims_for_mean" => Nsims, " expected_mean" => expected_avg)
5151"""
52- DNA negative feedback autoregulatory model. Protein acts as repressor.
52+ DNA negative feedback autoregulatory model. Protein acts as repressor.
5353"""
5454prob_jump_dnarepressor = JumpProblemNetwork (dna_rs, rates, tf, u0, prob, prob_data)
5555
@@ -65,7 +65,7 @@ Nsims = 16000
6565expected_avg = t -> rates[1 ] / rates[2 ] .* (1.0 - exp .(- rates[2 ] * t))
6666prob_data = Dict (" num_sims_for_mean" => Nsims, " expected_mean_at_t" => expected_avg)
6767"""
68- Simple birth-death process with constant production and degradation.
68+ Simple birth-death process with constant production and degradation.
6969"""
7070prob_jump_constproduct = JumpProblemNetwork (bd_rs, rates, tf, u0, prob, prob_data)
7171
@@ -84,7 +84,7 @@ Nsims = 32000
8484expected_avg = 84.876015624999994
8585prob_data = Dict (" num_sims_for_mean" => Nsims, " expected_mean" => expected_avg)
8686"""
87- Example with a mix of nonlinear reactions, including third order
87+ Example with a mix of nonlinear reactions, including third order
8888"""
8989prob_jump_nonlinrxs = JumpProblemNetwork (nonlin_rs, rates, tf, u0, prob, prob_data)
9090
@@ -105,7 +105,7 @@ u0 = [:X => 200.0, :Y => 60.0, :Z => 120.0, :R => 100.0, :S => 50.0, :SP => 50.0
105105tf = 4000.0
106106prob = DiscreteProblem (oscil_rs, u0, (0.0 , tf), eval_module = @__MODULE__ )
107107"""
108- Oscillatory system, uses a mixture of jump types.
108+ Oscillatory system, uses a mixture of jump types.
109109"""
110110prob_jump_osc_mixed_jumptypes = JumpProblemNetwork (oscil_rs, nothing , tf, u0, prob, nothing )
111111
@@ -153,10 +153,10 @@ u0 = [:S1 => params[1], :S2 => params[2], :S3 => params[3], :S4 => 0, :S5 => 0,
153153tf = 100.0
154154prob = DiscreteProblem (rs, u0, (0.0 , tf), rates, eval_module = @__MODULE__ )
155155"""
156- Multistate model from Gupta and Mendes,
157- "An Overview of Network-Based and -Free Approaches for Stochastic Simulation of Biochemical Systems",
158- Computation 2018, 6, 9; doi:10.3390/computation6010009
159- Translated from supplementary data file: Models/Multi-state/fixed_multistate.xml
156+ Multistate model from Gupta and Mendes,
157+ "An Overview of Network-Based and -Free Approaches for Stochastic Simulation of Biochemical Systems",
158+ Computation 2018, 6, 9; doi:10.3390/computation6010009
159+ Translated from supplementary data file: Models/Multi-state/fixed_multistate.xml
160160"""
161161prob_jump_multistate = JumpProblemNetwork (rs, rates, tf, u0, prob,
162162 Dict (" specs_to_sym_name" => specs_sym_to_name,
@@ -206,9 +206,9 @@ tf = 2000.0
206206prob = DiscreteProblem (rs, u0, (0.0 , tf), eval_module = @__MODULE__ )
207207
208208"""
209- Twenty-gene model from McCollum et al,
210- "The sorting direct method for stochastic simulation of biochemical systems with varying reaction execution behavior"
211- Comp. Bio. and Chem., 30, pg. 39-49 (2006).
209+ Twenty-gene model from McCollum et al,
210+ "The sorting direct method for stochastic simulation of biochemical systems with varying reaction execution behavior"
211+ Comp. Bio. and Chem., 30, pg. 39-49 (2006).
212212"""
213213prob_jump_twentygenes = JumpProblemNetwork (rs, nothing , tf, u0, prob, nothing )
214214
@@ -229,10 +229,10 @@ u0 = [:G => 1000, :M => 0, :P => 0, :P2 => 0, :P2G => 0]
229229tf = 4000.0
230230prob = DiscreteProblem (rn, u0, (0.0 , tf), rnpar, eval_module = @__MODULE__ )
231231"""
232- Negative feedback autoregulatory gene expression model. Dimer is the repressor.
233- Taken from Marchetti, Priami and Thanh,
234- "Simulation Algorithms for Comptuational Systems Biology",
235- Springer (2017).
232+ Negative feedback autoregulatory gene expression model. Dimer is the repressor.
233+ Taken from Marchetti, Priami and Thanh,
234+ "Simulation Algorithms for Comptuational Systems Biology",
235+ Springer (2017).
236236"""
237237prob_jump_dnadimer_repressor = JumpProblemNetwork (rn, rnpar, tf, u0, prob,
238238 Dict (" specs_names" => varlabels))
@@ -257,10 +257,10 @@ function getDiffu0(diffnetwork, N)
257257end
258258tf = 10.0
259259"""
260- Continuous time random walk (i.e. diffusion approximation) example.
261- Here the network in the JumpProblemNetwork is a function that returns a
262- network given the number of lattice sites.
263- u0 is a similar function that returns the initial condition vector.
260+ Continuous time random walk (i.e. diffusion approximation) example.
261+ Here the network in the JumpProblemNetwork is a function that returns a
262+ network given the number of lattice sites.
263+ u0 is a similar function that returns the initial condition vector.
264264"""
265265prob_jump_diffnetwork = JumpProblemNetwork (getDiffNetwork, params, tf, getDiffu0, nothing ,
266266 nothing )
0 commit comments