r - mgcv中的循环自适应样条

Question

我想使用循环预测器将 GAM 拟合到数据中，并且在循环的特定部分中大部分摆动。
有没有办法使循环样条(bs = 'cc' 或 'cp')自适应？ ...或者，等效地，使自适应样条 (bs = 'ad') 循环？

Answer 1

是的;这在 mgcv 中的自适应样条基础中已经允许。
自适应样条中的默认基础是 P 样条。您可以切换到该类型样条的循环版本或使用循环三次样条。
要使其工作，您必须将信息传递给 xt平滑函数的参数，同时离开 bs = "ad"为适应基础。
对于循环 P 样条，你会做

y ~ s(x, bs = "ad", xt = list(bs = "cp"))

对于循环三次样条，您将使用

y ~ s(x, bs = "ad", xt = list(bs = "cc"))

xt参数通常用于此类事情，其中​​基础具有可以配置的其他选项。 fs基础相似，其中 xt允许您控制用于随机平滑的基础。
要查看的另一个参数是 m ;哪里 k指定实际平滑的基础维度，您可以使用 m为自适应部分设置基础，更高 m表明在 x 范围内惩罚的更多潜在变化，就像 k允许更多的摆动 x .
这些细节在 ?smooth.construct.ad.smooth.spec 中讨论(或 ?adaptive.smooth 作为该页面的更简单快捷方式。)