R1.jl 374 Bytes
Newer Older
1
2

using MarkovProcesses
3
using Plots
4
5
6

load_model("ER")
load_automaton("automaton_F")
7
8
A_F_R1 = create_automaton_F(ER, 50.0, 75.0, 0.025, 0.05, "P")
sync_ER = A_F_R1*ER 
9
pm_sync_ER = ParametricModel(sync_ER, ("k3", Uniform(0.0, 100.0)))
10

11
12
13
14
r = automaton_abc(pm_sync_ER; nbr_particles = 1000)

histogram(r.mat_p_end', weights = r.weights, normalize = :density)
png("R1_hist.png")
15