time-lagged correlation coefficient between two time series
显示 更早的评论
I am trying to find the time-lagged correlation coefficient between two time series (two sea pressure time series at different points). I have two series of exactly the same length and with the same number of records, and I just want to see at what time lag the two series have the highest correlation. I obtained the correlation coeficient from corr.m but are looking for one that handles lag options to obtain the highest correlation at a given time.
Any help would be appreciated.
Thanks.
3 个评论
Kim
2013-2-8
I would also be very grateful for any input on this question. I have the same problem (but EEG data). I can see in the visual representation that there may be a time lag in the data of about 50 ms. I have tried xcorr with different lags, but am having trouble, 1. figuring out the significance of the return, 2. exactly how to interpret the stem plot. It should probably be obvious but I am uncertain that I have it right. Thanks!
Youssef Khmou
2013-2-8
hi, did you only use "corr.m"? or you also used "xcorr.m"?
wasswa peter
2020-7-26
Hello,I think better use R,here is the code I developed to find the time lag between meteorological drought and hydrological drought
library("ggpubr")
library("Kendall")
library("trend")
library("astsa")
my_data <- read.csv(file.choose())
my_data
names(my_data)
attach(my_data)
x <-SPI3
y <-SDI3
SPI3=ts (SPI3)
SDI3 = ts(SDI3)
ccfvalues = ccf(SPI3, SDI3)
ccfvalues
ccf (SPI3, SDI3, main = "SPI3 & SDI3 FOR STATION B",xlab = "Lag Time (months)",ylab= "Autocorrelation")
回答（2 个）
Naomi Krauzig
2019-12-12
1 个投票
Youssef Khmou
2013-2-8
编辑：Youssef Khmou
2013-2-8
hi,i propose two ways :
1)as you have to series P1, P2 of length N both:
y=xcorr(P1,P2); % Len= 2*N-1
plot(y)
2)Or logically, you can compute the quadratic error, as
( P1(i,j)-P2(i,j) )² for 1<=i,j <=N :
err=(P1-P2).^2;
f=1./err;
figure, plot(f), xlabel(' Measurments'), ylabel(' E^2'),
[Max,time_index]=max(f);
The "time_index" is when there is high correlation between p1 and p2 .
类别
在 帮助中心 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息
产品
2020-7-26
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
