Coherence Analysis for Continuous Variables

Coherence is a measure of the degree of relationship, as a function of frequency, between two time series.

Parameters

Parameter	Description
Reference	Specifies a reference variable.
Number of Fr. Values	Number of values in the spectrum.
Window Overlap (%)	Window overlap in percent (from 0 to 90).
Window Preprocessing	Preprocessing to be done for each window before calculating the spectrum of the window.
Windowing Function	Windowing Function to be applied to each window before calculating the spectrum of the window.
Discard Incomplete Windows	Option to discard the last window if it contains less than the expected number of data points (less than NumberOfFrValues*2 points).
Use Multitaper Algorithm	Option to use multi-taper algorithm.
Time-Bandwidth Product	Time-Bandwidth Product when using multi-taper algorithm.
Number of Tapers	Number of Tapers when using multi-taper algorithm.
Show Freq. From	An option to show a subset of frequencies. Specifies minimum frequency to be displayed.
Show Freq. To	An option to show a subset of frequencies. Specifies maximum frequency to be displayed.
Calculate	Specifies whether coherence values or coherence phases are calculated.
Confidence Level (%)	Confidence level (percent) for the coherence values. Can only be calculated if a single-taper algorithm is used. See Confidence Level Calculation below for details.
Draw confidence level	An option to draw the confidence level. Confidence level is not drawn if coherence phases are specified in Calculate parameter or if a multi-taper algorithm is used.
Conf. line style	Line style for drawing confidence level.
Select Data	If Select Data is From Time Range, only the data from the specified (by Select Data From and Select Data To parameters) time range will be used in analysis. See also Data Selection Options.
Select Data From	Start of the time range in seconds.
Select Data To	End of the time range in seconds.
Interval filter	Specifies the interval filter that will be used to preselect data before analysis. See also Data Selection Options.
Add to Results / Bin left	An option to add an additional vector (containing a left edge of each frequency bin) to the matrix of numerical results.
Add to Results / Bin middle	An option to add an additional vector (containing a center of each frequency bin) to the matrix of numerical results.
Add to Results / Bin right	An option to add an additional vector (containing a right edge of each frequency bin) to the matrix of numerical results.
Send to Matlab	An option to send the matrix of numerical results to Matlab. See also Matlab Options.
Matrix Name	Specifies the name of the results matrix in Matlab workspace.
Matlab command	Specifies a Matlab command that is executed after the numerical results are sent to Matlab.
Send to Excel	An option to send numerical results or summary of numerical results to Excel. See also Excel Options.
Sheet Name	The name of the worksheet in Excel where to copy the numerical results.
TopLeft	Specifies the Excel cell where the results are copied. Should be in the form CR where C is Excel column name, R is the row number. For example, A1 is the top-left cell in the worksheet.

Summary of Numerical Results

The following information is available in the Summary of Numerical Results

Column	Description
Variable	Variable name.
Y Min	Y axis minimum.
Y Max	Y axis maximum.
Filter Length	The length of the interval filter.
Confidence Level	Confidence level for the coherence values

Algorithm

For two variables X (reference) and Y (target) the following entities are calculated. FFTs of the data segments (after preprocessing) are calculated. Then, individual and cross-densities are calculated:

Pxx = FFT(X)*Conj(FFT(X)),

Pyy = FFT(Y)*Conj(FFT(Y))

Pxy = FFT(X)*Conj(FFT(Y)).

Here Conj(z) is a complex conjugate of z. Pxx, Pyy and Pxy values are averaged across all intervals and coherence values are calculated as

Mean(Pxy)*Mean(Pxy) / (Mean(Pxx)*Mean(Pyy)).

Coherence phase values are calculated as phase of Mean(Pxy).

Calculation of the Confidence Level

Confidence levels can only be calculated if a single-taper algorithm is used.

Confidence levels are calculated as described in Kattla and Lowery (2010) (see Reference blow). The confidence level Z is calculated as

Z = 1 - pow(alpha, 1/(w*L - 1))

where

pow(x,y) returns x to the power of y,

alpha = 1 - Confidence_Level*0.01,

L is the number of overlapped windows

w is the correction due to the Hanning window tapering (see eq. (5) and (6) in Kattla and Lowery (2010)).

Reference

Kattla S, Lowery MM. Fatigue related changes in electromyographic coherence between synergistic hand muscles. Exp Brain Res. 2010 Apr; 202(1):89-99. Epub 2009 Dec 12.