mgcv包使用张量积时，如何手动设置节点位置？

你的代码存在两个关键问题，导致报错：

1. 变量名错误：knots列表里误指定了响应变量y的节点，但张量积平滑针对的是x和z，应该把y=c(...)改为z=c(...)。
2. 节点数量与基函数数量不匹配：te()默认给每个变量的基函数数量k=5，但你只提供了3个节点，而cr样条要求节点数量必须等于指定的k值。

正确代码示例 #

set.seed(0)
n <- 500
x <- runif(n)
z <- runif(n)
y <- runif(n)
# 指定每个变量的k值与节点数量一致（3个节点对应k=3），同时修正knots里的变量名
gam(y ~ te(x, z, k = c(3, 3)), knots = list(x = c(0.2, 0.5, 0.8), z = c(0.2, 0.5, 0.8)))

补充说明 #

• 如果需要使用更多基函数（比如默认的k=5），你需要给每个变量提供对应数量的节点，例如：

  gam(y ~ te(x, z), knots = list(x = seq(0, 1, length.out = 5), z = seq(0, 1, length.out = 5)))
  

• 节点建议覆盖变量的取值范围（这里x和z都是0-1，包含边界值能减少样条外插误差）。

内容的提问来源于stack exchange，提问作者Lorraine