In this paper we present our study on the critical behavior of a stochastic anisotropic Bak-Sneppen (saBS) model, in which a parameter α is introduced to describe the interaction strength among nearest species. We estimate the threshold fitness fc and the critical exponent τr by numerically integrating a master equation for the distribution of avalanche spatial sizes. Other critical exponents are then evaluated from previously known scaling relations. The numerical results are in good agreement with the counterparts yielded by the Monte Carlo simulations. Our results indicate that all saBS models with nonzero interaction strength exhibit self-organized criticality, and fall into the same universality class, by sharing the universal critical exponents.PACS. 05.65.+b Criticality, self-organized -89.75.Fb Self-organization, complex systems -05.10.Gg Stochastic models in statistical physics and nonlinear dynamics