“…Accounting for the complexity of immune responses to a tumor, various mathematical models were developed to understand its dynamical mechanism in the past few decades. Among these models, most of them are deterministic (e.g., [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18] and references therein), and only a few of them are stochastic ( e.g., [19,20,21,22,23] and references therein). In the investigation of deterministic models, some researchers focus on the estimation of parameters and the analysis of dynamical properties, such as stability, bifurcation and its stability and direction [4,5,6,7,8,9,10,11], while others further extend the applications of these models by exploring the optimal control of cancer treatment [12,13,14,15,16,17,18].…”