I fit two Lorentzian peaks for each of the 15 power spectra. The allometric model was \(y=aP^b\), fit by linear regression in log space; uncertainties are standard errors on \(b\).

- Peak 1 intensity: \(b=0.9233\pm0.0121\)
- Peak 2 intensity: \(b=0.6835\pm0.0862\)
- Peak 1 position \(x_c\): \(b=0.00447\pm0.00008\)
- Peak 2 position \(x_c\): \(b=0.00229\pm0.00039\)

Files: [peak intensity plot](peak_intensity_vs_power.png), [peak position plot](peak_position_vs_power.png), and [fitted peak values for each power](lorentzian_peak_fits.csv).

RESULT: b = 0.9233 ± 0.0121