我有这个函数可以计算从 1 到 1600 的滞后的相关性。
t=min(length(x), length(y))
corf=function(i) {
x=x[c(1:(t-i+1))]
y=y[c(i:t)]
cor(x,y, use="pairwise.complete.obs")
}
其中 x 和 y 是每个包含大约 10,000 个值的数据集。
lg=c(1:1600)
我必须找到:
corf(lg)
为此,我正在使用:
mapply(corf, lg)
但这需要很长时间,因为我有 100 个数据集。
有没有mapply
更快的替代方案?
感谢您的阅读。