We present a set of 144 Galactic chemical evolution models applied to a Milky Way analogue, computed using four sets of low+intermediate star nucleosynthetic yields, six massive star yield compilations, and six functional forms for the initial mass function. A comparison is made between a grid of multiphase chemical evolution models computed with these yield combinations and empirical data drawn from the Milky Way's disc, including the solar neighbourhood. By means of a {chi}^2^ methodology, applied to the results of these multiphase models, the best combination of stellar yields and initial mass function capable of reproducing these observations is identified.