从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)
我希望输出实数。