In this talk, an numerical scheme called the linearly backward Milstein method with truncated Wiener process for the stochastic susceptible-infected-susceptible (SIS) epidemic model is proposed. It is proved that the numerical solution provided by this method remains in the initial region (0, N). The convergence and long time behavior of the numerical solution are analyzed. It is shown that the method not only has first-order mean-square convergence, but also maintains the extinction and persistence of the underlying SIS model. Finally, some numerical experiments will be provided to validate the theoretical results.