A bimodal lognormal distribution model for the prediction of COVID-19 deaths