An inverse Mean-Field Homogenization (MFH) process is developed to improve the computational efficiency of non-linear stochastic multiscale analyzes by relying on a micro-mechanics model. First full-field simulations of composite Stochastic Volume Element (SVE) realizations are performed to characterize the homogenized stochastic behavior. The uncertainties observed in the non-linear homogenized response, which result from the uncertainties of their micro-structures, are then translated to an incrementalsecant MFH formulation by defining the MFH input parameters as random effective properties. These effective input parameters, which correspond to the micro-structure geometrical information and to the material phases model parameters, are identified by conducting an inverse analysis from the full-field homogenized responses. Compared to the direct finite element analyzes on SVEs, the resulting stochastic MFH process reduces not only the computational cost, but also the order of uncertain parameters in the composite micro-structures, leading to a stochastic Mean-Field Reduced Order Model (MF-ROM). A data-driven stochastic model is then built in order to generate the random effective properties under the form of a random field used