We propose a fractional order model for human papillomavirus (HPV) dynamics, including the effects of vaccination and public health education on developing cervical cancer. First, we discuss the general structure of Caputo fractional derivatives and integrals. Next, we define the fractional HPV model using Caputo derivatives. The model equilibrium quantities, with their stability, are discussed based on the magnitude of the reproduction number. We compute and simulate numerical solutions of the presented fractional model using the Adams-Bashforth-Moulton scheme. Meanwhile, real data sourced from reports from the World Health Organization is used to establish the parameters and compute the basic reproduction number. We present figures of state variables for different fractional orders and the classical integer order. The impacts of vaccination and public health education are discussed through numerical simulations. From the results, we observe that an increase in both vaccination rates and public health education increases the quality of life, and thus, reduces disease burden and suffering in communities. The results also confirm that modeling HPV transmission dynamics using fractional derivatives includes history effects in the model, making the model further insightful and appropriate for studying HPV dynamics.
Citation: Simphiwe M. Simelane, Justin B. Munyakazi, Phumlani G. Dlamini, Oluwaseun F. Egbelowo. Projections of human papillomavirus vaccination and its impact on cervical cancer using the Caputo fractional derivative[J]. Mathematical Biosciences and Engineering, 2023, 20(7): 11605-11626. doi: 10.3934/mbe.2023515
[1] | Edoardo Beretta, Dimitri Breda . Discrete or distributed delay? Effects on stability of population growth. Mathematical Biosciences and Engineering, 2016, 13(1): 19-41. doi: 10.3934/mbe.2016.13.19 |
[2] | Cristeta U. Jamilla, Renier G. Mendoza, Victoria May P. Mendoza . Explicit solution of a Lotka-Sharpe-McKendrick system involving neutral delay differential equations using the r-Lambert W function. Mathematical Biosciences and Engineering, 2020, 17(5): 5686-5708. doi: 10.3934/mbe.2020306 |
[3] | Anthony Tongen, María Zubillaga, Jorge E. Rabinovich . A two-sex matrix population model to represent harem structure. Mathematical Biosciences and Engineering, 2016, 13(5): 1077-1092. doi: 10.3934/mbe.2016031 |
[4] | Azmy S. Ackleh, Keng Deng . Stability of a delay equation arising from a juvenile-adult model. Mathematical Biosciences and Engineering, 2010, 7(4): 729-737. doi: 10.3934/mbe.2010.7.729 |
[5] | J. M. Cushing, Simon Maccracken Stump . Darwinian dynamics of a juvenile-adult model. Mathematical Biosciences and Engineering, 2013, 10(4): 1017-1044. doi: 10.3934/mbe.2013.10.1017 |
[6] | John Cleveland . Basic stage structure measure valued evolutionary game model. Mathematical Biosciences and Engineering, 2015, 12(2): 291-310. doi: 10.3934/mbe.2015.12.291 |
[7] | Wei Feng, Michael T. Cowen, Xin Lu . Coexistence and asymptotic stability in stage-structured predator-prey models. Mathematical Biosciences and Engineering, 2014, 11(4): 823-839. doi: 10.3934/mbe.2014.11.823 |
[8] | Santanu Bhattacharya, Nandadulal Bairagi . Dynamic optimization of fishing tax and tourism fees for sustainable bioeconomic resource management. Mathematical Biosciences and Engineering, 2025, 22(7): 1751-1789. doi: 10.3934/mbe.2025064 |
[9] | Lorenzo Mari, Marino Gatto, Renato Casagrandi . Local resource competition and the skewness of the sex ratio: a demographic model. Mathematical Biosciences and Engineering, 2008, 5(4): 813-830. doi: 10.3934/mbe.2008.5.813 |
[10] | Hua Guo, Yuan Tian, Kaibiao Sun, Xinyu Song . Dynamic analysis of two fishery capture models with a variable search rate and fuzzy biological parameters. Mathematical Biosciences and Engineering, 2023, 20(12): 21049-21074. doi: 10.3934/mbe.2023931 |
We propose a fractional order model for human papillomavirus (HPV) dynamics, including the effects of vaccination and public health education on developing cervical cancer. First, we discuss the general structure of Caputo fractional derivatives and integrals. Next, we define the fractional HPV model using Caputo derivatives. The model equilibrium quantities, with their stability, are discussed based on the magnitude of the reproduction number. We compute and simulate numerical solutions of the presented fractional model using the Adams-Bashforth-Moulton scheme. Meanwhile, real data sourced from reports from the World Health Organization is used to establish the parameters and compute the basic reproduction number. We present figures of state variables for different fractional orders and the classical integer order. The impacts of vaccination and public health education are discussed through numerical simulations. From the results, we observe that an increase in both vaccination rates and public health education increases the quality of life, and thus, reduces disease burden and suffering in communities. The results also confirm that modeling HPV transmission dynamics using fractional derivatives includes history effects in the model, making the model further insightful and appropriate for studying HPV dynamics.
[1] | CDC fact sheet. Centers for Disease Control and Prevention, Genital HPV infection. Available from: http://www.cdc.gov/std/HPV/STDFact-HPV.htm. |
[2] |
S. L. Lee, A. M. Tameru, A mathematical model of Human Papillomavirus (HPV) in the United States and its impact on Cervical Cancer, J. Cancer, 3 (2012), 262–268. https://doi.org/10.7150/jca.4161 doi: 10.7150/jca.4161
![]() |
[3] | National Cancer Institute, HPV and Cancer. Available from: https://www.cancer.gov/about-cancer/causes-prevention/risk/infectious-agents/hpv-and-cancer. |
[4] | World Health Organization, Cervical Cancer Fact Sheet, Access date: 02 June 2022. Available from: https://www.who.int/news-room/fact-sheets/detail/cervical-cancer. |
[5] | K. P. Braaten, M. R. Laufer, Human Papillomavirus (HPV), HPV-related disease, and the HPV vaccine, Rev. Obstet. Gynecol., 1 (2008), 2–10. |
[6] | Centers for Disease Control and Prevention, Human Papillomavirus (HPV) Vaccination: What Everyone Should Know, Accesse date: 02 June 2022. https://www.cdc.gov/vaccines/vpd/hpv/public/index.html. |
[7] |
U. K. Nwajeri, A. B. Panle, A. Omame, M. C. Obi, C. P. Onyenegecha, On the fractional order model for HPV and Syphilis using non–singular kernel, Results Phys., 37 (2022), 105463. https://doi.org/10.1016/j.rinp.2022.105463 doi: 10.1016/j.rinp.2022.105463
![]() |
[8] | I. Podlubny, Fractional Differential Equations. An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of their Solution and Some of their Applications, Academic Press, San Diego, CA, 1999. |
[9] |
D. Baleanu, M. Hassan Abadi, A. Jajarmi, K. Zarghami Vahid, J. J. Nieto, A new comparative study on the general fractional model of COVID-19 with isolation and quarantine effects, Alexandria Eng. J., 61 (2022), 4779–4791. https://doi.org/10.1016/j.aej.2021.10.030 doi: 10.1016/j.aej.2021.10.030
![]() |
[10] |
P. A. Naik, J. Zu, K. M. Owolabi, Global dynamics of a fractional order model for the transmission of HIV epidemic with optimal control, Chaos, Solitons Fractals, 138 (2020), 109826. https://doi.org/10.1016/j.chaos.2020.109826 doi: 10.1016/j.chaos.2020.109826
![]() |
[11] |
A. Atangana, E. Alabaraoye, Solving a system of fractional partial differential equations arising in the model of HIV infection of CD4+ cells and attractor one-dimensional Keller-Segel equations, Adv. Differ. Equations, 94 (2013). https://doi.org/10.1186/1687-1847-2013-94 doi: 10.1186/1687-1847-2013-94
![]() |
[12] |
R. Scherer, S. L. Kalla, Y. Tang, J. Huang, The Grunwald-Letnikov method for fractional differential equations, Comput. Math. Appl., 62 (2011), 902–917. https://doi.org/10.1016/j.camwa.2011.03.054 doi: 10.1016/j.camwa.2011.03.054
![]() |
[13] |
U. K. Nwajeri, A. Omame, C. P. Onyenegecha, Analysis of a fractional order model for HPV and CT co-infection, Results Phys., 28 (2021), 104643. https://doi.org/10.1016/j.rinp.2021.104643 doi: 10.1016/j.rinp.2021.104643
![]() |
[14] | K. M. Owolabi, A. Atangana, Numerical Methods for Fractional Differentiation, Springer Singapore, 2019. https://doi.org/10.1007/978-981-15-0098-5 |
[15] |
U. K. Nwajeri, A. B. Panle, A. Omame, M. C. Obi, C. P. Onyenegecha, On the fractional order model for HPV and Syphilis using non–singular kernel, Results Phys., 37 (2022), 105463. https://doi.org/10.1016/j.rinp.2022.105463 doi: 10.1016/j.rinp.2022.105463
![]() |
[16] |
K. Diethelm, N. J. Ford, Analysis of fractional differential equations, J. Math. Anal. Appl., 265 (2002), 229–248. https://doi.org/10.1006/jmaa.2000.7194 doi: 10.1006/jmaa.2000.7194
![]() |
[17] |
S. M. Simelane, P. G. Dlamini, A fractional order differential equation model for Hepatitis B virus with saturated incidence, Results Phys., 24 (2021), 104114. https://doi.org/10.1016/j.rinp.2021.104114 doi: 10.1016/j.rinp.2021.104114
![]() |
[18] |
Z. M. Odibat, N. T. Shawagfeh, Generalized Taylor's formula, Appl. Math. Comput., 186 (2007), 286–293. https://doi.org/10.1016/j.amc.2006.07.102 doi: 10.1016/j.amc.2006.07.102
![]() |
[19] |
N. Ozalp, E. Demirci, A fractional order SEIR model with vertical transmission, Math. Comput. Model., 54 (2011), 1–6. https://doi.org/10.1016/j.mcm.2010.12.051 doi: 10.1016/j.mcm.2010.12.051
![]() |
[20] |
W. Lin, Global existence theory and chaos control of fractional differential equations, J. Math. Anal. Appl., 332 (2007), 709–726. https://doi.org/10.1016/j.jmaa.2006.10.040 doi: 10.1016/j.jmaa.2006.10.040
![]() |
[21] |
O. Diekmann, J. A. P. Heesterbeek, J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations, J. Math. Biol., 28 (1990), 365–382. https://doi.org/10.1007/BF00178324 doi: 10.1007/BF00178324
![]() |
[22] | J. Jia, J. Xiao, Stability analysis of a disease resistance seirs model with nonlinear incidence rate, Adv. Differ. Equations, 75 (2018). https://doi.org/10.1186/s13662-018-1494-1 |
[23] |
P. van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
![]() |
[24] | J. L. S. Allen, Introduction to Mathematical Biology, Pearson Education, Inc, Upper Saddle River, 2007. |
[25] | MATLAB. 9.7.0.1190202 (R2019b). Natick, Massachusetts: The MathWorks Inc, 2018. |
[26] | Wolfram Research, Inc, Mathematica, Version 9.0, Champaign, IL, 2012. |
[27] | Indian Council of Medical Research, Consensus Document for Management of Cancer Cervix, 2016. |
[28] | J. Ferlay, M. Ervik, F. Lam, M. Colombet, L. Mery, M. Pineros, et al., Global Cancer Observatory: Cancer Today. Lyon, France: International Agency for Research on Cancer, 2020. Available from: https://gco.iarc.fr/today. |
[29] |
E. Amponsah-Dacosta, N. Blose, V. V. Nkwinika, V. Chepkurui, Human papillomavirus vaccination in South Africa: Programmatic challenges and opportunities for integration with other adolescent health services, Front. Public Health, 10 (2022). https://doi.org/10.3389/fpubh.2022.799984 doi: 10.3389/fpubh.2022.799984
![]() |
[30] | United Nations, Department of Economic and Social Affairs, Population Division, World Population Prospects 2019, Online Edition. Rev. 1, Access date: 02 January 2021, Available from: https://population.un.org/wpp/Download/Standard/Population/. |
[31] | World Health Organization, Global strategy to accelerate the elimination of cervical cancer as a public health problem, Access date: 05 March 2021, Available from: https://www.who.int/publications/i/item/9789240014107. |
[32] |
N. Gupta, A. S. Chauhan, S. Prinja, A. K. Pandey, Impact of COVID-19 on outcomes for patients with cervical cancer in India, JCO Global Oncol., 7 (2021), 716–725. https://doi.org/10.1200/GO.20.00654 doi: 10.1200/GO.20.00654
![]() |
[33] | L. Bruni, G. Albero, B. Serrano, M. Mena, J. J. Collado, D. Gómez, et al., ICO/IARC Information Centre on HPV and Cancer (HPV Information Centre), Human Papillomavirus and Related Diseases in India, Summary Report, Access date: 18 June 2022. |
1. | Danielle A. Martin, Jane Hook, Sunny Wonsun Kim, Linda Larkey, Rebecca E. Lee, Mind–body therapy for cardiometabolic risk in U.S. middle-aged Black adults: a scoping review, 2025, 13, 2296-2565, 10.3389/fpubh.2025.1480369 |