Commit 40d24b5f authored by Peter Jentsch's avatar Peter Jentsch
Browse files

redo community intervals simulations

parent 7cca0dc9
...@@ -105,7 +105,7 @@ function community(particles) ...@@ -105,7 +105,7 @@ function community(particles)
priors_list = [ priors_list = [
Factored(poisson_priors...), Factored(poisson_priors...),
] ]
community_error(p,d) = err_hh(p,d) + err_rest(p,d) community_error(p,d) = err_ws(p,d) + err_rest(p,d)
bayesian_estimation("community",community_error,priors_list,dists,particles; alpha = 0.98) bayesian_estimation("community",community_error,priors_list,dists,particles; alpha = 0.98)
end end
......
...@@ -7,7 +7,7 @@ const color_palette_bright = palette(:seaborn_bright) #bright color theme for th ...@@ -7,7 +7,7 @@ const color_palette_bright = palette(:seaborn_bright) #bright color theme for th
# pgfplotsx() # pgfplotsx()
distribution_keys = [Poisson] distribution_keys = [Poisson]
function plot_all() function plot_all()
fnames = ["hh","ws","rest"] fnames = ["community","hh","ws","rest"]
map(plot_estimate,fnames) map(plot_estimate,fnames)
end end
function plot_estimate(fname) function plot_estimate(fname)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment