cput=2000:0:0
mem=256
N_PARTICLES=30;
%parpool(30);

num_iterations = 100;
model_A_path = '/mnt/bioth1/sethna/mouse_T/results/models_2015/aggregate_Bl6_WT_6Wk_thymus_unique_outframe_Mus_V_DJ_cuts_insertions_dinuc_ntbias_model/iteration_15/model.mat';
model_A_name = 'Bl6_WT_6Wk_thymus_outframe_aggregate';
model_B_path = '/mnt/bioth1/sethna/mouse_T/results/models_2015/aggregate_Bl6_WT_6Wk_thymus_unique_outframe_Mus_V_DJ_cuts_insertions_dinuc_ntbias_model/iteration_15/model.mat';
model_B_name = 'Bl6_WT_E17_thymus_outframe_aggregate';
alignments_file_base_name = 'results/alignments/aggregate_Bl6_WT_D4_blood_unique_outframe/aggregate_Bl6_WT_D4_blood_unique_outframe';
model_type = 'Mus_V_DJ_cuts_insertions_dinuc_ntbias_model';
%out_file = 'results/test.mat';

%tic; do_mixing_model_fitting(num_iterations, model_A_path, model_A_name, model_B_path, model_B_name, alignments_file_base_name, model_type, '-out_file', out_file); toc;

tic; do_mixing_model_fitting_local(num_iterations, model_A_path, model_A_name, model_B_path, model_B_name, alignments_file_base_name, model_type); toc;


%delete(gcp);
