-
Notifications
You must be signed in to change notification settings - Fork 0
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
ENH: lower-level functions for the adjustment of SNR #9
Conversation
tests/test_snr.py
Outdated
vertices = [[0, 1], [0, 1]] | ||
stc = prepare_stc(vertices) | ||
variance = get_sensor_space_variance(stc, fwd, filter=False) | ||
assert variance >= 0, "Variance should be non-negative" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think that this check is not sufficient - all values (stc and leadfield) are squared so the result should be positive regardless of which vertices are used. Please define the leadfield values yourself and calculate the expected variance with all sources and a subset of sources. Feel free to reduce the dimensionality to make the calculations easier (2x2 leadfield and 1-source stc should already be sufficient to check that the second source does not play a role).
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
MEG/EEG aspect is probably not so important here, both are numbers in the end
For review