I fitted two Lorentzian peaks in each spectrum, using a local linear baseline around each peak. The fitted amplitude is the peak intensity; \(x_c\) is the peak position.

The allometric fits use \(y=aP^b\). For the dominant peak, \(b=0.913\pm0.014\) for intensity and \(b=0.00457\pm0.00010\) for position. For the second peak, the corresponding exponents are \(0.853\pm0.067\) and \(0.00181\pm0.00056\).

Generated [the two power-dependence plots](power_dependence.png) and [the per-power peak fit results](lorentzian_peak_fits.csv).

RESULT: b = 0.913 ± 0.014