I use the model for growth that is used in Figure 4.
fig_pap_coli_11(1,0)
(and use internal variables from now on)
New metric:
p=calibrametrica(model_growth,flujosexp_finales,1000,100,Delta);
K>> hist(p)
K>> hist(log10(p))
K>> sum(p<.05)
ans =
2
Looks good, but I will repeat it with more repetitions.
Now I test the old metric:
p2=calibrametrica(model_growth,flujosexp_finales,1000,100,Delta,2);
K>> hist(p2)
K>> hist(log10(p2))
K>> sum(p2<.05)
ans =
0
K>> sum(p2<.1)
ans =
3
It seems there was a bias. I will also repeat it with more repetitions.