By employing a Bayesian approach for the statistical modelling of the Hb counts and of anaemia prevalence, our model specification via Markov chain Monte Carlo (MCMC) algorithms offers flexibility in fitting complex models and enables estimates for the whole distribution of the unknown parameters, including point and interval estimates, to be derived. This evidence concerns the gene GSTM1 and anemia (phenotype).