使用fda包执行平滑样条拟合（Smoothing Splines）并估计曲线的方差以及绘制置信区间。

使用fda包执行平滑样条拟合（Smoothing Splines）并估计曲线的方差以及绘制置信区间。

步骤 1: 加载fda库

首先，我们需要加载fda库。确保已安装fda包，如果没有，可以使用install.packages("fda")进行安装。

library(fda)

步骤 2: 定义积分用的数值积分点和权重

在这一步中，我们定义了用于数值积分的积分点 quadpts 和相应的权重 quadwts，采用了复合辛普森法则（Composite Simpson’s Rule）来计算曲线的积分。这些积分点和权重将在后续的计算中使用。

delta <- 0.02
quadpts <- seq(0, 1, delta)
nquadpts <- length(quadpts)
quadwts <- c(1, rep(c(4, 2), (nquadpts - 1) / 2))
quadwts <- quadwts * delta / 3

步骤 3: 计算二阶基函数在积分点上的值

这一步使用 eva