% file for running simulations
vfracvector=[0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1,0.05];
SNPvector=[87,176];

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

for k=1:length(SNPvector)
    scale=10;
    runs=2000;
    SNP=SNPvector(k);
    SNP_name=num2str(SNP);
    for i=1:length(vfracvector)
        vfrac=vfracvector(i);
        vfrac_name=num2str(vfrac*100);
        for j=26:60
            [ boot,stored_beta,accept,times_in_model,markers_in_model ] = human_MCMC_norm( runs, SNP, vfrac,scale);
            trial_name=num2str(j);
            new_filename = [filename1 SNP_name filename2 vfrac_name filename3 trial_name ];
            save(new_filename);
        end
    end
end

