Modelling and smoothing parameter estimation with multiple quadratic penalties