We perform automated bulge + disc decomposition on a sample of similar to 7500 galaxies from the Galaxy And Mass Assembly (GAMA) survey in the redshift range of 0.002 < z < 0.06 using Structural Investigation of Galaxies via Model Analysis, a wrapper around GALFIT3. To achieve robust profile measurements, we use a novel approach of repeatedly fitting the galaxies, varying the input parameters to sample a large fraction of the input parameter space. Using this method, we reduce the catastrophic failure rate significantly and verify the confidence in the fit independently of chi(2). Additionally, using the median of the final fitting values and the 16th and 84th percentile produces more realistic error estimates than those provided by GALFIT, which are known to be underestimated. We use the results of our decompositions to analyse the stellar mass - half-light radius relations of bulges, discs and spheroids. We further investigate the association of components with a parent disc or elliptical relation to provide definite z = 0 disc and spheroid M-star-R-e relations. We conclude by comparing our local disc and spheroid M-star - R-e to simulated data from EAGLE and high-redshift data from Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey-Ultra Deep Survey. We show the potential of using the M-star-R-e relation to study galaxy evolution in both cases but caution that for a fair comparison, all data sets need to be processed and analysed in the same manner.