PBjam analyzes the oscillation spectra of solar-like oscillators. The code performs two main tasks: identifying a set of modes of interest in a spectrum of oscillations, and accurately modeling those modes to measure their frequencies. Mode identification relies on a large set of previous observations of the model parameters, which are then used to construct a prior distribution to inform the sampling. PJjam models the modes using a nested sampling or MCMC algorithm, where Lorentzian profiles are fit to each of the identified modes.