Interpreting pigeons' standard output

During the execution of parallel tempering, interim diagnostics can be computed and printed to standard out at the end of every iteration (this can be disabled using show_report = false):

using Pigeons

pigeons(target = toy_mvn_target(100))

┌ Info: Neither traces, disk, nor online recorders included.
│    You may not have access to your samples (unless you are using a custom recorder, or maybe you just want log(Z)).
└    To add recorders, use e.g. pigeons(target = ..., record = [traces; record_default()])
────────────────────────────────────────────────────────────────────────────
  scans        Λ        time(s)    allc(B)  log(Z₁/Z₀)   min(α)     mean(α)
────────── ────────── ────────── ────────── ────────── ────────── ──────────
        2        7.5    2.3e-05    1.1e+04       -122   2.66e-15      0.167
        4       5.56   5.18e-05   1.18e+04       -119    2.6e-07      0.382
        8       5.82   5.88e-05   1.61e+04       -114    0.00137      0.353
       16       7.02   9.93e-05   2.15e+04       -116   0.000205      0.219
       32       7.14   0.000184   2.76e+04       -112     0.0379      0.206
       64       7.18   0.000357   3.96e+04       -116     0.0442      0.202
      128        7.2   0.000713   5.86e+04       -114      0.142        0.2
      256       7.21    0.00139   5.66e+04       -115      0.151      0.199
      512       7.14    0.00274   7.45e+04       -115       0.19      0.207
 1.02e+03       7.12    0.00536   7.02e+04       -115      0.173      0.209
────────────────────────────────────────────────────────────────────────────

The functions called to emit each of these can be found at all_reports(). Some key quantities:

• Λ: the global communication barrier, as described in Syed et al., 2021 and estimated using the sum of rejection estimator analyzed in the same reference. Syed et al., 2021 also developed a rule of thumb to configure the number of chains: PT should be set to roughly 2Λ.
• time and allc: the time (in second) and allocation (in bytes) used in each round.
• log(Z₁/Z₀): the stepping_stone() estimator for the log of the normalization constant, see the documentation page on approximation of the normalization constant.
• min(α) and mean(α): minimum and average swap acceptance rates over the PT chains.

Additional statistics can be shown when more recorders are added. For example, to accumulate other constant-memory summary statistics:

pigeons(target = toy_mvn_target(100), record = record_online(), explorer = AutoMALA())

───────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
  scans     restarts      Λ        time(s)    allc(B)  log(Z₁/Z₀)   min(α)     mean(α)    max|ρ|     mean|ρ|    min(αₑ)   mean(αₑ)
────────── ────────── ────────── ────────── ────────── ────────── ────────── ────────── ────────── ────────── ────────── ──────────
        2          0        6.2      0.145   1.21e+07       -131   3.84e-13      0.311          1          1      0.284      0.536
        4          0       6.21    0.00174    8.1e+04       -122    0.00023       0.31      0.827      0.345       0.41      0.539
        8          0       7.19    0.00318   8.67e+04       -117     0.0114      0.201      0.617      0.361      0.441      0.582
       16          0       7.14     0.0063   1.86e+05       -119    0.00252      0.207      0.586      0.268      0.459      0.523
       32          0        7.1     0.0124   2.55e+05       -111     0.0483      0.211       0.63       0.24      0.451      0.511
       64          0       7.09     0.0242    3.3e+05       -115      0.083      0.212      0.521      0.292      0.474      0.535
      128          0       7.27     0.0481   3.82e+05       -115      0.112      0.192      0.462      0.303      0.468      0.533
      256          1       7.34     0.0959   6.24e+05       -113      0.106      0.184      0.436      0.306      0.487      0.534
      512          3       7.08      0.191   9.53e+05       -115      0.164      0.214      0.449      0.332      0.474      0.532
 1.02e+03         11       7.18      0.381   1.57e+06       -115      0.161      0.202      0.472      0.332       0.48      0.534
───────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────

• max|ρ| and mean|ρ|: maximum and average (across chains) correlation of the random variables $L^t_i = V(X_i)$ and $L^{t+1}_i = V(X_i)$ where $V = \log \pi_N / \pi_1$, $X_i \sim \pi_{\beta_i}$, and $t, t+1$ are indices just before and after a call to step!().
• min(αₑ) and mean(αₑ): minimum and average (across chains) of the explorer's acceptance rates.

Programmatic access

The tables described above can also be accessed as a DataFrame via:

using Pigeons
pt = pigeons(target = toy_mvn_target(100))
pt.shared.reports.summary

Detailed statistics can be accessed via these DataFrames. For example to obtain mean swap acceptance probabilities for each round and pair of communicating chains, use:

first(pt.shared.reports.swap_prs, 20)

Creating a plot from this:

using StatsPlots
plotlyjs()
my_plot =  @df pt.shared.reports.swap_prs StatsPlots.plot(:round, :mean, group = :first)
StatsPlots.savefig(my_plot, "swap_prs.html");