Knowing the maintainability of a component or a system means that repair resource allocations, such as spare part procurement and maintenance training, can be planned and optimized more effectively. Repair data are often collected from multiple and distributed units in different operational conditions, which can introduce heterogeneity into the data. Part of such heterogeneity can be explained and isolated by the observable covariates, whose values and the way that they can affect the item’s maintainability are known. However, some factors which may affect maintainability are typically unknown (unobserved covariates), leading to unobserved heterogeneity. Nevertheless, many researchers have ignored the effect of observed and un-observed covariates, and this may lead to erroneous model selection, as well as wrong conclusions and decisions. Moreover, many authors have simplified their analysis by considering a complex system as a single item. In these studies, the assumption is that all repair data represent an identical repair process for the item. In practice, mechanical systems are composed of multiple parts, with various failure mechanisms, which need different repair processes (repair modes) to return to the operational phase; classical distribution, such as lognormal, which is only a function of time, may not be able to model such complexity. The paper utilizes the mixture frailty model (MFM) in the presence of some specific observed or unobserved covariates to predict maintainability more precisely. MFMs can model the effect of observed and unobserved covariates, as well as identifying different repair processes in the repair dataset. The application of the proposed model is demonstrated by a case study.