  The deviations from full or pure aging behavior, i.e. perfect $t/t_w$ scaling
of the correlation and response functions of aging glassy systems, are not well
understood theoretically. Recent experiments of Rodriguez et al. (Phys. Rev.
Lett. 91, 037203 (2003)) have shown that full aging applies to the
thermoremanent magnetization in the limit of infinite cooling rate during the
initial quench. In numerical models, instantaneous thermal quenches can be--and
usually are--applied. This has motivated the present numerical investigations
of the aging behavior of Edwards-Anderson Ising spin glass models with both
short- and long-range Gaussian interactions, respectively. We sample the
distribution of residence time $t$ in suitably defined metastable valleys
entered at age $t_w$, finding that the deviations from $t/t_w$ scaling are
small and decrease systematically as the system size grows and/or the
temperature decreases. Finally, the connection between this behavior and the
scaling of the correlation function itself is discussed.
