The functional response describes resource-dependent feeding rates of consumers. While models for feeding on a single resource species are well studied and supported by a large body of empirical research, consumers feeding on multiple resource species are ubiquitous in nature. However, feeding experiments designed for parameterizing multi-species functional responses (MSFR) are extremely rare, mainly due to logistical challenges and the non-trivial nature of their statistical analysis.Here, we describe how these models can be fitted to empirical data in a Bayesian framework. Specifically, we address the problem of prey depletion during experiments, which can be accounted for through dynamical modeling. In a comprehensive simulation study, we test the effects of experimental design, sample size and noise level on the identifiability of four distinct MSFR models. Additionally, we demonstrate the method’s versatility by applying it to a list of empirical datasets.We identify experimental designs for feeding trials that produce the most accurate parameter estimates in two- and three-prey scenarios. Although noise can introduce systematic bias in parameter estimates, model selection performs surprisingly well, almost always identifying the correct model even for small datasets.Our flexible method allows the simultaneous analysis of feeding experiments from both single- and multi-prey scenarios, either with or without prey depletion. This will help to improve our understanding of prey selectivity, prey switching and hence food web stability and biodiversity. Our approach equips researchers with the right statistical tools to improve the mechanistic understanding of feeding interactions in complex ecosystems.