How do I calculate a sliding correlation between two time series
10 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
I would like to compute the sliding Pearson's correlation coefficient between two columns of data which represent physiologic measurements over time, at 5-second intervals. For each point in time, I want to calculate the correlation between the two columns over the preceding 5 minutes.
0 comentarios
Respuestas (1)
dpb
el 23 de Mzo. de 2018
Presuming no missing observations and beginning at the minute, 5 min/[5 sec/sample/60 sec/min]--> 60 samples. Given the two signals, X, Y, then also assuming they cover an integer number of minutes
X=reshape(X,60,[]);
Y=reshape(Y,60,[]);
R=diag(corr(X,Y));
If mod(length(X),60) ~=0 then either discard the end effect elements or treat those partial sections individually.
The above is somewhat wasteful in that corr calculates the full pairwise correlations between columns of X, Y and discards everything excepting the diagonal elements, but unless the series are quite long performance is probably not an issue.
One can do the two columns in loop just pairwise otherwise, of course.
3 comentarios
dpb
el 23 de Mzo. de 2018
Editada: dpb
el 23 de Mzo. de 2018
That's long to look at in a text editor, granted, but not particularly large for data; I don't think compute time will be an issue.
As for MV's; hmmm...have to take that under some consideration with corr; it'd mess up the above to try to use straight logical addressing with isfinite because lengths would vary.
Do you have a preferred way to consider MVs? If were to infill, then the above would work; how much that might affect the result would depend heavily on how prevalent and just what the signals look like.
The other simple alternative I see at the moment would be to revert to the loop on paired columns that could look something like...
nc=size(X,2); % # columns (after reshape() )
R=zeros(nc,1); % preallocate for output
for i=1:nc % for columns
ix=all(isfinite([X(:,i) Y(:,i)]),2); % rows neither NaN
R(i)=corr(X(ix,i),Y(ix,i));
end
Ver también
Categorías
Más información sobre Spline Postprocessing en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!