We should change / integrate the following two things:

• When plot_1Dsmooth is called on a model with log link, it by default uses Odds Ratio as ylab. This should be changed to exp(Effect), since Odds Ratio is not true for all log link models.
• Allow the user to pass a numeric vector to the select argument. The function should then create multiple individual smooth effect plots, side by side.