R3.jl 396 Bytes
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

using MarkovProcesses
using Plots

load_model("ER")
load_automaton("automaton_F")
A_F_R3 = create_automaton_F(ER, 25.0, 50.0, 0.05, 0.075, "P")
sync_ER = A_F_R3*ER 
pm_sync_ER = ParametricModel(sync_ER, ("k3", Uniform(0.0, 100.0)))

r = automaton_abc(pm_sync_ER; nbr_particles = 1000)

histogram(r.mat_p_end', weights = r.weights, normalize = :density, xlims = (0.0, 100.0))
png("R3_hist.png")