Dynamic behavior of a stochastic epidemic model for skin sores: theoretical and computational perspectives

Why tiny skin infections matter

Small skin sores might seem like a minor nuisance, but in many remote communities they are a major and recurring health problem, especially for children. These infections can spread quickly, cause serious complications, and are difficult to monitor in places with limited medical resources. This study uses mathematics and computer simulations to better understand how skin sores spread through a community, how chance events can help or hinder an outbreak, and how to simulate these patterns reliably on a computer so that future control strategies can be planned more wisely.

Turning skin sores into a simple story

The authors start by describing the population with three groups: people who are at risk of infection, people currently infected with skin sores, and people who have recovered. They write equations that describe how people move between these groups as they are infected, recover, or die from natural causes. By studying these equations, they show that the model behaves sensibly: the numbers of people never become negative or explode to impossible values, and under some conditions the infection dies out while under others it can settle into a steady presence in the community. A key quantity in this analysis is a threshold number that compares how quickly the disease spreads with how quickly people recover or leave the population; values below this threshold mean skin sores fade away, while values above it mean the disease can persist.

Adding the role of chance

Real life is messy: contact patterns change day to day, environments vary, and outbreaks do not follow perfectly smooth curves. To reflect this, the authors extend their original model by adding random fluctuations to the infection and recovery processes. These random influences are represented by a mathematical object called Brownian motion, which captures unpredictable ups and downs over time. In this “stochastic” version of the model, the same disease can either take hold or fizzle out depending on how chance events line up. The researchers prove that, even with randomness, the model still keeps population sizes realistic and positive, and they identify a new threshold that combines the usual spread rate with the strength of random effects. If this adjusted threshold is below one, the infection almost surely disappears; if it is above one, it tends to persist.

How randomness can actually help

One striking conclusion is that randomness can sometimes help eliminate the disease. When the infection is already near the tipping point where it might or might not persist, the extra “noise” from changing contacts or environments can push the system toward extinction of the infection. In effect, random dips in transmission can outweigh bursts of spread. The analysis shows that this happens because the random swings reduce the average growth rate of the infected group. For public health, this means that improving hygiene, reducing crowding, or treating cases early—measures that make infections more fragile—can combine with natural variability to tip the balance toward disease elimination, even in settings where traditional models would predict ongoing transmission.

Putting the model to the test on a computer

To explore these ideas in practice, the authors compare several ways of simulating the model on a computer. They test standard approaches, like the basic Euler scheme and the more accurate Runge–Kutta method, alongside a more recent technique called a nonstandard finite difference (NSFD) method adapted to random effects. By running the same model with different time steps, they show that the NSFD method is especially good at keeping the simulated populations positive and stable, even when the time steps are very large. The Runge–Kutta method performs extremely well with small time steps but becomes less convenient when simulations must run over long periods. The simpler Euler method works only when the time steps are modest; with coarse steps it can distort the dynamics.

What this means for fighting skin sores

Overall, the study shows that a carefully built mathematical model can capture both the average spread of skin sores and the unpredictable twists that shape real outbreaks. The work demonstrates that random variation is not just noise to be ignored, but a crucial factor that can shift the long-term fate of the disease. It also highlights that the choice of computer method matters: the nonstandard finite difference approach gives reliable, stable simulations across a wide range of settings. For decision makers and researchers, this combination of theory and computation offers a more trustworthy way to explore how skin sores might spread, when they are likely to die out, and which control measures are most likely to push communities toward a future with fewer infections.

Citation: Raza, A., Lampart, M., Rocha, E.M. et al. Dynamic behavior of a stochastic epidemic model for skin sores: theoretical and computational perspectives. Sci Rep 16, 14091 (2026). https://doi.org/10.1038/s41598-026-43704-4

Keywords: skin sores, stochastic epidemic model, numerical simulation, disease transmission, remote communities