我在R中使用了CollocInfer软件包,它在很大程度上依赖于fda软件包来估计一个常微分方程的参数。我能够成功地实现广义轮廓方法来获得这些参数估计,但我无法知道如何使用Profile.LS()函数的输出来计算ODE在t0上的估计解决方案,这将根据数据给出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)在yfd1上计算x=0值的(最简单)方法是什么?
发布于 2013-01-31 21:51:32
可以使用eval.fd以指定的参数值计算函数数据对象。
eval.fd(0,yfd1)
rep1
[1,] 3.27759https://stackoverflow.com/questions/14635105
复制相似问题