psiphy — Package for Statistical Inference of Physics

A Python package for cosmological parameter inference, providing tools for simulation-based inference (SBI), MCMC sampling, Fisher forecasting, and posterior diagnostics.

Module	Contents
psiphy.mcmc	MHSampler (emcee), NestedSampler (nautilus or dynesty, with Bayesian evidence), importance sampling, sequential_monte_carlo
psiphy.utils.prior	Uniform, log-uniform, Gaussian and mixture priors with one interface for all samplers
psiphy.sbi	Likelihood-free inference: BOLFI, LFIRE, rejection ABC
psiphy.forecasting	FisherMatrix: numerical Fisher forecasts
psiphy.plotting	Triangle plots, credible intervals, posterior diagnostics
psiphy.toy_models	Simulators with known answers for testing

pip install psiphy

New here? Start with Installation and Quickstart, then follow the Tutorials.

Getting started

• Installation
• Quickstart

Tutorials

• Tutorials
  • Toy Models Overview
  • MCMC sampling and priors
  • Nested sampling and Bayesian evidence
  • Importance Sampling and SMC
  • Likelihood-free inference
  • Fisher Matrix Forecasting with psiphy
  • Posterior Diagnostics Tutorial

API Reference

• psiphy.utils
  • kl_divergence()
  • js_distance()
  • Kernel Density Estimation
  • bandwidth_kde_silverman()
  • bandwidth_kde_scott()
  • bandwidth_kdeCV()
  • cross_val_loss_kdeCV_loo()
  • cross_val_loss_kdeCV_kFold()
  • loading_verbose()
  • spherical_to_cartesian()
  • save_model()
  • load_model()
  • LH_sampling()
  • MC_sampling()
  • MCS_nsphere()
  • LHS_nsphere()
  • BasePrior
  • UniformPrior
  • LogUniformPrior
  • GaussianPrior
  • MixturePrior
  • expected_improvement()
  • probability_of_improvement()
  • confidence_bound()
  • get_acquisition()
  • propose_location()
  • check_inside_nsphere()
  • draw_nsphere_surface()
• psiphy.sbi
  • tqdm()
  • rABC
  • emceeSRE
  • LFIRE
  • BOLFI
  • tqdm()
  • ABC_gpL
• psiphy.mcmc
  • tqdm()
  • eff_sample_size()
  • importance_sampling()
  • sequential_importance_sampling()
  • sequential_monte_carlo()
  • SMC()
  • MHSampler
  • NestedSampler
• psiphy.forecasting
  • tqdm
  • FisherMatrix
• psiphy.plotting
  • credible_limit()
  • plot_lfire()
  • plot_1Dmarginal()
  • plot_2Dmarginal()
  • plot_1Dmarginal_lfire()
  • plot_2Dmarginal_lfire()
  • walk_parameter()
  • corner_density()
  • print_CI_samples()
  • density_1D()
  • density_2D()
  • get_CI_HDR_2D()
  • get_CI_HDR_1D()
  • find_intervals()
  • plot_dist_corner_getdist()
  • walk_parameter()
  • plot_corner_dist()
  • plot_corner_multiple_dist()
  • print_CI_samples()
  • density_1D()
  • contour_2D()
  • density_2D()
  • plot_2d_power()
  • get_CI_HDR_2D()
  • get_CI_HDR_1D()
  • find_intervals()
• psiphy.toy_models
  • GaussianSignal
  • gaussian_signal
  • NoisyLine
  • noisy_line
  • MA2Model
  • PowerLawSpectrum
  • LotkaVolterra

Project

• Contributing and support

Indices and tables

• Index

• Module Index