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.

Identifying patterning behavior in a plant infestation of insect pests.

In this study, we developed a mechanistic model formulated as a system of reaction-diffusion equations (RDE) to explore the spatiotemporal dynamics of a theoretical pest with a tillering host plant in a controlled rectangular plant field. Local perturbation analysis, a recently developed method of analysis for wave propagation, was utilized to determine patterning regimes resulting from the local and global behaviors of the slow and fast diffusing components of the RDE system, respectively. Turing analysis was done to show that the RDE system does not exhibit Turing patterns. With bug mortality as the bifurcation parameter, regions with oscillations and stable coexistence of the pest and tillers were identified. Numerical simulations illustrate the patterning regimes in 1D and 2D settings. The oscillations suggest that recurrences in pest infestation is possible. Moreover, simulations showed that patterns produced in the model are strongly influenced by the pests' homogeneous dynamics inside the controlled environment.

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.

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.

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.

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.