In this paper, based on Euler–Marryama method and theory of stochastic processes, a stochastic discrete SIVS epidemic model with general nonlinear incidence and vaccination is proposed by adding random perturbation and then discretizing the corresponding stochastic differential equation model. Firstly, the basic properties of continuous and discrete deterministic SIVS epidemic models are obtained. Then a criterion on the asymptotic mean-square stability of zero solution for a general linear stochastic difference system is established. As the applications of this criterion, the sufficient conditions on the stability in probability of the disease-free and endemic equilibria for the stochastic discrete SIVS epidemic model are obtained. The numerical simulations are given to illustrate the theoretical results.