Dear R users, I am using R-function gam(gam) to fit additive models, more specifically partially linear models. Say, Y~x1+x2+s(z,df=d). Is there a way to extract the smoother matrix for the smoothing term z, or is there a way to extract the overall hat matrix H such that Y_hat=HY. Thank you in advance! Sincerely, Song ******************************************************************** Song Liu School of Statistics 313 FordH 224 Church St. SE Minneapolis,MN 55455