compute_beam_stats function

Compute bootstrap model coefficients for BEAM