有没有办法可以用 R 在数据帧上创建循环？

Question

我正在尝试创建一个循环，该循环将对我所有感兴趣的变量应用平滑函数。我使用的数据集是 R 中的 aemet 数据集，重点是温度。为了在第 10 列上应用平滑，我运行了以下代码，效果很好。

library(fda)
library(fda.usc)
data(aemet,package = "fda.usc")
tt = aemet$temp$argvals
temp = as.data.frame(aemet$temp$data,row.names=F)
range.tt = aemet$temp$rangeval
inv.temp = data.frame(t(aemet$temp$data)) # 365 x 73 matrix
names(inv.temp) = aemet$df$name # Oviedo is the 10th column of inv.temp
ovibasis5 = create.fourier.basis(rangeval = range(tt),nbasis = 5)
ovifourier5.fd = smooth.basis(argvals = tt, y = inv.temp[,i],fdParobj = ovibasis5)$fd
ovi5 = eval.fd(tt,ovifourier5.fd)
ovi5

接下来，我想对数据集中的所有 73 列进行平滑处理，所以我使用如下循环函数

for(i in 1: inv.temp){
ovibasis5 = create.fourier.basis(rangeval = range(tt),nbasis = 5)
ovifourier5.fd = smooth.basis(argvals = tt, y = inv.temp[,i],fdParobj = ovibasis5)$fd
ovi5[,i] = eval.fd(tt,ovifourier5.fd)
}
ovi5

但后来我收到此错误消息 “1:inv.temp 中的错误：NA/NaN 参数 另外：警告信息： 在 1:inv.temp 中：数值表达式有 73 个元素：只使用第一个 > > ovi5"

我做错了什么？

Answer 1

您可以将ovibasis5 放在循环之外，因为它是恒定的并且不会改变它在循环中的值。使用与inv.temp 相同的行数和列数初始化输出矩阵（ovi5）。

library(fda)
library(fda.usc)

ovibasis5 = create.fourier.basis(rangeval = range(tt),nbasis = 5)
ovi5 <- matrix(nrow = nrow(inv.temp), ncol = ncol(inv.temp))

for(i in seq_along(inv.temp)){
  ovifourier5.fd = smooth.basis(argvals = tt, y = inv.temp[,i],fdParobj = ovibasis5)$fd
  ovi5[,i] = eval.fd(tt,ovifourier5.fd)
}
#Sample output
ovi5[1:5, 1:5]

#         [,1]     [,2]     [,3]     [,4]     [,5]
#[1,] 10.99621 9.378179 7.839638 4.930860 5.370917
#[2,] 10.97734 9.360540 7.827852 4.913171 5.365705
#[3,] 10.95973 9.344325 7.817700 4.897670 5.362975
#[4,] 10.94339 9.329520 7.809155 4.884335 5.362678
#[5,] 10.92829 9.316107 7.802192 4.873138 5.364765

#complete output 
ovi5