  We present a global analysis of the inclusive quasielastic electron
scattering data with a superscaling approach with relativistic effective mass.
The SuSAM* model exploits the approximation of factorization of the scaling
function $f^*(\psi^*)$ out of the cross section under quasifree conditions. Our
approach is based on the relativistic mean field theory of nuclear matter where
a relativistic effective mass for the nucleon encodes the dynamics of nucleons
moving in presence of scalar and vector potentials. Both the scaling variable
$\psi^*$ and the single nucleon cross sections include the effective mass as a
parameter to be fitted to the data alongside the Fermi momentum $k_F$. Several
methods to extract the scaling function and its uncertainty from the data are
proposed and compared. The model predictions for the quasielastic cross section
and the theoretical error bands are presented and discussed for nuclei along
the periodic table from $A=2$ to $A=238$: $^2$H, $^3$H, $^3$He, $^4$He,
$^{12}$C, $^{6}$Li, $^{9}$Be, $^{24}$Mg, $^{59}$Ni,
  $^{89}$Y, $^{119}$Sn, $^{181}$Ta, $^{186}$W, $^{197}$Au, $^{16}$O, $^{27}$Al,
$^{40}$Ca, $^{48}$Ca, $^{56}$Fe, $^{208}$Pb, and $^{238}$U.
  We find that more than 9000 of the total $\sim 20000$ data fall within the
quasielastic theoretical bands. Predictions for $^{48}$Ti and $^{40}$Ar are
also provided for the kinematics of interest to neutrino experiments.
