Treatment effects vary across different patients, and estimation of this variability is essential for clinical decisionâmaking. We aimed to develop a model estimating the benefit of alternative treatment options for individual patients, extending a risk modeling approach in a network metaâanalysis framework. We propose a twoâstage prediction model for heterogeneous treatment effects by combining prognosis research and network metaâanalysis methods where individual patient data are available. In the first stage, a prognostic model to predict the baseline risk of the outcome. In the second stage, we use the baseline risk score from the first stage as a single prognostic factor and effect modifier in a network metaâregression model. We apply the approach to a network metaâanalysis of three randomized clinical trials comparing the relapses in Natalizumab, Glatiramer Acetate, and Dimethyl Fumarate, including 3590 patients diagnosed with relapsingâremitting multiple sclerosis. We find that the baseline risk score modifies the relative and absolute treatment effects. Several patient characteristics, such as age and disability status, impact the baseline risk of relapse, which in turn moderates the benefit expected for each of the treatments. For highârisk patients, the treatment that minimizes the risk of relapse in 2âyears is Natalizumab, whereas Dimethyl Fumarate might be a better option for lowârisk patients. Our approach can be easily extended to all outcomes of interest and has the potential to inform a personalized treatment approach.