Returns annual and hourly wages at the mean and the 10th, 25th, median,
75th, and 90th percentiles (25th/75th are the usual junior/senior
anchors), plus estimated burdened hourly rates using the specified
burden multiplier range. Hourly figures are BLS's published hourly
wages (annual / 2080 only where BLS publishes no hourly wage). Also
returns the cell's employment and relative standard errors
(reliability) and flags thin estimates.

Wages are point-in-time estimates for the OEWS reference month
(wage_period in the response), not current rates: escalate them to the
period of performance before pricing.

BLS wages are BASE wages (no fringe, overhead, G&A, or profit).
Multiply by a burden factor to estimate fully-loaded rates:
- 1.5x-1.7x: lean contractor
- 1.8x-2.2x: mid-range professional services (default)
- 2.0x-2.5x: large contractor with clearance overhead
- 2.5x-3.0x: high-overhead (SCIF, deployed)

The burdened range should roughly align with GSA CALC+ ceiling rates
for comparable labor categories. If CALC+ >> burdened BLS, the role
may require specialized skills or clearance overhead. Document the gap.

Burden multipliers must be positive and burden_low <= burden_high.
Reasonable range: 1.3 (lean) to 4.0 (high-overhead/clearance).