model_prior.jl 583 Bytes
Newer Older
1
2
3
4
5
6
7

using MarkovProcesses
load_model("ER")

test_all = true

load_model("ER")
8
9
prior1 = ParametricModel(ER, (["k2"], Uniform(2.0, 4.0)))
draw_model!(prior1)
10
11
test_all = test_all && 2.0 <= ER["k2"] <= 4.0

12
13
prior2 = ParametricModel(ER, (["k3","k2"], Product(Uniform.([2.5,6.0], [3.5,7.0]))))
draw_model!(prior2)
14
15
test_all = test_all && 2.5 <= ER["k3"] <= 3.5 && 6.0 <= ER["k2"] <= 7.0

16
17
prior3 = ParametricModel(ER, (["k3"], Uniform(10.0, 11.0)), (["k2"], Uniform(13.0, 14.0)))
draw_model!(prior3)
18
19
20
21
test_all = test_all && 10.0 <= ER["k3"] <= 11.0 && 13.0 <= ER["k2"] <= 14.0

return test_all