ObjectivesMulti-parametric MRI techniques, including intravoxel incoherence motion (IVIM), iterative decomposition of water and fat with echo asymmetry and least squares estimation (IDEAL IQ), T2* mapping and T2 mapping, were employed to develop and validate a predictive model for non-alcoholic steatohepatitis (NASH) diagnosis and liver fibrosis (LF) staging in rats. The combined model was interpreted using SHapley Additive exPlanations (SHAP) values for model understanding. Materials and methods160 healthy Sprague-Dawley (SD) rats were divided into control (n = 24) and experimental (n = 136) groups, and the 12-week and 16-week groups were injected intraperitoneally with carbon tetrachloride (CCl4) for 4 weeks, one month before the final feeding period. All rats were subjected to pathological examination to determine LF stage. Upon the study’s completion, 147 SD rats were assessed for liver fibrosis. Results84 SD rats were diagnosed with NASH and 31, 10, and 43 rats were histologically diagnosed with no fibrosis (F0), early LF (F1-F2), and late LF (F3-F4). For diagnosis of NASH and staging of liver fibrosis associated with NASH, a combined multi-parameter magnetic resonance imaging prediction model has a higher area under the ROC curve (AUC) than individual diffusion parameters, especially in advanced stages of fibrosis, with an AUC of 0.929 for the combined model. In SHAP, the free fluid(FF) value contributes most to the model for diagnosing NASH and advanced liver fibrosis, while the T2 value contributes most for diagnosing liver fibrosis and the apparent diffusion coefficient (ADC) value contributes most for diagnosing liver cirrhosis. ConclusionsThe integration of several magnetic resonance imaging methods along with SHAP value analysis provides a more comprehensive and theoretical investigation into the mechanisms and factors contributing to the progression of liver fibrosis in NASH, offering insights from various angles.