pith. sign in

arxiv: 2502.01093 · v1 · submitted 2025-02-03 · 🌀 gr-qc · astro-ph.HE· astro-ph.IM

A practical Bayesian method for gravitational-wave ringdown analysis with multiple modes

Pith reviewed 2026-05-23 03:51 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HEastro-ph.IM
keywords gravitational wavesringdown analysisquasi-normal modesBayesian inferenceimportance samplingblack hole spectroscopycomputational accelerationFIREFLY
0
0 comments X

The pith

FIREFLY analytically marginalizes amplitudes and phases of multiple quasi-normal modes to cut ringdown inference time from hours to minutes while keeping posteriors and evidence consistent.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces an algorithm called FIREFLY to make Bayesian analysis of black-hole ringdown signals feasible when multiple quasi-normal modes must be included. It achieves this by analytically integrating out the amplitude and phase parameters of each mode, drawing on the structure of the F-statistic, then using importance sampling to recover the full posterior over the remaining parameters. The result is a reduction in computation time from hours to minutes that grows with the number of modes, while the recovered posteriors and evidence values stay consistent with those from standard nested sampling or MCMC runs. A reader would care because upcoming detectors are expected to resolve several modes in a single event, turning ringdown data into a practical tool for measuring black-hole properties at high precision.

Core claim

FIREFLY analytically marginalizes the amplitude and phase parameters of QNMs to reduce the computational cost and speed up the full-parameter inference from hours to minutes, while achieving consistent posterior and evidence. Rigorously based on the principle of Bayesian inference and importance sampling, the method is statistically interpretable, flexible in prior choice, and compatible with various advanced sampling techniques.

What carries the argument

The FIREFLY algorithm, which analytically marginalizes the amplitude and phase parameters of quasi-normal modes and applies importance sampling to the remaining parameters.

If this is right

  • The computational speedup grows larger as the number of included quasi-normal modes increases.
  • The method remains compatible with any advanced sampling technique used on the reduced parameter space.
  • Prior choices on the remaining parameters stay fully flexible without altering the marginalization step.
  • The approach supplies a statistically interpretable route to faster ringdown analyses that scales with detector sensitivity.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same marginalization step could be adapted to other gravitational-wave signal models that separate linear amplitudes from nonlinear parameters.
  • If the consistency holds under realistic noise, the method could support near-real-time ringdown parameter estimation during observing runs.
  • Extending the importance-sampling weights to include detector calibration uncertainties would test whether the speedup survives additional systematic parameters.

Load-bearing premise

That analytical marginalization of amplitudes and phases followed by importance sampling produces posteriors and evidences that match full nested sampling or MCMC across the full parameter space and all prior choices.

What would settle it

A side-by-side run of FIREFLY and full nested sampling on the same simulated ringdown signal that contains three or more quasi-normal modes, checking whether the recovered log-evidence and marginal posteriors on the mode frequencies and damping times agree within sampling noise.

read the original abstract

Gravitational-wave (GW) ringdown signals from black holes (BHs) encode crucial information about the gravitational dynamics in the strong-field regime, which offers unique insights into BH properties. In the future, the improving sensitivity of GW detectors is to enable the extraction of multiple quasi-normal modes (QNMs) from ringdown signals. However, incorporating multiple modes drastically enlarges the parameter space, posing computational challenges to data analysis. Inspired by the $F$-statistic method in the continuous GW searches, we develope an algorithm, dubbed as FIREFLY, for accelerating the ringdown signal analysis. FIREFLY analytically marginalizes the amplitude and phase parameters of QNMs to reduce the computational cost and speed up the full-parameter inference from hours to minutes, while achieving consistent posterior and evidence. The acceleration becomes more significant when more QNMs are considered. Rigorously based on the principle of Bayesian inference and importance sampling, our method is statistically interpretable, flexible in prior choice, and compatible with various advanced sampling techniques, providing a new perspective for accelerating future GW data analysis.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The manuscript introduces FIREFLY, an algorithm for Bayesian ringdown analysis of gravitational-wave signals that analytically marginalizes the linear amplitude and phase parameters of multiple quasi-normal modes (QNMs) and employs importance sampling over the nonlinear parameters (frequencies, damping times, etc.). It claims this reduces full-parameter inference from hours to minutes while producing posteriors and evidences consistent with standard nested sampling or MCMC, with the speedup increasing for larger numbers of modes. The method is presented as statistically interpretable, prior-flexible, and compatible with existing samplers, building on F-statistic ideas.

Significance. If the consistency of posteriors and evidences holds under the stated conditions, the approach offers a practical route to multi-mode ringdown inference at scale, which will become essential as detector sensitivity improves and higher-order QNMs become measurable. The analytic marginalization step is a standard consequence of the Gaussian likelihood for stationary noise, and the importance-sampling wrapper preserves the Bayesian interpretation without introducing new fitted quantities.

major comments (2)
  1. [Abstract] Abstract: the central claim that the method 'achieves consistent posterior and evidence' is load-bearing yet unsupported by any reported diagnostics (effective sample size, weight histograms, KL divergence to reference posteriors, or direct comparison tables) across mode counts, prior widths, or frequency separations. Without these, it is impossible to assess whether the importance-sampling approximation remains reliable when modes are closely spaced or priors are broad.
  2. [Abstract] The manuscript does not specify how the proposal distribution for importance sampling is constructed or adapted when the number of QNMs increases; if the proposal is derived from a lower-mode run or from a fixed approximation, the variance of the weights can grow rapidly, undermining the claimed consistency for the full multi-mode case.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive report. The two major comments identify important gaps in the validation and documentation of the importance-sampling step. We address each point below and will revise the manuscript accordingly.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the central claim that the method 'achieves consistent posterior and evidence' is load-bearing yet unsupported by any reported diagnostics (effective sample size, weight histograms, KL divergence to reference posteriors, or direct comparison tables) across mode counts, prior widths, or frequency separations. Without these, it is impossible to assess whether the importance-sampling approximation remains reliable when modes are closely spaced or priors are broad.

    Authors: We agree that the current manuscript does not report the quantitative diagnostics listed. In the revised version we will add (i) effective sample sizes and weight histograms for representative cases, (ii) KL-divergence values between the importance-sampled and nested-sampling posteriors, and (iii) summary tables comparing posterior means, credible intervals and log-evidence differences across mode counts, prior widths and frequency separations. These additions will be placed in a new subsection of the results and referenced from the abstract. revision: yes

  2. Referee: [Abstract] The manuscript does not specify how the proposal distribution for importance sampling is constructed or adapted when the number of QNMs increases; if the proposal is derived from a lower-mode run or from a fixed approximation, the variance of the weights can grow rapidly, undermining the claimed consistency for the full multi-mode case.

    Authors: The referee correctly notes that the construction and adaptation of the proposal distribution are not described. We will revise the methods section to specify (a) how the proposal is obtained (typically from a preliminary run with fewer modes or from an analytic Gaussian approximation to the marginalized likelihood), (b) the procedure used to update the proposal when additional modes are introduced, and (c) any safeguards (e.g., weight truncation or adaptive resampling) employed to control weight variance. We will also add a brief discussion of the scaling of weight variance with mode number. revision: yes

Circularity Check

0 steps flagged

No circularity: standard analytic marginalization + importance sampling

full rationale

The paper derives FIREFLY from the standard Bayesian marginalization of linear amplitude/phase parameters under Gaussian noise (inspired by the external F-statistic) followed by importance sampling over nonlinear parameters. This is a direct application of established inference principles with no self-definitional reduction, no fitted input renamed as prediction, and no load-bearing self-citation chain. The consistency claim is presented as a consequence of the method rather than an input; the derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based on abstract only; the method relies on standard domain assumptions of linear superposition of quasi-normal modes and the validity of importance sampling for recovering the full posterior.

axioms (1)
  • domain assumption Ringdown signals are linear superpositions of quasi-normal modes with independent amplitudes and phases.
    Standard modeling assumption in gravitational-wave ringdown analysis invoked to justify marginalization.

pith-pipeline@v0.9.0 · 5735 in / 1207 out tokens · 50803 ms · 2026-05-23T03:51:16.276752+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Cracking Gravitational Wave Multiple Ringdown Modes in Space

    gr-qc 2026-04 unverdicted novelty 6.0

    FIREFLY algorithm enables 200-fold faster multi-mode ringdown analysis for space-borne gravitational wave detectors while remaining compatible with time-delay interferometry.

  2. FluxMC: Rapid and High-Fidelity Inference for Space-Based Gravitational-Wave Observations

    astro-ph.IM 2026-04 unverdicted novelty 6.0

    FluxMC integrates flow matching with parallel tempering MCMC to converge in under five hours on high-fidelity IMRPhenomHM waveforms for massive black hole binaries, where standard methods fail after hundreds of hours ...

  3. A Robust and Efficient F-statistic-based Framework for Consistent Bayesian Inference of Compact Binary Coalescences

    gr-qc 2025-09 conditional novelty 6.0

    F-statistic framework analytically maximizes over distance and polarization to enable faster Bayesian inference of compact binary coalescences with a new evidence formulation that matches full frequency-domain results...

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages · cited by 3 Pith papers · 23 internal anchors

  1. [1]

    Observation of Gravitational Waves from a Binary Black Hole Merger

    Abbott, B.P., et al. : Observation of Gravitational Waves from a Binary Black Hole Merger. Phys. Rev. Lett. 116(6), 061102 (2016) https://doi.org/10.1103/PhysRevLett.116.061102 arXiv:1602.03837 [gr-qc]

  2. [2]

    Tests of general relativity with GW150914

    Abbott, B.P., et al. : Tests of general relativity with GW150914. Phys. Rev. Lett. 116(22), 221101 (2016) https://doi.org/10.1103/PhysRevLett.116.221101 arXiv:1602.03841 [gr-qc]. [Erratum: Phys.Rev.Lett. 121, 129902 (2018)]

  3. [3]

    Abbott, R., et al.: Tests of General Relativity with GWTC-3 (2021) arXiv:2112.06861 [gr-qc]

  4. [4]

    A gravitational-wave standard siren measurement of the Hubble constant

    Abbott, B.P., et al. : A gravitational-wave standard siren measurement of the Hubble constant. Nature 551(7678), 85–88 (2017) https://doi.org/10.1038/nature24471 arXiv:1710.05835 [astro-ph.CO]

  5. [5]

    GW170817: Measurements of Neutron Star Radii and Equation of State

    Abbott, B.P., et al. : GW170817: Measurements of neutron star radii and equation of state. Phys. Rev. Lett. 121(16), 161101 (2018) https://doi.org/10.1103/PhysRevLett.121.161101 arXiv:1805.11581 [gr-qc]

  6. [6]

    The population of merging compact binaries inferred using gravitational waves through GWTC-3

    Abbott, R., et al. : Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3. Phys. Rev. X 13(1), 011048 (2023) https://doi.org/10.1103/PhysRevX.13.011048 arXiv:2111.03634 [astro-ph.HE]

  7. [7]

    Isi, M., Giesler, M., Farr, W.M., Scheel, M.A., Teukolsky, S.A.: Testing the no-hair theorem with GW150914. Phys. Rev. Lett. 123(11), 111102 (2019) https://doi.org/10.1103/PhysRevLett.123.111102 arXiv:1905.00869 [gr-qc]

  8. [8]

    Isi, M., Farr, W.M., Giesler, M., Scheel, M.A., Teukolsky, S.A.: Testing the Black-Hole Area Law with GW150914. Phys. Rev. Lett. 127(1), 011103 (2021) https://doi.org/10.1103/PhysRevLett.127.011103 arXiv:2012.04486 [gr-qc]

  9. [9]

    Cardoso, V., Franzin, E., Pani, P.: Is the gravitational-wave ringdown a probe of the event horizon? Phys. Rev. Lett. 116(17), 171101 (2016) https://doi.org/10.1103/PhysRevLett.116.171101 arXiv:1602.07309 [gr-qc]. [Erratum: Phys.Rev.Lett. 117, 089902 (2016)]

  10. [10]

    Teukolsky, S.A.: Perturbations of a rotating black hole. 1. Fundamental equations for gravitational elec- tromagnetic and neutrino field perturbations. Astrophys. J. 185, 635–647 (1973) https://doi.org/10.1086/ 152444

  11. [11]

    Berti, E., Cardoso, J., Cardoso, V., Cavaglia, M.: Matched-filtering and parameter estimation of ringdown waveforms. Phys. Rev. D 76, 104044 (2007) https://doi.org/10.1103/PhysRevD.76.104044 arXiv:0707.1202 [gr-qc]

  12. [12]

    Dreyer, O., Kelly, B.J., Krishnan, B., Finn, L.S., Garrison, D., Lopez-Aleman, R.: Black hole spectroscopy: Testing general relativity through gravitational wave observations. Class. Quant. Grav.21, 787–804 (2004) https://doi.org/10.1088/0264-9381/21/4/003 arXiv:gr-qc/0309007

  13. [13]

    Christensen, N., Meyer, R.: Markov chain Monte Carlo methods for Bayesian gravitational radiation data analysis. Phys. Rev. D 58, 082001 (1998) https://doi.org/10.1103/PhysRevD.58.082001

  14. [14]

    Sharma, S.: Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy. Ann. Rev. Astron. Astrophys. 55, 213–259 (2017) https://doi.org/10.1146/annurev-astro-082214-122339 arXiv:1706.01629 [astro-ph.IM]

  15. [15]

    2004, in American Institute of Physics Conference Series, Vol

    Skilling, J.: Nested Sampling. AIP Conf. Proc. 735(1), 395 (2004) https://doi.org/10.1063/1.1835238

  16. [16]

    Bayesian Analysis 1(4), 833–859 (2006) https://doi.org/10.1214/06-BA127

    Skilling, J.: Nested sampling for general Bayesian computation. Bayesian Analysis 1(4), 833–859 (2006) https://doi.org/10.1214/06-BA127

  17. [17]

    Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO

    Reitze, D., et al. : Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO. Bull. Am. Astron. Soc. 51(7), 035 (2019) arXiv:1907.04833 [astro-ph.IM]

  18. [18]

    The US Program in Ground-Based Gravitational Wave Science: Contribution from the LIGO Laboratory

    Reitze, D., et al. : The US Program in Ground-Based Gravitational Wave Science: Contribution from the LIGO Laboratory. Bull. Am. Astron. Soc. 51, 141 (2019) arXiv:1903.04615 [astro-ph.IM] 19

  19. [19]

    : The Einstein Telescope: A third-generation gravitational wave observatory

    Punturo, M., et al. : The Einstein Telescope: A third-generation gravitational wave observatory. Class. Quant. Grav. 27, 194002 (2010) https://doi.org/10.1088/0264-9381/27/19/194002

  20. [20]

    Hild, S., et al.: Sensitivity Studies for Third-Generation Gravitational Wave Observatories. Class. Quant. Grav. 28, 094013 (2011) https://doi.org/10.1088/0264-9381/28/9/094013 arXiv:1012.0908 [gr-qc]

  21. [21]

    Laser Interferometer Space Antenna

    Amaro-Seoane, P., et al. : Laser Interferometer Space Antenna. e-prints (2017) https://doi.org/10.48550/ arXiv.1702.00786 [astro-ph.IM]

  22. [22]

    Hu, W.-R., Wu, Y.-L.: The Taiji Program in Space for gravitational wave physics and the nature of gravity. Natl. Sci. Rev. 4(5), 685–686 (2017) https://doi.org/10.1093/nsr/nwx116

  23. [23]

    TianQin: a space-borne gravitational wave detector

    Luo, J., et al. : TianQin: a space-borne gravitational wave detector. Class. Quant. Grav. 33(3), 035010 (2016) https://doi.org/10.1088/0264-9381/33/3/035010 arXiv:1512.02076 [astro-ph.IM]

  24. [24]

    Bhagwat, S., Pacilio, C., Barausse, E., Pani, P.: Landscape of massive black-hole spectroscopy with LISA and the Einstein Telescope. Phys. Rev. D 105(12), 124063 (2022) https://doi.org/10.1103/PhysRevD. 105.124063 arXiv:2201.00023 [gr-qc]

  25. [25]

    Bhagwat, S., Pacilio, C., Pani, P., Mapelli, M.: Landscape of stellar-mass black-hole spectroscopy with third-generation gravitational-wave detectors. Phys. Rev. D 108(4), 043019 (2023) https://doi.org/10. 1103/PhysRevD.108.043019 arXiv:2304.02283 [gr-qc]

  26. [26]

    Pitte, C., Baghi, Q., Besan¸ con, M., Petiteau, A.: Exploring tests of the no-hair theorem with LISA. Phys. Rev. D 110(10), 104003 (2024) https://doi.org/10.1103/PhysRevD.110.104003 arXiv:2406.14552 [gr-qc]

  27. [27]

    Berti, E., Cardoso, V., Starinets, A.O.: Quasinormal modes of black holes and black branes. Class. Quant. Grav. 26, 163001 (2009) https://doi.org/10.1088/0264-9381/26/16/163001 arXiv:0905.2975 [gr-qc]

  28. [28]

    Jaranowski, P., Krolak, A., Schutz, B.F.: Data analysis of gravitational - wave signals from spinning neutron stars. 1. The Signal and its detection. Phys. Rev. D 58, 063001 (1998) https://doi.org/10.1103/ PhysRevD.58.063001 arXiv:gr-qc/9804014

  29. [29]

    Ashok, A., Covas, P.B., Prix, R., Papa, M.A.: Bayesian F-statistic-based parameter estimation of contin- uous gravitational waves from known pulsars. Phys. Rev. D 109(10), 104002 (2024) https://doi.org/10. 1103/PhysRevD.109.104002 arXiv:2401.17025 [gr-qc]

  30. [30]

    Astrophys

    Wang, H.-T., Yim, G., Chen, X., Shao, L.: Gravitational Wave Ringdown Analysis Using the F -statistic. Astrophys. J. 974(2), 230 (2024) https://doi.org/10.3847/1538-4357/ad7096 arXiv:2409.00970 [gr-qc]

  31. [31]

    Prix, R., Krishnan, B.: Targeted search for continuous gravitational waves: Bayesian versus maximum- likelihood statistics. Class. Quant. Grav. 26, 204013 (2009) https://doi.org/10.1088/0264-9381/26/20/ 204013 arXiv:0907.2569 [gr-qc]

  32. [32]

    LIGO Document, https://dcc.ligo.org/LIGO-T1500618/public

    Prix, R.: Bayesian QNM search on black hole ringdown modes (applied to GW150914). LIGO Document, https://dcc.ligo.org/LIGO-T1500618/public

  33. [33]

    Nature Phys.18(1), 112–117 (2022) https://doi.org/10.1038/s41567-021-01425-7 arXiv:1909.06296 [astro-ph.IM]

    Gabbard, H., Messenger, C., Heng, I.S., Tonolini, F., Murray-Smith, R.: Bayesian parameter estimation using conditional variational autoencoders for gravitational-wave astronomy. Nature Phys.18(1), 112–117 (2022) https://doi.org/10.1038/s41567-021-01425-7 arXiv:1909.06296 [astro-ph.IM]

  34. [34]

    Dax, M., Green, S.R., Gair, J., Macke, J.H., Buonanno, A., Sch¨ olkopf, B.: Real-Time Gravitational Wave Science with Neural Posterior Estimation. Phys. Rev. Lett. 127(24), 241103 (2021) https://doi.org/10. 1103/PhysRevLett.127.241103 arXiv:2106.12594 [gr-qc]

  35. [35]

    Dax, M., Green, S.R., Gair, J., P¨ urrer, M., Wildberger, J., Macke, J.H., Buonanno, A., Sch¨ olkopf, B.: Neu- ral Importance Sampling for Rapid and Reliable Gravitational-Wave Inference. Phys. Rev. Lett. 130(17), 171403 (2023) https://doi.org/10.1103/PhysRevLett.130.171403 arXiv:2210.05686 [gr-qc]

  36. [36]

    Pacilio, C., Bhagwat, S., Cotesta, R.: Simulation-based inference of black hole ringdowns in the time domain. Phys. Rev. D 110(8), 083010 (2024) https://doi.org/10.1103/PhysRevD.110.083010 arXiv:2404.11373 [gr-qc] 20

  37. [37]

    Springer (2004), https://link.springer.com/ book/10.1007/978-1-4757-4145-2

    Robert, C.P., Casella, G.: Monte carlo statistical methods. Springer (2004), https://link.springer.com/ book/10.1007/978-1-4757-4145-2

  38. [38]

    : The SXS Collaboration catalog of binary black hole simulations

    Boyle, M., et al. : The SXS Collaboration catalog of binary black hole simulations. Class. Quant. Grav. 36(19), 195006 (2019) https://doi.org/10.1088/1361-6382/ab34e2 arXiv:1904.04831 [gr-qc]

  39. [39]

    Giesler, M., Isi, M., Scheel, M.A., Teukolsky, S.: Black Hole Ringdown: The Importance of Overtones. Phys. Rev. X 9(4), 041060 (2019) https://doi.org/10.1103/PhysRevX.9.041060 arXiv:1903.08284 [gr-qc]

  40. [40]

    Statistics and Computing 29(5), 891 (2019) https: //doi.org/10.1007/s11222-018-9844-0

    Higson, E., Handley, W., Hobson, M., Lasenby, A.: Dynamic nested sampling: an improved algorithm for parameter estimation and evidence calculation. Statistics and Computing 29(5), 891 (2019) https: //doi.org/10.1007/s11222-018-9844-0

  41. [41]

    Computational Optimal Transport: With Ap- plications to Data Science

    Peyr´ e, G., Cuturi, M.: Computational optimal transport: With applications to data science. Foundations and Trends in Machine Learning 11(5-6), 355–607 (2019) https://doi.org/10.1561/2200000073

  42. [42]

    Wang, Y., Shang, Y., Babak, S.: EMRI data analysis with a phenomenological waveform. Phys. Rev. D 86, 104050 (2012) https://doi.org/10.1103/PhysRevD.86.104050 arXiv:1207.4956 [gr-qc]

  43. [43]

    Abbott, R

    Hu, Q., Veitch, J.: Rapid Premerger Localization of Binary Neutron Stars in Third-generation Gravitational-wave Detectors. Astrophys. J. Lett. 958(2), 43 (2023) https://doi.org/10.3847/2041-8213/ ad0ed4 arXiv:2309.00970 [gr-qc]

  44. [44]

    Berti, E., Klein, A.: Mixing of spherical and spheroidal modes in perturbed Kerr black holes. Phys. Rev. D 90(6), 064012 (2014) https://doi.org/10.1103/PhysRevD.90.064012 arXiv:1408.1860 [gr-qc]

  45. [45]

    Finn, L.S.: Detection, measurement and gravitational radiation. Phys. Rev. D 46, 5236–5249 (1992) https://doi.org/10.1103/PhysRevD.46.5236 arXiv:gr-qc/9209010

  46. [46]

    Frontiers in Astronomy and Space Sciences 11, 1326926 (2024) https://doi.org/10.3389/fspas.2024.1326926 arXiv:2406.18905 [stat.ME]

    Loredo, T.J., Wolpert, R.L.: Bayesian inference: more than Bayes’s theorem. Frontiers in Astronomy and Space Sciences 11, 1326926 (2024) https://doi.org/10.3389/fspas.2024.1326926 arXiv:2406.18905 [stat.ME]

  47. [47]

    Thrane, E., Talbot, C.: An introduction to Bayesian inference in gravitational-wave astronomy: parameter estimation, model selection, and hierarchical models. Publ. Astron. Soc. Austral. 36, 010 (2019) https: //doi.org/10.1017/pasa.2019.2 arXiv:1809.02293 [astro-ph.IM]. [Erratum: Publ.Astron.Soc.Austral. 37, e036 (2020)]

  48. [48]

    WIREs Computational Statistics 2(1), 54–60 (2010) https://doi.org/10.1002/wics.56

    Tokdar, S.T., Kass, R.E.: Importance sampling: a review. WIREs Computational Statistics 2(1), 54–60 (2010) https://doi.org/10.1002/wics.56

  49. [49]

    Bilby: A user-friendly Bayesian inference library for gravitational-wave astronomy

    Ashton, G., et al. : BILBY: A user-friendly Bayesian inference library for gravitational-wave astron- omy. Astrophys. J. Suppl. 241(2), 27 (2019) https://doi.org/10.3847/1538-4365/ab06fc arXiv:1811.02042 [astro-ph.IM] 21