The classical Chin-Hosford-Mendorf (CHM) criterion for deformation twinning is used widely in polycrystal plasticity computer simulations. It assumes that twin nuclei are abundant, and that twin propagation and growth control the realisation of deformation twins. However, recent experimental studies reported in the literature have found that nominally unfavourable twin systems are activated in grains. This has been attributed to twin nucleation in some of the unfavourable twinning systems with low Schmid factors, and its suppression in some of the nominally more favourable ones. Presently, this explanation is quantitatively examined. Full and relaxed constraint versions of the Taylor, ALAMEL, and binary tree based models, all implementing the CHM criterion are used to simulate uniaxial compression of a zirconium billet. Model predictions are compared with experimentally measured twinning statistics reported in the literature for a Zr polycrystal. The ALAMEL and binary tree models, which explicitly represent intergranular interactions, are found to capture the twinning statistics well. These observations suggest that the CHM criterion is adequate to capture twinning in Zr, provided intergranular interactions are represented in the model used to interpret the experiments. © 2017 Elsevier Ltd. All rights reserved.