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
より速い代替手段はありますか?
読んでくれてありがとう。