R code for modelling the biomarker panel