First pass (though broken) at model estimation