The role of isospin in quasielastic electron scattering and charge-changing neutrino reactions is investigated in the relativistic impulse approximation. We analyze proton and neutron scaling functions making use of various theoretical descriptions for the final-state interactions, focusing on the effects introduced by the presence of strong scalar and vector terms in the relativistic mean field approach. An explanation for the differences observed in the scaling functions evaluated from (e, e ) and (ν, μ) reactions is provided by invoking the differences in isoscalar and isovector contributions. © 2007 Elsevier B.V. All rights reserved. PACS: 25.30.Pt; 25.30.Fj; 24.10.Jv Extensive analyses of quasielastic (QE) inclusive electron scattering data performed in recent years [1][2][3][4] have clearly demonstrated the quality of the behavior known as scaling. These analyses are based on the so-called superscaling function, f (ψ ), obtained by dividing the cross section by an appropriate function which contains the single nucleon physics, and plotting the result against the scaling variable ψ (q, ω) (see, e.g., [5]). One then studies the dependences upon the momentum transfer q and the specific nucleus chosen. From (e, e ) world data one concludes that scaling of the first kind (no dependence on q) is reasonably respected at energies below the QE peak, whereas scaling of the second kind (no dependence on the nuclear species) is fulfilled very well in the same region. The simultaneous occurrence of both kinds of scaling is called superscaling. At energies above the QE peak, breaking of both * Corresponding author.E-mail address: jac@us.es (J.A. Caballero).kinds of scaling is observed, residing mostly in the transverse channel, and likely due to effects beyond the impulse approximation.Experimental data lead to a scaling function with a characteristic asymmetric shape, having a tail that extends to high values of the transferred energy ω (positive values of the scaling variable ψ (q, ω)). The asymmetric shape of the scaling function is largely absent in non-relativistic (NR) models based on a mean field approach. In contrast, the study presented in [6,7] has shown that the asymmetry can in fact be obtained within the relativistic impulse approximation (RIA), given that a description of final-state interactions (FSI) using strong relativistic mean field (RMF) potentials is assumed. Recently, we have shown [8] that an asymmetrical scaling function can be also obtained within the framework of a semi-relativistic (SR) model, based on improved NR expansions of the on-shell electromagnetic current, provided that the FSI are described by the Dirac equation-based (DEB) potential [9] derived from the RMF. Note that, in the SR model, the nonlocalities arising from the 0370-2693/$ -see front matter