% file for running simulations
vfracvector=[0.45,0.4,0.35,0.3,0.25,0.2,0.15,0.1,0.05];
runvector=400;

filename1='mousesnp_';
filename2='_variance_fraction_';
filename3='_trial_';

for k=1:1
    scale=10;
    runs=runvector(k);
    
    for i=1:length(vfracvector)
        
        SNP=25;%60;
        run_name=num2str(SNP);
        vfrac=vfracvector(i);
        vfrac_name=num2str(vfrac*100);
        for j=1:20
            [qinfo,stored_beta,accept,qvarmonitor,log_R_track, times_in_model,markers_in_model,frequentist_log_likelihood ] = mouse_MCMC_function_new( runs, SNP, vfrac,scale );
            trial_name=num2str(j);
            new_filename = [filename1 run_name filename2 vfrac_name filename3 trial_name ];
            save(new_filename);
        end
    end
end

