Poly(lactic-co-glycolic acid) is one of the most used polymers for drug delivery systems (DDSs). It shows excellent biocompatibility, biodegradability, and allows spatio-temporal control of the release of a drug by altering its chemistry. In spite of this, few formulations have reached the market. To characterize and optimize the drug release process, mathematical models offer a good alternative as they allow interpreting and predicting experimental findings, saving time and money. However, there is no general model that describes all types of drug release of polymeric DDSs. This study aims to perform a statistical comparison of several mathematical models commonly used in order to find which of them best describes the drug release profile from PLGA particles synthesized by nanoprecipitation method. For this purpose, 40 datasets extracted from scientific articles published since 2016 were collected. Each set was fitted by the models: order zero to fifth order polynomials, Korsmeyer-Peppas, Weibull and Hyperbolic Tangent Function. Some data sets had few observations that do not allow to apply statistic test, thus bootstrap resampling technique was performed. Statistic evidence showed that Hyperbolic Tangent Function model is the one that best fit most of the data.