{ "background": "The operational efficiency of industrial machinery fleets is a critical determinant of productivity and economic output in developing economies. Current diagnostic methods often rely on aggregated metrics, failing to account for heterogeneous machine types, operational conditions, and site-specific factors, leading to imprecise efficiency estimates and suboptimal maintenance strategies.", "purpose and objectives": "This working paper develops and evaluates a novel Bayesian hierarchical modelling framework for the efficiency diagnostics of heterogeneous machinery fleets. The primary objective is to provide a robust, probabilistic method for quantifying efficiency gains at multiple operational levels, from individual units to entire fleets, while formally incorporating uncertainty.", "methodology": "We propose a three-level hierarchical model where the operational output $y{ij} \\sim \\text{Normal}(\\mu{ij}, \\sigma^2)$ for machine $i$ at site $j$ is modelled with a mean $\\mu{ij} = \\alpha + \\betaj + \\gammai + \\mathbf{X}{ij}^T\\boldsymbol{\\delta}$. Here, $\\betaj \\sim \\text{Normal}(0, \\tau{\\text{site}}^2)$ and $\\gammai \\sim \\text{Normal}(0, \\tau{\\text{machine}}^2)$ are random effects for site and machine type, respectively. Parameters are estimated using Hamiltonian Monte Carlo, with posterior credible intervals used for inference on efficiency parameters.", "findings": "Application to a case-study fleet demonstrates the model's utility in identifying latent inefficiencies. The posterior distribution for the site-level variance parameter $\\tau_{\\text{site}}^2$ indicated that approximately 40% of the variation in output was attributable to location-specific operational practices, a factor obscured in conventional analyses. Credible intervals for machine-type effects revealed that certain categories were operating at significantly lower efficiency than fleet-wide averages.", "conclusion": "The Bayesian hierarchical model offers a superior diagnostic tool for fleet management by providing granular, probabilistic efficiency estimates that explicitly model variability across operational