pycbc icon indicating copy to clipboard operation
pycbc copied to clipboard

Chain / acf Plots for Ensembles Should Default to Mean Over Walkers

Open farr opened this issue 7 years ago • 2 comments

It is clearer to interpret and easier to see trends when plots like chain plots (ACF plots, etc) show the trend of the mean position of the ensemble versus time rather than an individual track for each walker. (For the ACF, the ACF of the mean walker position is equivalent to the mean of the ACFs for each walker, and will have much reduced noise if the samples are really decorrelated.)

This is true for all ensemble samples (kombine, emcee, emcee_pt).

It is, of course, reasonable to add an option like --show-all-walkers to plot individual walker trajectories when the user wants to check for, e.g. stuck walkers, but it shouldn't be the default.

farr avatar Jul 26 '17 17:07 farr

PR #1832 addressed the first half of this (for plot_acf). We still need the ability to take the mean over walkers for plot_samples. @farr @cmbiwer if you'd like to take a stab at that sometime the next couple of weeks, feel free. Otherwise I'll try to get it after I get back from the PAX conference in two weeks.

cdcapano avatar Aug 08 '17 14:08 cdcapano

@cdcapano Is this still an issue?

ahnitz avatar Sep 21 '21 18:09 ahnitz