This is an existing R function I modified for my research. It allows for estimating the distribution of data, bootstrapping that data, then displaying the cluster center of those bootstrapped values for the most accurate estimation of the underlying distribution.
Original R function I modified authored by: Marie-Laure Delignette-Muller [aut], Christophe Dutang [aut], Regis Pouillot [ctb], Jean-Baptiste Denis [ctb], Aurelie Siberchicot [aut, cre], https://cran.r-project.org/web/packages/fitdistrplus/index.html.