Only in very simple situations can we describe the correlation between two time series as a single number. Often it varies both as a function of time, and of wavelength. Datasets may be uncorrelated on short timescales due to noise, but strongly correlated on larger wavelengths (eg annual).
This program will compute (using wavelets) the correlation as a function of both time and wavelength. Anticorrelation (s=-1) is displayed in blue, zero correlation (s=0) in green, and positive correlation (s=+1) in red.
The image shows the correlation between the oil price and the gold price, which is mostly red (implying strong positive correlation) except for a period between 1985 and 1995.
Note that it is valid to compare datasets that have been measured in different units using this method.
A good code, thanks
Like the quick visual results
A very nice method for comparison. In your Computers & Geosciences paper from 2008 you included dot product plots for some of the time series but the code for these plots is not included above. Could you give any pointers on implementing this code?
Clever method. Hope to see future enhancements.
Create scripts with code, output, and formatted text in a single executable document.