This paper introduces a practical data-driven method to discriminate among large-scale kinetic reaction models. The approach centers around a computable measure of model/data mismatch. We introduce two provably convergent algorithms that were developed to accommodate large ranges of uncertainty in the model parameters. The algorithms are demonstrated on a simple toy example and a methane combustion model with more than 100 uncertain parameters. They are subsequently used to discriminate between two models for a contemporarily studied biological signaling network.