Open
Description
I perform cross-correlation from one light curve like this:
corr = CrossCorrelation(lc, lc, norm='variance')
corr is much greater than 1 at time_lags = 0.
The code in crosscorrelation.py (lines 200-201) may have a bug:
variance1 = np.var(lc1.counts) - np.mean(lc1.counts_err) ** 2
variance2 = np.var(lc2.counts) - np.mean(lc2.counts_err) ** 2
Should it be like this:
variance1 = np.var(lc1.counts) - np.var(lc1.counts_err)
variance2 = np.var(lc2.counts) - np.var(lc2.counts_err)
Metadata
Metadata
Assignees
Labels
No labels