The values of [ypred]j are calculated using the routine ‘ode15s’ in the Matlab and the values of [yexp]j are estimated based on clinical observations (specifically, MM-cell density is estimated to increase 5-fold, IL-6 concentration to increase 10-fold, RANKL concentration to increase 4-fold and OPG concentration to decrease 0.7-fold, at the end of the simulations (see Table 6)). Here, TNFRSF11B is linked to Miyoshi myopathy.