SmpCouts.m

function [MAP]=SmpCouts(Seq,samples,Mle,N)
labels={'μ','K','Alpha','c','p'};
minA=min(samples);
maxA=max(samples);
if  nargin < 4
    N=50;
end
dM=(maxA-minA)/N;
bin=minA:dM:maxA+dM;
nbin=length(bin);
for i=1:nbin
    num=histcounts(samples,bin);
end

plot(bin(1:end-1),num,'LineWidth',1.7,'Color','k');
hold on
ind=find(num==max(num),1);
MAP=bin(ind);
fprintf('The Mostlikely Value of parameter %d is %.4f \n',Seq,MAP)
xline(MAP,'LineWidth',2,'Color','r','LineStyle','-.');
if ~isempty(Mle)
    hold on
    Mle=Mle(Seq);
    xline(Mle,'LineWidth',2,'Color',[0.5020    0.5020    0.5020],'LineStyle','--');
end
Fun_defaultAxes; box on;
title((labels(Seq)));
xlim('tight');
end