Build the design and the penalty matrix for models involving penalised splines based on a formula and a data set