Neutron reflectometry is the foremost technique for in situ determination of the volume fraction profiles of polymer brushes at planar interfaces. However, the subtle features in the reflectometry data produced by these diffuse interfaces challenge data interpretation. Historically, data analyses have used least-squares approaches that do not adequately quantify the uncertainty of the modeled profile and ignore the possibility of other structures that also match the collected data (multimodality). Here, a Bayesian statistical approach is used that permits the structural uncertainty and multimodality to be quantified for polymer brush systems. A free-form model is used to describe the volume fraction profile, minimizing assumptions regarding brush structure, while only allowing physically reasonable profiles to be produced. The model allows the total volume of polymer and the profile monotonicity to be constrained. The rigor of the approach is demonstrated via a round-trip analysis of a simulated system, before it is applied to real data examining the well characterized collapse of a thermoresponsive brush. It is shown that, while failure to constrain the interfacial volume and consider multimodality may result in erroneous structures being derived, carefully constraining the model allows for robust determination of polymer brush compositional profiles. This work highlights that an appropriate combination of flexibility and constraint must be used with polymer brush systems to ensure the veracity of the analysis. The code used in this analysis is provided, enabling the reproduction of the results and the application of the method to similar problems.