  We consider the recent data sets of quasi-periodic oscillations from eight
different low mass X-ray binaries. We here interpret their physical features in
the context of given regular black hole solutions and verify their
applicability to neutron star configurations. We evaluate the numerical
constraints over the free parameters of Bardeen, Hayward and Dymnikova regular
solutions by performing a set of Markov chain Monte Carlo analyses, based on
the Metropolis algorithm. For each source, we evaluate the best-fit parameters,
among which mass and magnetic charge, and compare and contrast them with the
current literature. We also infer the corresponding innermost stable circular
orbit radii and the radial extents of the accretion disks. Focusing on how to
identify discrepancies among theoretical models and observations, our results
show that, in most of the cases, regular black holes, in particular the Bardeen
and Hayward spacetimes are slightly more suitable to describe neutron stars
than Schwarzschild geometry, whereas the Dymnikova metric is ruled out.
