  We provide a probabilistic SIRD model for the COVID-19 pandemic in Italy,
where we allow the infection, recovery and death rates to be random. In
particular, the underlying random factor is driven by a fractional Brownian
motion. Our model is simple and needs only some few parameters to be
calibrated.
