The H mode of confinement in Tokamaks is characterized by a thin region of high gradients, located at the edge of the plasma and called the Edge Transport Barrier. Even if various theoretical models have been proposed for the interpretation of the edge physics, the main empirical scaling laws of the plasma confinement time are expressed in terms of global plasma parameters and they do not discriminate between the edge and core regions. Moreover all the scaling laws are assumed to be power law monomials. In the present paper, a new methodology is proposed to investigate the validity of both assumptions. The approach is based on Symbolic Regression via Genetic Programming and allows first the extraction of the most statistically reliable models from the available experimental data in the ITPA database. Non linear fitting is then applied to the mathematical expressions found by Symbolic regression. The obtained scaling laws are compared with the traditional scalings in power law form.