Elemental abundances of stars are the result of the complex enrichment history of their galaxy. The interpretation of observed abundances requires flexible modeling tools to explore and quantify the information about Galactic chemical evolution (GCE) stored in such data. I will present Chempy, a newly developed code for GCE modeling, representing a parametrized open one-zone model within a Bayesian framework. Unlike established approaches, Chempy can sample the posterior probability distribution in the full model parameter space and test data-model matches for different nucleosynthetic yield sets. As an illustrative application, I will show that interesting parameter constraints result from only the ages and elemental abundances of Sun, Arcturus and the present-day interstellar medium (ISM).