In the statistical framework of likelihood-free inference, the posterior distribution of model parameters is explored via simulation rather than direct evaluation of the likelihood function, permitting inference in situations where this function is analytically intractable. We consider the problem of estimating cosmological parameters using measurements of the weak gravitational lensing of galaxies; specifically, we propose the use a likelihood-free approach to investigate the posterior distribution of some parameters in the ΛCDM model upon observing a large number of sheared galaxies. The choice of summary statistic used when comparing observed data and simulated data in the likelihood-free inference framework is critical, so we work toward a principled method of choosing the summary statistic, aiming for dimension reduction while seeking a statistic that is as close as possible to being sufficient for the parameters of interest.