r - 一般化された加算モデルから滑らかな関数の表現を抽出する方法

Question

一般化された加法モデルで適合した滑らかな関数の式を抽出する際に問題があります。

次のような加法的モデルを適合させました。

library(mgcv)
model <- gam(y ~ x1 + x2 + s(x3, bs = "cr"), data = newdata)

の近似平滑関数をプロットすることはできますがs(x3)、その正確な式を取得し、その導関数を取得したいと考えています。どうすればそれを達成できますか？