Chapter 5: binomial rate measurement

S=31;

F=19;

N=S+F;

f=linspace(0,1,1201)';

%%%%%%%%%%%%%%%%%%%%%%%%%

%%% 1d f distribution %%%

Lf=(S-.5)*log(f)+(F-.5)*log(1-f); Lf=Lf-max(Lf);

ciF=CIp(f,Lf);

figure(1); clf

subplot(2,1,1); hold on

plot(f,exp(Lf),'k--','LineWidth',1.75); dciF=diff(ciF(1:2)); axis([ciF(1)-.2*dciF ciF(2)+.2*dciF 0 1])

for n=1:2, plot(ciF(n)*[1 1],[0 1],'-','Color',.65*[1 1 1],'LineWidth',.7); end