model_prior.jl 550 Bytes
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21

using MarkovProcesses
load_model("ER")

test_all = true

load_model("ER")
prior1 = ModelPrior(ER, (["k2"], Uniform(2.0, 4.0)))
draw!(prior1)
test_all = test_all && 2.0 <= ER["k2"] <= 4.0

prior2 = ModelPrior(ER, (["k3","k2"], Product(Uniform.([2.5,6.0], [3.5,7.0]))))
draw!(prior2)
test_all = test_all && 2.5 <= ER["k3"] <= 3.5 && 6.0 <= ER["k2"] <= 7.0

prior3 = ModelPrior(ER, (["k3"], Uniform(10.0, 11.0)), (["k2"], Uniform(13.0, 14.0)))
draw!(prior3)
test_all = test_all && 10.0 <= ER["k3"] <= 11.0 && 13.0 <= ER["k2"] <= 14.0

return test_all