<h3>Background</h3> There is expectation of using biomarkers to personalize treatment. Yet, a successful treatment selection cannot be confirmed before 5 or 10 years of progression-free survival (PFS). Treatment individualization based on the probability of an individual patient to achieve undetectable MRD with a singular regimen, could represent a new model towards personalized treatment with fast assessment of its success. This idea has not been investigated previously. <h3>Methods</h3> We sought to define a machine learning model to predict undetectable MRD in newly-diagnosed transplant-eligible MM patients (n=278). The training (n=152) and internal validation cohort (n=60) consisted of 212 active MM patients enrolled in the GEM2012MENOS65 clinical trial. The external validation cohort was defined by 66 high-risk smoldering MM patients enrolled in the GEM-CESAR clinical trial, which treatment differed only by the substitution of bortezomib by carfilzomib during induction and consolidation. Patients were included in the study based on data availability of 17 parameters (p≤.05) associated with MRD outcomes. <h3>Results</h3> We started by investigating patients' MRD status after VRD induction, HDT/ASCT and VRD consolidation (GEM2012MENOS65) according to their ISS and Revised-ISS, LDH levels, and cytogenetic alterations. High LDH levels and del(17p13), two features relatively infrequent at diagnosis, were the only parameters associated with lower rates of undetectable MRD. The ISS and R-ISS were not predictive. Therefore, we aimed to evaluate other disease features associated with MRD outcomes and develop more effective models based on machine learning logistic regression. The most effective one resulted from integrating cytogenetic [t(4;14) and/or del(17p13)], tumor burden (plasma cell clonality in bone marrow and circulating tumor cells in peripheral blood) and immune related (myeloid precursors, mature B cells, intermediate neutrophils, eosinophils, CD27negCD38pos T cells and CD56brightCD27neg NK cells) biomarkers. Data obtained for an individual patient can be substituted into our formula, which results in a numerical probability of achieving undetectable (>0.5) vs persistent (0.685 or <0.365 (observed in 102/212 patients), MRD outcomes are respectively predicted with higher confidence. Standard-confidence, high-confidence, and external validation predictions were accurate in 152/212 (71.7%), 85/102 (83.3%), and 48/66 (72.7%) patients respectively. Patients predicted to achieve undetectable MRD using standard and high-confidence values showed longer PFS and OS than those with probability of persistent MRD. <h3>Conclusions</h3> We demonstrated that selecting a regimen based on probable MRD outcomes, and confirming soon after if that probability was accurate, is a possible new approach towards individualized treatment in MM. The model is available at www.MRDpredictor.com to facilitate its use in clinical practice.