如何从R中的fdsmooth对象提取内部结的子集？

Question

假设我测量了从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, ]之类的操作，这会“破坏”对象。