beta_estimation (pyleoclim.utils.wavelet.beta_estimation)

pyleoclim.utils.wavelet.beta_estimation(psd, freq, fmin=None, fmax=None, logf_binning_step='max', verbose=False)[source]

Estimate the power slope of a 1/f^beta process.

Parameters

  • psd (array) – the power spectral density

  • freq (array) – the frequency vector

  • fmin (float) – the min of frequency range for beta estimation

  • fmax (float) – the max of frequency range for beta estimation

  • verbose (bool) – if True, will print out debug information

Returns

  • beta (float) – the estimated slope

  • f_binned (array) – binned frequency vector

  • psd_binned (array) – binned power spectral density

  • Y_reg (array) – prediction based on linear regression