Blazars optical emission is generally dominated by relativistic jets, although the host galaxy, accretion disc, and broad-line region (BLR) may also contribute significantly. Disentangling their contributions has been challenging for years due to the dominance of the jet. To quantify the contributions to the spectral variability, we use the statistical technique for dimensionality reduction non-negative matrix factorization on a spectroscopic data set of 26 γ-ray blazars. This technique allows to model large numbers of spectra in terms of a reduced number of components. We use a priori knowledge to obtain components associated with meaningful physical processes. The sources are classified according to their optical spectrum as host-galaxy dominated BL Lac objects (BL Lacs), BL Lacs, or flat spectrum radio quasars (FSRQs). Host-galaxy sources show less variability, as expected, and bluer-when-brighter (BWB) trends, as the other BL Lacs. For FSRQs, more complicated colour-flux behaviours are observed: redder-when-brighter for low states saturating above a certain level and, in some cases, turning to BWB. We are able to reproduce the variability observed during 10 yr using only two to four components, depending on the type. The simplest scenario corresponds to host-galaxy blazars, whose spectra are reconstructed using the stellar population and a power law (PL) for the jet. BL Lac spectra are reproduced using from two to four PLs. Different components can be associated with acceleration/cooling processes taking place in the jet. The reconstruction of FSRQs also incorporates a QSO-like component to account for the BLR, plus a very steep PL, associated with the accretion disc.