Skip to content

Test `sokal_factor` in `lib.math.correlation_time`

The sokal_factor in lib.math.correlation_time is used to make the correlation more robust against noise in a timeseries. See function docs for more details. However, we do not test that the factor is working correctly. I propose two tests to ensure this:

  1. Test that method="sokal" and sokal_factor=0 is the same as method="chodera" for any timeseries.
  2. Extract the correlation time with method="chodera" for a timeseries with a smooth decaying self correlation function. Add noise to the timeseries/correlation function and extract the correlation time with method="sokal". Hopefully these two agree.