Background/aim: The COVID-19 pandemic originated in Wuhan, China, in December 2019 and became one of the worst global health crises ever. While struggling with the unknown nature of this novel coronavirus, many researchers and groups attempted to project the progress of the pandemic using empirical or mechanistic models, each one having its drawbacks. The first confirmed cases were announced early in March, and since then, serious containment measures have taken place in Turkey. Materials and methods: Here, we present a different approach, a Bayesian negative binomial multilevel model with mixed effects, for the projection of the COVID-19 pandemic and we apply this model to the Turkish case. The model source code is available at https:// github.com/kansil/covid-19. We predicted the confirmed daily cases and cumulative numbers from June 6th to June 26th with 80%, 95%, and 99% prediction intervals (PI). Results: Our projections showed that if we continued to comply with the measures and no drastic changes were seen in diagnosis or management protocols, the epidemic curve would tend to decrease in this time interval. Also, the predictive validity analysis suggests that the proposed model projections should have a PI around 95% for the first 12 days of the projections. Conclusion: We expect that drastic changes in the course of COVID-19 in Turkey will cause the model to suffer in predictive validity, and this can be used to monitor the epidemic. We hope that the discussion on these projections and the limitations of the epidemiological forecasting will be beneficial to the medical community, and policy makers.Key words: COVID-19, pandemic, epidemiology, Bayesian regression, Turkey