function [ merge_all_drugs,drugs_to_include,drug_threshold,report_skew ] = set_model_options( model_number )
% Set the model parameters based on the model number
switch model_number
    case 1
        merge_all_drugs = 0;
        drugs_to_include = [1 1 1];
        drug_threshold = 1;
        report_skew = 0;
    case 2
        merge_all_drugs = 1;
        drugs_to_include = [0 0 0];
        drug_threshold = 1;
        report_skew = 0;
    case 3
        merge_all_drugs = 1;
        drugs_to_include = [1 1 1];
        drug_threshold = 1;
        report_skew = 0;
    case 4
        merge_all_drugs = 0;
        drugs_to_include = [1 0 0];
        drug_threshold = 1;
        report_skew = 0;
    case 5
        merge_all_drugs = 0;
        drugs_to_include = [0 1 0];
        drug_threshold = 1;
        report_skew = 0;
    case 6
        merge_all_drugs = 0;
        drugs_to_include = [0 0 1];
        drug_threshold = 1;
        report_skew = 0;
    case 7
        merge_all_drugs = 0;
        drugs_to_include = [1 1 1];
        drug_threshold = 5;
        report_skew = 0;
    case 8
        merge_all_drugs = 0;
        drugs_to_include = [1 1 1];
        drug_threshold = 1;
        report_skew = 1;
    case 9
        merge_all_drugs = 0;
        drugs_to_include = [0 0 0];
        drug_threshold = 1;
        report_skew = 1;
    case 10
        merge_all_drugs = 1;
        drugs_to_include = [1 1 1];
        drug_threshold = 1;
        report_skew = 1;
    otherwise
        disp('error with model number');
end

end