A multivariate statistical approach to predict COVID-19 count data with epidemiological interpretation and uncertainty quantification