plot an ensemble of data

plot_ensemble.m

function ave_Y = plot_ensemble(gamma,Y)
% input: gamma(N); Y(m,N)
% plot average shading for the N experiments
% return the mean

alpha_confidence = 0.1;              % significance level
if length(gamma) ~= size(Y,2)
    error('Y has the wrong size');
    return
end

[muhat,sigmahat,muci,sigmaci] = normfit(Y,alpha_confidence);
fill(vertcat(gamma(:),flipud(gamma(:))),...
    vertcat(muci(2,:)',flipud(muci(1,:)')),'r','LineStyle','none');
ave_Y = muhat;