对加拿大气象数据进行函数数据分析（FDA） 完整步骤

加拿大气象数据进行函数数据分析（FDA）的多个步骤，包括温度数据的函数估计、相关性分析、主成分分析和可视化。以下是对每个部分的教程：

步骤 1: 创建基函数和定义约束

首先，我们创建了Fourier基函数 daybasis365，以用于估计温度函数。然后，我们定义了一个光滑度惩罚操作符 harmLfd，用于约束函数的光滑度。

library(fda)
# 创建Fourier基函数
daybasis365 = create.fourier.basis(c(0,365),365)

# 定义光滑度惩罚操作符
harmLfd = vec2Lfd(c(0,(2*pi/365)^2,0), c(0, 365))

步骤 2: 估计温度函数

接下来，我们使用函数 smooth.basis 对每一天的温度数据进行函数估计，并将结果存储在 tempfd 中。

tempfdPar = fdPar(daybasis365,harmLfd,1e4)
tempfd = smooth.basis(1:365,daily$tempav,tempfdPar)