% file for running simulations
vfracvector=[0.5,0.45,0.4,0.35,0.3,0.25,0.2,0.15,0.1,0.09,0.08,0.07,0.06,0.05,0.04,0.03,0.02,0.01];



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


for k=1:1%length(runvector)
    scale=10;
    runs=1000000;
    
    for i=1:length(vfracvector)
        %[25,37,46,60]
        SNP=46;
        
        snp_name=num2str(SNP)
        
        vfrac=vfracvector(i);
        
        vfrac_name=num2str(vfrac*100);
        for j=1:1
            [qinfo,accept, times_in_model,markers_in_model,frequentist_log_likelihood ] = mouse_MCMC_function_new( 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

