A MATHEMATICAL MODEL OF DIPHTHERIA TRANSMISSION DYNAMICS WITH HETEROGENEOUS SUSCEPTIBILITY
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
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.
Copyright (c) 2026 Mohamad Tafrikan, Fatmawati Fatmawati, Windarto Windarto, Chinwendu E. Madubueze

This work is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.
Authors who publish with this Journal agree to the following terms:
- Author retain copyright and grant the journal right of first publication with the work simultaneously licensed under a creative commons attribution license that allow others to share the work within an acknowledgement of the work’s authorship and initial publication of this journal.
- Authors are able to enter into separate, additional contractual arrangement for the non-exclusive distribution of the journal’s published version of the work (e.g. acknowledgement of its initial publication in this journal).
- Authors are permitted and encouraged to post their work online (e.g. in institutional repositories or on their websites) prior to and during the submission process, as it can lead to productive exchanges, as well as earlier and greater citation of published works.




1.gif)


