MRG, BUG: Fix biased cov estimation #7369
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
We have an inconsistency where MNE-C and MNE-Python with
method='empirical'
would compute an unbiased estimator of the covariance, but any other mode would compute a biased estimator (the maximum likelihood estimator). This PR makes us always compute the unbiased one.We could alternatively add a
bias
kwarg tocompute_covariance
andcompute_raw_covariance
with a deprecation cycle for people. But I'm not sure it's worth it because in practice the scale factor should be quite small as it isN / (N - 1)
whereN
is the number of (time) samples used to compute the covariance, which is typically quite large.