%sort and display the predicted values. sorted_post_pred(:,count) = sort(post_pred_sample(:,count)); upper(count)=sorted_post_pred(int32(MM*.025),count); lower(count)=sorted_post_pred(int32(MM*.975),count); post_mean = mean(post_pred_sample); %post_v = mean(post_var_sample); post_prediction(count+1) = post_mean(count+1); %post_var(count+1) = post_v(count+1) display('#failed in month I+1, #predicted, 95%upper, 95%lower') [obs(1,(count+1)) post_mean(count+1) upper(count) lower(count)]