Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

The calculated GC value has a negative number. Is it normal? #89

Open
neurochenhe opened this issue Dec 19, 2024 · 2 comments
Open

The calculated GC value has a negative number. Is it normal? #89

neurochenhe opened this issue Dec 19, 2024 · 2 comments

Comments

@neurochenhe
Copy link

Dear developer, thank you for such a great toolkit for brain region connectivity calculation.

I am using pipeline_connectivity.m and calculating the GC in the example dataset, I found that the value stored in the EEG.roi.GC variable after the calculation has a negative value, is this an anomaly, and should the Granger causality value be greater than 0, and should not have a negative value.

When plotting the schematic via pop_roi_connectplot(EEG, 'measure', 'GC', 'plotmatrix', 'on'); you can see that the colorbar starts with a negative value. Thank you for your answer!

image

@nguyen-td
Copy link
Collaborator

Hey there,

Thanks for using our toolbox!

@stefanhaufe Please correct me if I am wrong. GC is a directional connectivity metric, meaning that it indicates a flow from region A to B. A negative value then indicates a flow from B to A.

@stefanhaufe
Copy link
Collaborator

It is correct that GC is a non-negative quantity. However, we are not looking at original GC but so-called net GC, which is the difference GC_i->j - GC_j->i . This quantity is always antisymmetric, that is netGC_i->j = - netGC_j->i and this is what is shown in the plot.

netGC is already much more robust wrt. source mixing than original GC. To be fully robust, we recommend using time reversed GC (TRGC), which is the difference netGC_i->j - netGCR_i->j, where netGCR_i->j is netGC estimated on time reversed signals.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants