Add the standard mgcv approach with bs='re', and gam.vcomp
Add the standard mgcv approach with bs='re', and gam.vcomp