Data come from calculations in post “Merging Bayes results from coarse and detailed calculations”
>> [m,ind_alfa]=min(abs(log10(.05)-logalfa_merge))
m =
0.0010
ind_alfa =
18
>> [m,ind_beta]=min(abs(log10(1.5)-logbeta_merge))
m =
0.0761
ind_beta =
16
>> close all
>> plot(pot_merge,prob_merge(:,ind_alfa,ind_beta))
>> figure
>> plot(pot_merge,log10(prob_merge(:,ind_alfa,ind_beta)))
>> prob_linea=prob_merge(:,ind_alfa,ind_beta)/sum(prob_merge(:,ind_alfa,ind_beta));
>> expected=sum(prob_linea’.*pot_merge)
expected =
1.3382
>> sd=sqrt(sum(prob_linea’.*(pot_merge-expected).^2))
sd =
0.0866
Leave a Reply