Ferroic domain walls are currently investigated by several state-of-the art techniques in order to get a better understanding of their distinct, functional properties. Here, Principal Component Analysis (PCA) of Raman maps is used to study ferroelectric domain walls (DWs) in LiNbO3 and ferroelastic DWs in NdGaO3. It is shown that PCA allows to quickly and reliably identify small Raman peak variations at ferroelectric DWs and that the value of a peak shift can be deduced - accurately and without a-priori - from a first order Taylor expansion of the spectra. The ability of PCA to separate the contribution of ferroelastic domains and DWs to Raman spectra is emphasized. More generally, our results provide a novel route for the statistical analysis of any property mapped across a DW.