R
中创建一个基矩阵,其第三个条目是
(i,j)
,其中B_j(x_i) - B_j(0)
是第三立方体B-Spline基函数。这可以实施吗?
从
B_j(.)
软件包中使用
j
函数为我带来了带有条目
bs
的基矩阵。请注意,实际值
splines
取决于我的数据,但我不是在寻找此值。我只希望我的基础函数B_j(x_i)
评估to在
B_j(0)
中,我试图通过用
B_j(t)
.替换来做到这一点。
对于更多上下文,请检查Https://stats.stackexchange.com/q/662217/464154
.
this创建基础,然后从每个列中减去第一列(尽管这似乎并没有实际将值设置为每个列的第一个元素的值...)
0
您可以创建一个执行此操作的包装函数,然后删除第一列(现在是全零):我不知道您是否要进行进一步的调整(例如,请致电t = 0
B_j(t)
the thingB_j(t) - B_j(0)
x <- 0:100
sp <- splines::bs(x, 5)
sp[] <- sweep(sp, MARGIN = 1, FUN = "-", STATS = sp[,1])