Get predicted values from GAM models

I used to rely on the default plot produced  by the GAM command; now I realized that it is quite simple to produce and plot predicted predicted values using the "predict.gam()" function in the "mcgv" package.