Unreinforced masonry (URM) walls are commonly found in historic and legacy buildings around the world. The structural resistance of URM walls under in-plane (IP) and out-of-plane (OOP) loads is of primary concern to engineers, as their failure is generally sudden, with catastrophic loss of strength and structural integrity. Due to the complex behavior and inherent uncertainties of the masonry material, engineers opt for the use of low-fidelity (LF) resistance models with limited accuracy, such as design code models and other simplified analytical models in the literature. Models with enhanced prediction accuracy have attracted growing attention, particularly when uncertainty analysis (e.g., reliability evaluation) is needed. As such, high-fidelity (HF) models, such as nonlinear finite element models based on advanced computational mechanics, have been developed and used to characterize the structural behaviors and failure modes of URM walls, particularly the resistance, with remarkable success in terms of accuracy. However, their direct use for resistance prediction and uncertainty analysis is scarce due to the computational burden and technical complications. To address this issue, this study takes an efficient multifidelity (MF) approach that leverages both HF and LF models via information fusion to enhance LF models with only a few HF model evaluations for URM walls. The main research thrust is to develop an MF surrogate model to facilitate uncertainty analysis in the IP and OOP resistance of URM walls. The analysis results indicate that the MF surrogate models developed are capable of achieving significant improvements in terms of accuracy and efficiency in predicting the IP and OOP resistances of URM walls both deterministically and probabilistically, compared with the LF model and the surrogate model developed only based on a limited number of HF model runs.