There is an increasing need for shape statistics in medical imaging to provide quantitative measures to aid in diagnosis, prognosis and therapy planning. In view of this, we describe methods for computing such statistics by utilizing a well-posed framework for representing the shape of surfaces as currents. Given this representation we can compute an atlas as a mean representation of the population and the main modes of variation around this mean. The modes are computed using principal component analysis (PCA) and applying standard correlation analysis to these allows to correlate shape features with clinical indices. Beyond this, we can compute a generative model of growth using partial least squares regression (PLS) and canonical correlation analysis (CCA). In this chapter, we investigate a clinical application of these statistical techniques on the shape of the heart for patients with repaired Tetralogy of Fallot (rToF), a severe congenital heard defect that requires surgical repair early in infancy. We relate the shape to the severity of the pathology and we build a bi-ventricular growth model of the rToF heart from cross-sectional data which gives insights about the evolution of the disease.