Давайте предположим, что я измерил сигнал от t = 1 до t = 100. Однако для анализа меня интересует только период от t = 10 до t = 90. Я хочу избегать ввода y[10:90]
, потому что я хочу иметь «истинные» производные значения, а не те, которые применяются естественными граничными условиями.
library(fda)
t = 1:100
y = t^2 + rnorm(100)
basis = create.bspline.basis(c(1, 100), breaks = t, norder = 4)
fun = smooth.basisPar(t, y, basis, lambda = 0.001)$fd
Если я просто сделаю что-то вроде fun$coefs = fun$coefs[10:90, ]
, это "сломает" объект.