I am trying to fit a GAM regression to my data. However I have a question. So for instance if I have
Y = S(A) + S(B) as my GAM model,
Can I include a linear model of one variables, say B to the model. So that I get
Y = S(A) + B +S(B).
I do this because I want to obtain point estimates for the variable B. In that case how would this affect the interpretation of the smooth function. Can I say my smooth function S(B) shows the deviation from my average point estimates?
I would appreciate your support