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

filename1='cc_SNP_';
filenamea='_and_'
filename2='_variance_fraction_';
filename3='_trial_';

for k=1:1
    scale=10;
    runs=2000;
    SNP1=SNPvector(1);
    SNP2=SNPvector(2);
    SNP_name1=num2str(SNP1);
    SNP_name2=num2str(SNP2);
    for i=1:length(vfracvector)
        vfrac=vfracvector(i);
        vfrac_name=num2str(vfrac*100);
        for j=1:40
            [ boot,stored_beta,accept,times_in_model,markers_in_model ] = human_MCMC_norm( runs, SNP1,SNP2, vfrac,scale);
            trial_name=num2str(j);
            new_filename = [filename1 SNP_name1 filenamea SNP_name2 filename2 vfrac_name filename3 trial_name ];
            save(new_filename);
        end
    end
end

for k=1:1
    scale=10;
    runs=2000;
    SNP1=SNPvector(1);
    SNP2=SNPvector(2);
    SNP_name1=num2str(SNP1);
    SNP_name2=num2str(SNP2);
    for i=1:1
        vfrac=0.15;
        vfrac_name=num2str(vfrac*100);
        for j=31:40
            [ boot,stored_beta,accept,times_in_model,markers_in_model ] = human_MCMC_norm( runs, SNP1,SNP2, vfrac,scale);
            trial_name=num2str(j);
            new_filename = [filename1 SNP_name1 filenamea SNP_name2 filename2 vfrac_name filename3 trial_name ];
            save(new_filename);
        end
    end
end
