Then, by evaluating V −1 , the control reproduction number is defined by the spectral radius of F V −1 and given by R c = max{R s c , R r 0 }, where R s c and R r 0 are expressed in (1) and (2) of the main text, respectively. 