function [mu_dist, mu_cdf] = PostMuVals(preds_in,data) preds = preds_in(2:end); %exclude lowest linewidth to make sure preds and dpres are same length NN = length(data); %number of data points %mean and mean square of log data Xbar = sum(log(data))./NN; X2bar = sum(log(data).^2)./NN; %paramters of t-distribution mumu = Xbar; numu = NN-1; sigmu = sqrt((X2bar - Xbar.^2)./(NN-1)); mumu %calculate pdf values from t-dist pd = makedist('tLocationScale','mu',mumu,'sigma',sigmu,'nu',numu); mu_dist = pdf(pd,log(preds)); %calculate cdf logdpreds = log(preds_in(2:end)) - log(preds_in(1:end-1)); mu_cdf = cumsum(mu_dist.*logdpreds); end