% file for running simulations


filename1='laplace_';
filenamea='single_';
scale_match=icdf('Normal',0.75,0,sqrt(10));

    scale=-scale_match/log(0.5);
    runs=3000;
    SNP1=1;
    SNP2=2;
   
  
        vfrac=1;
        vfrac_name=num2str(vfrac*100);
        for j=1:50
            [ boot,stored_beta,accept,times_in_model,markers_in_model ] = laplace_single( runs, SNP1,SNP2, vfrac,scale);
            trial_name=num2str(j);
            new_filename = [filename1 filenamea trial_name ];
            save(new_filename);
        end
 
