mgcvパッケージを使用して、次の方法で多項式スプラインをデータに適合させています。
x.gam <- gam(cts ~ s(time, bs = "ad"), data = x.dd,
family = poisson(link = "log"))
はめあいの機能的な形を抽出しようとしています。x.gam
でありgamObject
、ドキュメントを読んでいますが、適合関数を手動で再構築するのに十分な情報が見つかりませんでした。
x.gam$smooth
結び目が配置されているかどうかに関する情報が含まれています。x.gam$coefficients
スプライン係数を与えますが、多項式スプラインが使用される順序がわからず、コードを調べても何も明らかになりませんでした。
フィットを手動で再構築できるように、使用されるノット、係数、および基底を抽出する適切な方法はありますか?