Modeling the Impact of Vaccination and Post-Treatment on Rabies Transmission
Article Information
Abstract
Rabies remains a serious public health concern, as dog bites account for the majority of human cases. In this study, we develop a comprehensive mathematical model to investigate the dynamics of rabies transmission by incorporating two key intervention strategies: an asymptotic class (\(P\)) and a booster vaccination class (\(B\)). The basic reproduction number (\(R_0\)) is derived as a threshold parameter that governs whether the disease spreads or dies out, based on a system of nonlinear differential equations. A sensitivity analysis of \(R_0\) is conducted to identify the most influential parameters affecting disease transmission. The results indicate that the transmission rate (\(\beta\)) has a significant impact on \(R_0\), emphasizing the importance of reducing contact between susceptible and infected populations. Stability analysis of both the disease-free and endemic equilibria reveals that the disease can be eradicated when \(R_0 < 1\), whereas it persists when \(R_0 > 1\). Numerical simulations, performed using the classical Runge–Kutta method, illustrate the comparative effects of vaccination and treatment interventions. The results demonstrate that the inclusion of the booster vaccination class (\(B\)) substantially reduces the number of infected individuals compared to the scenario without vaccination. Moreover, the therapeutic class (\(P\)) further accelerates recovery and alleviates the overall disease burden. Analysis of varying transmission rates (\(\beta\)) shows that higher values lead to a rapid rise in infection levels, underscoring the necessity for effective intervention strategies in high-contact environments. Overall, the model highlights that integrating booster vaccination and treatment measures within rabies control programs can significantly decrease disease prevalence. The findings of this study provide valuable insights for public health authorities aiming to design efficient strategies for the control and eventual eradication of rabies.
Graphical Abstract
Keywords
Data Availability Statement
Funding
Conflicts of Interest
Ethical Approval and Consent to Participate
References
- Panjeti, V. G., & Real, L. A. (2011). Mathematical models for rabies. Advances in Virus Research, 79, 377–395.
[CrossRef] [Google Scholar] - Layan, M., Dellicour, S., Baele, G., Cauchemez, S., & Bourhy, H. (2021). Mathematical modelling and phylodynamics for the study of dog rabies dynamics and control: A scoping review. PLOS Neglected Tropical Diseases, 15(5), e0009449.
[CrossRef] [Google Scholar] - Hankins, D. G., & Rosekrans, J. A. (2004). Overview, prevention, and treatment of rabies. Mayo Clinic Proceedings, 79(5), 671–676.
[CrossRef] [Google Scholar] - Demirci, E. (2014). A new mathematical approach for rabies endemy. Applied Mathematical Sciences, 8(2), 59–67. http://dx.doi.org/10.12988/ams.2014.39525
[Google Scholar] - Tulu, A. M. (2017). Mathematical Modeling on the Impact of Infective Immigrants on the Spread of Dog Rabies (Doctoral dissertation).
[Google Scholar] - Barecha, C. B., Girzaw, F., Kandi, R. V., & Pal, M. (2017). Epidemiology and public health significance of rabies. Perspectives in Clinical Research, 5(1), 55–67.
[Google Scholar] - Mumbu, A. J. (2022). Modelling dynamics of dog rabies disease with vaccination and treatment in dog population [Doctoral dissertation, University of Dodoma].
[Google Scholar] - Hampson, K., Coudeville, L., Lembo, T., Sambo, M., Kieffer, A., Attlan, M., \ldots (2015). Estimating the global burden of endemic canine rabies. PLOS Neglected Tropical Diseases, 9(4), e0003709.
[CrossRef] [Google Scholar] - Changalucha, J., Steenson, R., Grieve, E., Cleaveland, S., Lembo, T., Lushasi, K., Mchau, G., Mtema, Z., Sambo, M., Nanai, A., \ldots (2019). The need to improve access to rabies post-exposure vaccines: Lessons from Tanzania. Vaccine, 37(Suppl. 1), A45–A53.
[CrossRef] [Google Scholar] - Warrell, D. A. (2007). Rabies and envenomings: a neglected public health issue: report of a consultative meeting, World Health Organization, Geneva, 10 January 2007. World Health Organization.
[Google Scholar] - WHO Rabies Modelling Consortium. (2020). Zero human deaths from dog-mediated rabies by 2030: perspectives from quantitative and mathematical modelling. Gates Open Research, 3, 1564.
[CrossRef] [Google Scholar] - Fagbo, S., & Asiri, A. (2022). Sand cat rabies in Saudi Arabia: Leveraging to improve spillover surveillance and one health impact. International Journal of Infectious Diseases, 116, S106-S107.
[CrossRef] [Google Scholar] - Ruan, S. (2017). Spatiotemporal epidemic models for rabies among animals. Infectious Disease Modelling, 2(3), 277–287.
[CrossRef] [Google Scholar] - Kanankege, K. S., Errecaborde, K. M., Wiratsudakul, A., Wongnak, P., Yoopatthanawong, C., Thanapongtharm, W., ... & Perez, A. (2022). Identifying high-risk areas for dog-mediated rabies using Bayesian spatial regression. One Health, 15, 100411.
[CrossRef] [Google Scholar] - Huang, Y., & Li, M. (2020). Application of a mathematical model in determining the spread of the rabies virus: Simulation study. JMIR Medical Informatics, 8(5), e18627.
[CrossRef] [Google Scholar] - Coyne, M. J., Smith, G., & McAllister, F. E. (1989). Mathematic model for the population biology of rabies in raccoons in the mid-Atlantic states. American Journal of Veterinary Research, 50(12), 2148–2154.
[Google Scholar] - Eze, O. C., Mbah, G. E., Nnaji, D. U., & Onyiaji, N. E. (2020). Mathematical modelling of transmission dynamics of rabies virus. International Journal of Mathematics Trends and Technology, 66(6), 126–134.
[CrossRef] [Google Scholar] - Agaba, G. O. (2014). Modelling of the spread of rabies with pre-exposure vaccination of humans. Mathematical Theory and Modeling, 4(8), 57–63.
[Google Scholar] - Ruan, S. (2017). Modeling the transmission dynamics and control of rabies in China. Mathematical Biosciences, 286, 65–93.
[CrossRef] [Google Scholar] - Ndendya, J. Z., Leandry, L., & Kipingu, A. M. (2023). A next-generation matrix approach using Routh–Hurwitz criterion and quadratic Lyapunov function for modeling animal rabies with infective immigrants. Healthcare Analytics, 4, 100260.
[CrossRef] [Google Scholar] - Omame, A., Isah, M. E., & Abbas, M. (2023). An optimal control model for COVID-19, zika, dengue, and chikungunya co-dynamics with reinfection. Optimal Control Applications and Methods, 44(1), 170–204.
[CrossRef] [Google Scholar] - Agbata, C. (2025). A Mathematical Model for the Control of Ebola Virus Disease with Vaccination Effect. Journal of Mathematical Epidemiology, 1(1), 26-48.
[CrossRef] [Google Scholar] - Liu, J., Jia, Y., & Zhang, T. (2017). Analysis of a rabies transmission model with population dispersal. Nonlinear Analysis: Real World Applications, 35, 229–249.
[CrossRef] [Google Scholar] - Gonzalez-Roldán, J. F., Undurraga, E. A., Meltzer, M. I., Atkins, C., Vargas-Pino, F., Gutiérrez-Cedillo, V., & Hernández-Pérez, J. R. (2021). Cost-effectiveness of the national dog rabies prevention and control program in Mexico, 1990–2015. PLOS Neglected Tropical Diseases, 15(3), e0009130.
[CrossRef] [Google Scholar] - Kanda, K., Jayasinghe, A., Jayasinghe, C., & Yoshida, T. (2022). A regional analysis of the progress of current dog-mediated rabies control and prevention. Pathogens, 11(10), 1130.
[CrossRef] [Google Scholar] - Lushasi, K., Hayes, S., Ferguson, E. A., Changalucha, J., Cleaveland, S., Govella, N. J., & Hampson, K. (2021). Reservoir dynamics of rabies in south-east Tanzania and the roles of cross-species transmission and domestic dog vaccination. Journal of Applied Ecology, 58(11), 2673–2685.
[CrossRef] [Google Scholar] - Bornaa, C. S., Seidu, B., & Daabo, M. I. (2020). Mathematical analysis of rabies infection. Journal of Applied Mathematics, 2020(1), 1804270.
[CrossRef] [Google Scholar] - Ega, T. T., Luboobi, L. S., & Kuznetsov, D. (2015). Modeling the dynamics of rabies transmission with vaccination and stability analysis. Journal of Applied Mathematics, 2015, 1–12.
[Google Scholar] - Ndii, M. Z., Amarti, Z., Wiraningsih, E. D., & Supriatna, A. K. (2018). Rabies epidemic model with uncertainty in parameters: Crisp and fuzzy approaches. IOP Conference Series: Materials Science and Engineering, 332(1), 012031.
[CrossRef] [Google Scholar] - Jehangir, H., Ali, N., Ahmad, I., Younas, H., Ahmad, Z., & El-Agamy, R. F. (2024). Global Dynamics and Numerical Simulation of a Vaccinated Mathematical Model for Rabies Disease. Global Journal of Sciences, 1(1), 14-27.
[CrossRef] [Google Scholar] - Younas, H., Ahmad, I., Ali, N., Haq, I. U., Albalwi, M. D., Muhammad, S., & Shuaib, M. (2024). Modeling rabies dynamics: the impact of vaccination and infectious immigrants on public health. Contemp. Math, 5(3), 3255.
[CrossRef] [Google Scholar]
Cite This Article
TY - JOUR AU - Ullah, Hanif AU - Ali, Amjid AU - Younas, Hazrat PY - 2025 DA - 2025/12/11 TI - Modeling the Impact of Vaccination and Post-Treatment on Rabies Transmission JO - Journal of Numerical Simulations in Physics and Mathematics T2 - Journal of Numerical Simulations in Physics and Mathematics JF - Journal of Numerical Simulations in Physics and Mathematics VL - 1 IS - 2 SP - 84 EP - 97 DO - 10.62762/JNSPM.2025.826101 UR - https://www.icck.org/article/abs/JNSPM.2025.826101 KW - rabies KW - mathematical model KW - stability KW - bifurcation KW - rabies model AB - Rabies remains a serious public health concern, as dog bites account for the majority of human cases. In this study, we develop a comprehensive mathematical model to investigate the dynamics of rabies transmission by incorporating two key intervention strategies: an asymptotic class (\(P\)) and a booster vaccination class (\(B\)). The basic reproduction number (\(R_0\)) is derived as a threshold parameter that governs whether the disease spreads or dies out, based on a system of nonlinear differential equations. A sensitivity analysis of \(R_0\) is conducted to identify the most influential parameters affecting disease transmission. The results indicate that the transmission rate (\(\beta\)) has a significant impact on \(R_0\), emphasizing the importance of reducing contact between susceptible and infected populations. Stability analysis of both the disease-free and endemic equilibria reveals that the disease can be eradicated when \(R_0 < 1\), whereas it persists when \(R_0 > 1\). Numerical simulations, performed using the classical Runge–Kutta method, illustrate the comparative effects of vaccination and treatment interventions. The results demonstrate that the inclusion of the booster vaccination class (\(B\)) substantially reduces the number of infected individuals compared to the scenario without vaccination. Moreover, the therapeutic class (\(P\)) further accelerates recovery and alleviates the overall disease burden. Analysis of varying transmission rates (\(\beta\)) shows that higher values lead to a rapid rise in infection levels, underscoring the necessity for effective intervention strategies in high-contact environments. Overall, the model highlights that integrating booster vaccination and treatment measures within rabies control programs can significantly decrease disease prevalence. The findings of this study provide valuable insights for public health authorities aiming to design efficient strategies for the control and eventual eradication of rabies. SN - 3068-9082 PB - Institute of Central Computation and Knowledge LA - English ER -
@article{Ullah2025Modeling,
author = {Hanif Ullah and Amjid Ali and Hazrat Younas},
title = {Modeling the Impact of Vaccination and Post-Treatment on Rabies Transmission},
journal = {Journal of Numerical Simulations in Physics and Mathematics},
year = {2025},
volume = {1},
number = {2},
pages = {84-97},
doi = {10.62762/JNSPM.2025.826101},
url = {https://www.icck.org/article/abs/JNSPM.2025.826101},
abstract = {Rabies remains a serious public health concern, as dog bites account for the majority of human cases. In this study, we develop a comprehensive mathematical model to investigate the dynamics of rabies transmission by incorporating two key intervention strategies: an asymptotic class (\(P\)) and a booster vaccination class (\(B\)). The basic reproduction number (\(R\_0\)) is derived as a threshold parameter that governs whether the disease spreads or dies out, based on a system of nonlinear differential equations. A sensitivity analysis of \(R\_0\) is conducted to identify the most influential parameters affecting disease transmission. The results indicate that the transmission rate (\(\beta\)) has a significant impact on \(R\_0\), emphasizing the importance of reducing contact between susceptible and infected populations. Stability analysis of both the disease-free and endemic equilibria reveals that the disease can be eradicated when \(R\_0 < 1\), whereas it persists when \(R\_0 > 1\). Numerical simulations, performed using the classical Runge–Kutta method, illustrate the comparative effects of vaccination and treatment interventions. The results demonstrate that the inclusion of the booster vaccination class (\(B\)) substantially reduces the number of infected individuals compared to the scenario without vaccination. Moreover, the therapeutic class (\(P\)) further accelerates recovery and alleviates the overall disease burden. Analysis of varying transmission rates (\(\beta\)) shows that higher values lead to a rapid rise in infection levels, underscoring the necessity for effective intervention strategies in high-contact environments. Overall, the model highlights that integrating booster vaccination and treatment measures within rabies control programs can significantly decrease disease prevalence. The findings of this study provide valuable insights for public health authorities aiming to design efficient strategies for the control and eventual eradication of rabies.},
keywords = {rabies, mathematical model, stability, bifurcation, rabies model},
issn = {3068-9082},
publisher = {Institute of Central Computation and Knowledge}
}
Article Metrics
Publisher's Note
ICCK stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and Permissions
Copyright © 2025 by the Author(s). Published by Institute of Central Computation and Knowledge. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
Portico