Count data often exhibit discrepancies in the frequencies of zeros, which commonly occur across various application domains. These data may include excess zeros (zero inflation) or, less frequently, a scarcity of zeros (zero deflation). In regression models, both situations can arise at different levels of covariates. The zero-modified power series regression model provides an effective framework for modeling such count data, as it does not require prior knowledge of the type of zero modification, whether zero inflation or zero deflation, and can accommodate overdispersion, equidispersion, or underdispersion present in the data. This paper proposes a Bayesian estimation procedure based on the stochastic gradient Hamiltonian Monte Carlo algorithm, effectively addressing many challenges associated with estimating the model parameters. Additionally, we introduce a measure of Bayesian efficiency to evaluate the impact of prior information on parameter estimation. The practical utility of the proposed method is demonstrated through both simulated and real data across different types of zero modification.