1

我有这个函数可以计算从 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更快的替代方案?

感谢您的阅读。

4

0 回答 0