A MATHEMATICAL MODEL OF DIPHTHERIA TRANSMISSION DYNAMICS WITH HETEROGENEOUS SUSCEPTIBILITY

Keywords: Diphtheria disease, Mathematical model, Stability analysis, Parameter estimation, Vaccination

Abstract

Despite the availability of vaccines, diphtheria continues to pose a public health risk in Indonesia due to uneven vaccination coverage across regions. Previous models have not distinguished between highly susceptible (unvaccinated) and susceptible (vaccinated) populations, nor have they been calibrated with actual Indonesian epidemiological data. To address this gap, this study develops a five-compartment diphtheria transmission model: Highly susceptible (unvaccinated)-Susceptible (vaccinated)-Exposed-Infectious-Recovered (S_1 S_2 EIR), which incorporates two levels of susceptibility based on vaccination status, using empirical diphtheria case data in Indonesia from 2012 to 2023. The analysis begins by proving the positivity, boundedness, and uniqueness of solutions, followed by the calculation of the basic reproduction number using the Next-Generation Matrix method. The analysis shows that the disease-free equilibrium (DFE) is locally and globally asymptotically stable when R₀<1, while the endemic equilibrium (EE) is globally stable when R₀>1. Simulations indicate that the interaction parameter for the unvaccinated group  η₁, strongly accelerates epidemic growth, leading to a higher and earlier infection peak, whereas increased vaccination coverage and recovery rates effectively suppress transmission. This model can be used because the Mean Absolute Percentage Error (MAPE) between the data and the model solution for diphtheria cases in Indonesia is 8.77%.  These results highlight the importance of interventions focused on highly susceptible groups to prevent more severe outbreaks. Therefore, this study is significant in strengthening the theoretical understanding of diphtheria transmission, while also providing data-driven insights as recommendations for policymakers to implement effective and efficient outbreak control measures.

Downloads

Download data is not yet available.

References

WHO, “DIPHTHERIA.” Accessed: Jun. 10, 2024. [Online]. Available: https://www.who.int/news-room/questions-and-answers/item/diphtheria

Ministry of Health of the Republic of Indonesia, “GUIDELINES FOR PREVENTION AND CONTROL OF DIPHTHERIA.” Accessed: Jun. 11, 2024. [Online]. Available: https://www.cdc.gov/diphtheria/about/?CDC_AAref_Val

Anderson, “WHAT TO KNOW ABOUT DIPHTHERIA.” Accessed: Jun. 10, 2024. [Online]. Available: https://www.webmd.com/a-to-z-guides/what-to-know-diphtheria-causes.

National Health Service, “DIPHTHERIA.” Accessed: Jun. 10, 2024. [Online]. Available: https://www.nhs.uk/conditions/diphtheria/.

Centers for Disease Control and Prevention, “DIPHTHERIA.” [Online]. Available: https://www.cdc.gov/diphtheria/about/?CDC_AAref_Val=https://www.cdc.gov/diphtheria/about/causes-transmission.html

K. Rogers, “DIPHTHERIA,” The Editors of Encyclopaedia Britannica. Accessed: Jul. 10, 2024. [Online]. Available: https://www.britannica.com/science/diphtheria.

S. R. Lamichhane, “DIPHTHERIA,” StatPearls. Accessed: Sep. 15, 2025. [Online]. Available: https://www.ncbi.nlm.nih.gov/books/NBK560911/.

Ministry of Health of the Republic of Indonesia, “INDONESIAN HEALTH PROFILE 2023.” Accessed: Jun. 12, 2024. [Online]. Available: https://kemkes.go.id/id/category-download/profil-kesehatan

T. Witelski and M. Bowen, “METHODS OF MATHEMATICAL MODELLING: CONTINUOUS SYSTEMS AND DIFFERENTIAL EQUATIONS,” Contin. Syst. Differ. Equations, pp. 23–45, 2015. doi : https://doi.org/10.1007/978-3-319-23042-9_2

F. Brauer, C. Castillo-Chavez, F. Brauer, and C. Castillo-Chavez, “EPIDEMIC MODELS,” Math. Model. Popul. Biol. Epidemiol., pp. 345–409, 2008.doi : https://doi.org/10.1007/978-1-4614-1686-9_9

W. Kermack and A. McKendrick, “A CONTRIBUTION TO THE MATHEMATICAL THEORY OF EPIDEMICS,” in Proceedings of the Royal Society of London 115, 1927. doi : https://doi.org/10.1098/rspa.1927.0118

H. Hethcote, “THREE BASIC EPIDEMIOLOGICAL MODELS,” Appl. Math. Ecol., pp. 119–144, 1989. doi : https://doi.org/10.1007/978-3-642-61317-3_5

G. Puspita, M. Kharis, and Supriyono, “PEMODELAN MATEMATIKA PADA PENYEBARAN PENYAKIT DIFTERI DENGAN PENGARUH KARANTINA DAN VAKSINASI,” Unnes J. Math, vol. 6, no. 1, p. 26, 2017.

M. Torrea, J. L. Torrea, and D. Ortega, “A MODELING OF A DIPHTHERIA EPIDEMIC IN THE REFUGEES CAMPS”, bioRxiv, p. 208835, 2017. [Online]. Available: https://www.biorxiv.org/content/10.1101/208835v1, doi : https://doi.org/10.1101/208835

K. Sornbundit, W. Triampo, and C. Modchang, “MATHEMATICAL MODELING OF DIPHTHERIA TRANSMISSION IN THAILAND,” Comput. Biol. Med., vol. 87, pp. 162–168, 2017, doi: 10.1016/j.compbiomed.2017.05.031. doi : https://doi.org/10.1016/j.compbiomed.2017.05.031

F. Ilahi and A. Widiana, “THE EFFECTIVENESS OF VACCINE IN THE OUTBREAK OF DIPHTHERIA: MATHEMATICAL MODEL AND SIMULATION,” IOP Conf. Ser. Mater. Sci. Eng., vol. 434, no. 1, 2018, doi: 10.1088/1757-899X/434/1/012006. doi : https://doi.org/10.1088/1757-899X/434/1/012006

I. S. Fauzi et al., “ASSESSING THE IMPACT OF BOOSTER VACCINATION ON DIPHTHERIA TRANSMISSION: MATHEMATICAL MODELING AND RISK ZONE MAPPING,” Infect. Dis. Model., vol. 9, no. 1, pp. 245–262, 2024, doi: 10.1016/j.idm.2024.01.004. doi : https://doi.org/10.1016/j.idm.2024.01.004

S. O. Adewale, S. O. Ajao, I. A. Olopade, G. A. Adeniran, S. O. Ajao, and I. T. Mohammed, “MATHEMATICAL ANALYSIS OF QUARANTINE ON THE DYNAMICAL TRANSMISSION OF DIPHTHERIA DISEASE SOLVING RICCATI EQUATION USING ADOMIAN DECOMPOSITION METHOD VIEW PROJECT COVID PROJECT VIEW PROJECT IDAYAT TEMILADE MOHAMMED MATHEMATICAL ANALYSIS OF QUARANTINE ON THE,” Int. J. Sci. Eng. Investig., vol. 6, no. 64, pp. 8–17, 2017, [Online]. Available: www.IJSEI.com

N. Izzati, A. Andriani, and R. Robi’aqolbi, “OPTIMAL CONTROL OF DIPHTHERIA EPIDEMIC MODEL WITH PREVENTION AND TREATMENT,” J. Phys. Conf. Ser., vol. 1663, no. 1, 2020, doi: 10.1088/1742-6596/1663/1/012042.doi : https://doi.org/10.1088/1742-6596/1663/1/012042

N. Rahmi and M. I. Pratama, “MODEL ANALYSIS OF DIPHTHERIA DISEASE TRANSMISSION WITH VACCINATION, QUARANTINE, AND HAND-WASHING BEHAVIOR,” JTAM (Jurnal Teor. dan Apl. Mat., vol. 7, no. 2, p. 462, 2023, doi: 10.31764/jtam.v7i2.13466.doi : https://doi.org/10.31764/jtam.v7i2.13466

K. P. Putra and M. Rosha, “MODEL MATEMATIKA TIPE SIQR PENYEBARAN PENYAKIT DIFTERI DENGAN PENGARUH VAKSINASI,” J. Math. UNP, vol. 7, no. 4, pp. 84–93, 2022. doi : https://doi.org/10.24036/unpjomath.v7i4.13993

M. Sato, R. Ratianingsih, and Hajar, “MEMBANGUN MODEL MATEMATIKA PENYEBARAN PENYAKIT DIFTERI,” J. Ilm. Mat. dan Ter., vol. 18, no. 2, pp. 221–229, 2021.doi : https://doi.org/10.22487/2540766X.2021.v18.i2.15705

M. Sholeh, “STABILITY ANALYSIS OF THE SIQR MODEL OF DIPHTHERIA DISEASE SPREAD AND MIGRATION IMPACT,” BAREKENGJournal Math. Its Appl., vol. 19, no. 1, pp. 0173–0184, 2025. doi : https://doi.org/10.30598/barekengvol19iss1pp173-184

J. Xu and C. Xue, “UNIQUENESS AND EXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF FUNCTIONAL DIFFERENTIAL EQUATIONS,” AIMS Math., vol. 8, no. 1, pp. 676–690, 2023. doi : https://doi.org/10.3934/math.2023032

E. Sinaga, “PENERAPAN METODE LEAST SQUARES METHOD DALAM ESTIMASI PENJUALAN PRODUK ELEKTRONIK,” J. Comput. Informatics Res., vol. 2, no. 2, pp. 44–48, 2023. doi : https://doi.org/10.47065/comforch.v2i2.380

H. Motulsky and A. Christopoulos, Using Global Fitting to Test Treatment Effects. Oxford: Oxford University Press, 2023.

R. Prasetyo, “IMPLEMENTATION OF THE LEAST SQUARE METHOD FOR POPULATION GROWTH RATE PREDICTION IN AIR SUGIHAN DISTRICT,” J. Inf. Syst. Informatics, vol. 4, no. 2, pp. 282–299, 2022.doi : https://doi.org/10.51519/journalisi.v4i2.252

Central Statistics Agency of the Republic of Indonesia, “POPULATION BY AGE GROUP AND GENDER.” Accessed: Apr. 29, 2025. [Online]. Available: https://www.bps.go.id/id/statisticstable/3/WVc0MGEyMXBkVFUxY25KeE9HdDZkbTQzWkVkb1p6MDkjMw==/jumlah-penduduk-menurut-kelompok-umur-dan-jenis-kelamin--ribu-jiwa---2018.html?year=2023

K. W. D. Nugraha, Setiaji, F. Sibuea, and B. Hardhana, Profil Kesehatan Indonesia. Jakarta: Kementrian Kesehatan Republik Indonesia, 2021.

C. S. A. of the R. of Indonesia, “LIFE EXPECTANCY BY PROVINCE AND GENDER.” Accessed: Apr. 29, 2025. [Online]. Available: https://www.bps.go.id/id/statistics-table/2/NTAxIzI=/angka-harapan-hidup-ahh-menurut-provinsi-dan-jenis-kelamin.html

W. Y. Yang, W. Cao, J. Kim, K. W. Park, H. H. Park, J. Joung, Jong-Suk Ro, H. L. Lee, C. Hong, T. Im. Applied Numerical Methods Using MATLAB. John Wiley & Sons, 2020.doi : https://doi.org/10.1002/9781119626879

B. Wei, “STABILITY ANALYSIS OF EQUILIBRIUM POINT AND LIMIT CYCLE OF TWO-DIMENSIONAL NONLINEAR DYNAMICAL SYSTEMS—A TUTORIAL,” Appl. Sci., vol. 13, no. 2, 2023. doi : https://doi.org/10.3390/app13021136

R. Mahardika, Widowati, and Y. D. Sumanto, “ROUTH-HURWITZ CRITERION AND BIFURCATION METHOD FOR STABILITY ANALYSIS OF TUBERCULOSIS TRANSMISSION MODEL,” in Journal of Physics: Conference Series, Journal of Physics: Conference Series, 2019. doi : https://doi.org/10.1088/1742-6596/1217/1/012056

G. J. Olsder and J. W. van der Woude, MATHEMATICAL SYSTEMS THEORY, vol. 4. VSSD Netherlands, 2005.

R. Sundari and E. Apriliani, “CONSTRUCTION OF LYAPUNOV FUNCTIONS TO DETERMINE STABILITY,” J. Sains dan Seni ITS, vol. 6, no. 1, pp. 1–5, 2017. doi : https://doi.org/10.12962/j23373520.v6i1.22862

K. T. Alligood, T. D. Sauer, J. A. Yorke, and D. Chillingworth, “CHAOS: AN INTRODUCTION TO DYNAMICAL SYSTEMS,” SIAM Rev., vol. 40, no. 3, p. 732, 1998.

Published
2026-04-08
How to Cite
[1]
M. Tafrikan, F. Fatmawati, W. Windarto, and C. E. Madubueze, “A MATHEMATICAL MODEL OF DIPHTHERIA TRANSMISSION DYNAMICS WITH HETEROGENEOUS SUSCEPTIBILITY”, BAREKENG: J. Math. & App., vol. 20, no. 3, pp. 1967-1984, Apr. 2026.