A model-based strategy for the COVID-19 vaccine roll-out in the Philippines.

COVID-19 has affected millions of people worldwide, causing illness and death, and disrupting daily life while imposing a significant social and economic burden. Vaccination is an important control measure that significantly reduces mortality if properly and efficiently distributed. In this work, an age-structured model of COVID-19 transmission, incorporating an unreported infectious compartment, is developed. Three age groups are considered: young (0-19 years), adult (20-64 years), and elderly (65+ years). The transmission rate and reporting rate are determined for each group by utilizing the number of COVID-19 cases in the National Capital Region in the Philippines. Optimal control theory is employed to identify the best vaccine allocation to different age groups. Further, three different vaccination periods are considered to reflect phases of vaccination priority groups: the first, second, and third account for the inoculation of the elderly, adult and elderly, and all three age groups, respectively. This study could guide in making informed decisions in mitigating a population-structured disease transmission under limited resources.

Risk estimation of lifted mask mandates and emerging variants using mathematical model.

Background: More than half of the population in Korea had a prior COVID-19 infection. In 2022, most nonpharmaceutical interventions, except mask-wearing indoors, had been lifted. And in 2023, the indoor mask mandates were eased. Methods: We developed an age-structured compartmental model that distinguishes vaccination history, prior infection, and medical staff from the rest of the population. Contact patterns among hosts were separated based on age and location. We simulated scenarios with the lifting of the mask mandate all at once or sequentially according to the locations. Furthermore, we investigated the impact of a new variant assuming that it has higher transmissibility and risk of breakthrough infection. Results: We found that the peak size of administered severe patients may not exceed 1100 when the mask mandate is lifted everywhere, and 800 if the mask mandate only remains in the hospital. If the mask mandate is lifted in a sequence (except hospital), then the peak size of administered severe patients may not exceed 650. Moreover, if the new variant has both higher transmissibility and immune reduction, the effective reproductive number of the new variant is approximately 3 times higher than that of the current variant, and additional interventions may be needed to keep the administered severe patients from exceeding 2,000, which is the critical level we set. Conclusion: Our findings showed that the lifting of the mask mandate, except in hospitals, would be more manageable if implemented sequentially. Considering a new variant, we found that depending on the population immunity and transmissibility of the variant, wearing masks and other interventions may be necessary for controlling the disease.

Modelling the effects of social distancing, antiviral therapy, and booster shots on mitigating Omicron spread.

As the COVID-19 situation changes because of emerging variants and updated vaccines, an elaborate mathematical model is essential in crafting proactive and effective control strategies. We propose a COVID-19 mathematical model considering variants, booster shots, waning, and antiviral drugs. We quantify the effects of social distancing in the Republic of Korea by estimating the reduction in transmission induced by government policies from February 26, 2021 to February 3, 2022. Simulations show that the next epidemic peak can be estimated by investigating the effects of waning immunity. This research emphasizes that booster vaccination should be administered right before the next epidemic wave, which follows the increasing waned population. Policymakers are recommended to monitor the waning population immunity using mathematical models or other predictive methods. Moreover, our simulations considering a new variant's transmissibility, severity, and vaccine evasion suggest intervention measures that can reduce the severity of COVID-19.

Estimation of monkeypox spread in a nonendemic country considering contact tracing and self-reporting: A stochastic modeling study.

In May 2022, monkeypox started to spread in nonendemic countries. To investigate contact tracing and self-reporting of the primary case in the local community, a stochastic model is developed. An algorithm based on Gillespie's stochastic chemical kinetics is used to quantify the number of infections, contacts, and duration from the arrival of the primary case to the detection of the index case (or until there are no more local infections). Different scenarios were set considering the delay in contact tracing and behavior of infectors. We found that the self-reporting behavior of a primary case is the most significant factor affecting outbreak size and duration. Scenarios with a self-reporting primary case have an 86% reduction in infections (average: 5-7, in a population of 10 000) and contacts (average: 27-72) compared with scenarios with a non-self-reporting primary case (average number of infections and contacts: 27-72 and 197-537, respectively). Doubling the number of close contacts per day is less impactful compared with the self-reporting behavior of the primary case as it could only increase the number of infections by 45%. Our study emphasizes the importance of the prompt detection of the primary case.

Optimal selection of COVID-19 vaccination sites in the Philippines at the municipal level.

In this work, we present an approach to determine the optimal location of coronavirus disease 2019 (COVID-19) vaccination sites at the municipal level. We assume that each municipality is subdivided into smaller administrative units, which we refer to as barangays. The proposed method solves a minimization problem arising from a facility location problem, which is formulated based on the proximity of the vaccination sites to the barangays, the number of COVID-19 cases, and the population densities of the barangays. These objectives are formulated as a single optimization problem. As an alternative decision support tool, we develop a bi-objective optimization problem that considers distance and population coverage. Lastly, we propose a dynamic optimization approach that recalculates the optimal vaccination sites to account for the changes in the population of the barangays that have completed their vaccination program. A numerical scheme that solves the optimization problems is presented and the detailed description of the algorithms, which are coded in Python and MATLAB, are uploaded to a public repository. As an illustration, we apply our method to determine the optimal location of vaccination sites in San Juan, a municipality in the province of Batangas, in the Philippines. We hope that this study may guide the local government units in coming up with strategic and accessible plans for vaccine administration.

Multi-Faceted Analysis of COVID-19 Epidemic in Korea Considering Omicron Variant: Mathematical Modeling-Based Study.

BACKGROUND: The most recent variant of concern, omicron (B.1.1.529), has caused numerous cases worldwide including the Republic of Korea due to its fast transmission and reduced vaccine effectiveness. METHODS: A mathematical model considering age-structure, vaccine, antiviral drugs, and influx of the omicron variant was developed. We estimated transmission rates among age groups using maximum likelihood estimation for the age-structured model. The impact of non-pharmaceutical interventions (NPIs; in community and border), quantified by a parameter µ in the force of infection, and vaccination were examined through a multi-faceted analysis. A theory-based endemic equilibrium study was performed to find the manageable number of cases according to omicron- and healthcare-related factors. RESULTS: By fitting the model to the available data, the estimated values of µ ranged from 0.31 to 0.73, representing the intensity of NPIs such as social distancing level. If µ < 0.55 and 300,000 booster shots were administered daily from February 3, 2022, the number of severe cases was forecasted to exceed the severe bed capacity. Moreover, the number of daily cases is reduced as the timing of screening measures is delayed. If screening measure was intensified as early as November 24, 2021 and the number of overseas entrant cases was contained to 1 case per 10 days, simulations showed that the daily incidence by February 3, 2022 could have been reduced by 87%. Furthermore, we found that the incidence number in mid-December 2021 exceeded the theory-driven manageable number of daily cases. CONCLUSION: NPIs, vaccination, and antiviral drugs influence the spread of omicron and number of severe cases in the Republic of Korea. Intensive and early screening measures during the emergence of a new variant is key in controlling the epidemic size. Using the endemic equilibrium of the model, a formula for the manageable daily cases depending on the severity rate and average length of hospital stay was derived so that the number of severe cases does not surpass the severe bed capacity.

Optimal placement of tsunami sensors with depth constraint.

Tsunamis are destructive natural disasters that can cause severe damage to property and the loss of many lives. To mitigate the damage and casualties, tsunami warning systems are implemented in coastal areas, especially in locations with high seismic activity. This study presents a method to identify the placement of near-shore detection sensors by minimizing the tsunami detection time, obtained by solving the two-dimensional shallow water equations (SWE). Several benchmark tests were done to establish the robustness of the SWE model, which is solved using a staggered finite volume method. The optimization problem is solved using particle swarm optimization (PSO). The proposed method is applied to different test problems. As an application, the method is used to find the optimal location of a detection sensor using data from the 2018 Palu tsunami. Our findings show that detection time can be significantly reduced through the strategic placement of tsunami sensors.

A Baybayin word recognition system.

Baybayin is a pre-Hispanic Philippine writing system used in Luzon island. With the effort in reintroducing the script, in 2018, the Committee on Basic Education and Culture of the Philippine Congress approved House Bill 1022 or the "National Writing System Act," which declares the Baybayin script as the Philippines' national writing system. Since then, Baybayin OCR has become a field of research interest. Numerous works have proposed different techniques in recognizing Baybayin scripts. However, all those studies anchored on the classification and recognition at the character level. In this work, we propose an algorithm that provides the Latin transliteration of a Baybayin word in an image. The proposed system relies on a Baybayin character classifier generated using the Support Vector Machine (SVM). The method involves isolation of each Baybayin character, then classifying each character according to its equivalent syllable in Latin script, and finally concatenate each result to form the transliterated word. The system was tested using a novel dataset of Baybayin word images and achieved a competitive 97.9% recognition accuracy. Based on our review of the literature, this is the first work that recognizes Baybayin scripts at the word level. The proposed system can be used in automated transliterations of Baybayin texts transcribed in old books, tattoos, signage, graphic designs, and documents, among others.

A mathematical model of the dynamics of lymphatic filariasis in Caraga Region, the Philippines.

Despite being one of the first countries to implement mass drug administration (MDA) for elimination of lymphatic filariasis (LF) in 2001 after a pilot study in 2000, the Philippines is yet to eliminate the disease as a public health problem with 6 out of the 46 endemic provinces still implementing MDA for LF as of 2018. In this work, we propose a mathematical model of the transmission dynamics of LF in the Philippines and a control strategy for its elimination using MDA. Sensitivity analysis using the Latin hypercube sampling and partial rank correlation coefficient methods suggests that the infected human population is most sensitive to the treatment parameters. Using the available LF data in Caraga Region from the Philippine Department of Health, we estimate the treatment rates r 1 and r 2 using the least-squares parameter estimation technique. Parameter bootstrapping showed small variability in the parameter estimates. Finally, we apply optimal control theory with the objective of minimizing the infected human population and the corresponding implementation cost of MDA, using the treatment coverage γ as the control parameter. Simulation results highlight the importance of maintaining a high MDA coverage per year to effectively minimize the infected population by the year 2030.

Optical character recognition system for Baybayin scripts using support vector machine.

In 2018, the Philippine Congress signed House Bill 1022 declaring the Baybayin script as the Philippines' national writing system. In this regard, it is highly probable that the Baybayin and Latin scripts would appear in a single document. In this work, we propose a system that discriminates the characters of both scripts. The proposed system considers the normalization of an individual character to identify if it belongs to Baybayin or Latin script and further classify them as to what unit they represent. This gives us four classification problems, namely: (1) Baybayin and Latin script recognition, (2) Baybayin character classification, (3) Latin character classification, and (4) Baybayin diacritical marks classification. To the best of our knowledge, this is the first study that makes use of Support Vector Machine (SVM) for Baybayin script recognition. This work also provides a new dataset for Baybayin, its diacritics, and Latin characters. Classification problems (1) and (4) use binary SVM while (2) and (3) apply the multiclass SVM classification. On average, our numerical experiments yield satisfactory results: (1) has 98.5% accuracy, 98.5% precision, 98.49% recall, and 98.5% F1 Score; (2) has 96.51% accuracy, 95.62% precision, 95.61% recall, and 95.62% F1 Score; (3) has 95.8% accuracy, 95.85% precision, 95.8% recall, and 95.83% F1 Score; and (4) has 100% accuracy, 100% precision, 100% recall, and 100% F1 Score.

Explicit solution of a Lotka-Sharpe-McKendrick system involving neutral delay differential equations using the r-Lambert W function.

Structured population models, which account for the state of individuals given features such as age, gender, and size, are widely used in the fields of ecology and biology. In this paper, we consider an age-structured population model describing the population of adults and juveniles. The model consists of a system of ordinary and neutral delay differential equations. We present an explicit solution to the model using a generalization of the Lambert W function called the r-Lambert W function. Numerical simulations with varying parameters and initial conditions are done to illustrate the obtained solution. The proposed method is also applied to an insect population model with long larval and short adult phases.

Application of particle swarm optimization in optimal placement of tsunami sensors.

Rapid detection and early warning systems demonstrate crucial significance in tsunami risk reduction measures. So far, several tsunami observation networks have been deployed in tsunamigenic regions to issue effective local response. However, guidance on where to station these sensors are limited. In this article, we address the problem of determining the placement of tsunami sensors with the least possible tsunami detection time. We use the solutions of the 2D nonlinear shallow water equations to compute the wave travel time. The optimization problem is solved by implementing the particle swarm optimization algorithm. We apply our model to a simple test problem with varying depths. We also use our proposed method to determine the placement of sensors for early tsunami detection in Cotabato Trench, Philippines.