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