从R的FDA / CollocInfer包中的平滑数据中检索值

Question

我在R中使用了很大程度上依赖于fda包的CollocInfer包来估计常微分方程的参数。我能够成功实现广义分析方法来获得这些参数估计，但我无法弄清楚如何使用Profile.LS（）函数的输出来评估在t0处对ODE的估计解，这将给出估计基于数据的ODE的初始条件。

编辑我之前的示例太简单了，并没有捕获到问题。以下是修改后的示例。

简单示例假设我有数据x，y，我使用B样条基础平滑

require(fda)
x <- seq(0, 25)
y <- 3 + sin(x) + rnorm(1, 0, 1)
b.basis <- create.bspline.basis(range(x), nbasis=8, norder=4)
smooth.y <- smooth.basis(x, y, b.basis)
yfd <- smooth.y$fd
yfdPar <- fdPar(yfd, 1, 1)
yfd1 <- smooth.fd(yfd, yfdPar)

在x = 0时计算yfd1值的（最简单）方法是什么？

Answer 1

您可以使用eval.fd

评估指定参数值的功能数据对象

eval.fd(0,yfd1)
        rep1
[1,] 3.27759