我有这个函数来计算滞后的相关性,从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
有哪些替代方案更快?
感谢您的阅读。