是否有用于分离平滑功能数据的R函数?

时间:2019-05-14 18:29:36

标签: r

从1971年到2000年,我绘制了两条平滑的函数曲线,分别用8和349等距节的B样条近似。现在,我想将这些函数分成较短的间隔,例如每年,并计算这些功能之间的半指标。 Rstudio中有此任务的功能吗?

我试图在fd对象中更改rangeval但在功能上 semimetric.basis(fd_1,fd_2,nderiv = 0)发生错误

  

“ bsplineS中的错误(evalarg,breaks,norder,nderiv,returnMatrix):   结没有增加”。

basis <- create.bspline.basis(rangeval = c(1971, 2000),
                              nbasis = 355,
                              norder = 8) 
Lfdobj_f <- int2Lfd(4)
Sm_f <- smooth.basisPar(argvals = seq(from = 1971, to = 2000,length  = 349),              
                        y = cbind(func_1, func_2),
                        fdobj = basis,
                        Lfdobj = Lfdobj_f,
                        lambda = 1e-9)

Sm_f[1, 0]$fd$basis$rangeval <- c(1971, 1972)
Sm_f[2, 0]$fd$basis$rangeval <- c(1971, 1972)

output <- semimetric.basis(Sm_f[1, 0]$fd, Sm_f[2, 0]$fd, nderiv = 0)

我希望输出实数。

0 个答案:

没有答案