Pioran is a Julia package to estimate the power spectrum of any time series using bending power laws. This method is based on Gaussian process regression with the fast algorithm of Foreman-Mackey, et al. 2017 and by approximating the power spectrum model with basis functions. The method is described in Lefkir, et al. 2023 (in prep).
using Pkg; Pkg.add("Pioran")See the documentation at https://www.mehdylefkir.fr/Pioran.jl.
Example scripts are provided in the examples directory. To infer the parameters of the power spectrum, I use either Turing.jl for Hamiltonian Monte Carlo or the Python library ultranest for nested sampling. The scripts are written in a way that you can use either of these libraries.