function outline_func(one, two, three, kkk, mmm) % straight forward outline of running of scripts clc display('outline') display('basic parameters') y=basics_func(one,two,three,kkk,mmm) n=y(1); m=y(2); m2=y(3); p=y(4); n_sim=y(5); r=y(6); rel=y(7); pf=y(8); c=y(9); KK=y(10); MM=y(11); prior_l=y(12); prior_x_al=y(13); prior_var_al=y(14); I=10; display('simulating fails per month') fails_per_month_f display('sampling from the prior distribution') obs(1:I) = K(1,1:I); %this gives I observed #failures in months 1:I. Only one obs for each month. % note that the first obs is 0 since 0 fails in first % month. sampling_f display('samping from the posterior using importance sampling') post_sampling_f display('sort and display results') sort_disp display('posterior predictive') post_pred display('cross validation') for i=c:I i recreated(i) = cross_val(i,obs); end hold off plot(c:I,recreated(c:I),'r') hold on plot(1:I,obs)