Crystal structure prediction is a well-suited application for a statistical modelling approach because of the large numbers of crystal structures, all with identical chemical composition and differing only in the arrangement of their constituent molecules. For the three molecules studied here, each crystal structure landscape contains between 388 and 526 distinct crystal structures in the energy range studied and we show that PBE0 lattice energies can be predicted to an accuracy of 1 kJ mol -1 with between 4.2 and 6.8% of the computational cost of the full PBE0 calculations. This cost reduction is achieved by using a larger set of training structures calculated at a lower level of theory, here the GGA functional PBE, after which a very small number of PBE0 calculations are required for the second level of the model. As the cost of GGA and hybrid calculations scale roughly as O(n 3 e ) and O(n 4 e ), respectively, this feature is of particular importance for larger molecules and could bring substantial cost savings.