Model building using bi-level optimization