Ecologists have long sought a way to explain how the remarkable biodiversity observed in nature is maintained. On the one hand, simple models of interacting competitors cannot produce the stable persistence of very large ecological communities(1-5). On the other hand, neutral models(6-9), in which species do not interact and diversity is maintained by immigration and speciation, yield unrealistically small fluctuations in population abundance(10), and a strong positive correlation between a species' abundance and its age(11), contrary to empirical evidence. Models allowing for the robust persistence of large communities of interacting competitors are lacking. Here we show that very diverse communities could persist thanks to the stabilizing role of higher-order interactions(12,13), in which the presence of a species influences the interaction between other species. Although higher-order interactions have been studied for decades(14-16), their role in shaping ecological communities is still unclear(5). The inclusion of higher-order interactions in competitive network models stabilizes dynamics, making species coexistence robust to the perturbation of both population abundance and parameter values. We show that higher-order interactions have strong effects in models of closed ecological communities, as well as of open communities in which new species are constantly introduced. In our framework, higher-order interactions are completely defined by pairwise interactions, facilitating empirical parameterization and validation of our models.