% file for running simulations
vfracvector1=[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];
vfracvector2=[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_multi_eq_';
filename1i='_and_';
filename2='_variance_fraction_';
filename3='_trial_';


for k=1:1%length(runvector)
    scale=10;
    runs=2000;
    run_name=num2str(runs);
    for i=1:length(vfracvector1)
        
        SNP1=37;
        SNP2=46;
        snp_name1=num2str(SNP1)
        snp_name2=num2str(SNP2)
        vfrac1=vfracvector1(i);
        vfrac2=vfracvector2(i);
        vfrac_name=num2str(vfrac1*100);
        for j=1:25
            [qinfo,stored_beta,accept,qvarmonitor,log_R_track, times_in_model,markers_in_model,qtlsim_a,qtlsim_b ] = mouse_MCMC_function_multi( runs, SNP1,SNP2, vfrac1, vfrac2,scale );
            trial_name=num2str(j);
            new_filename = [filename1 snp_name1  filename1i snp_name2 filename2 vfrac_name filename3 trial_name ];
            save(new_filename);
        end
    end
end

