We study the low-energy behavior of the vertex function of a single Anderson impurity away from half filling for finite magnetic fields, using the Ward identities with careful consideration of the antisymmetry and analytic properties. The asymptotic form of the vertex function ${\mathrm{\ensuremath{\Gamma}}}_{\ensuremath{\sigma}{\ensuremath{\sigma}}^{\ensuremath{'}};{\ensuremath{\sigma}}^{\ensuremath{'}}\ensuremath{\sigma}}^{}(i\ensuremath{\omega},i{\ensuremath{\omega}}^{\ensuremath{'}};i{\ensuremath{\omega}}^{\ensuremath{'}},i\ensuremath{\omega})$ is determined up to terms of linear order with respect to the two frequencies $\ensuremath{\omega}$ and ${\ensuremath{\omega}}^{\ensuremath{'}}$, as well as the ${\ensuremath{\omega}}^{2}$ contribution for antiparallel spins ${\ensuremath{\sigma}}^{\ensuremath{'}}\ensuremath{\ne}\ensuremath{\sigma}$ at ${\ensuremath{\omega}}^{\ensuremath{'}}=0$. From these results, we also obtain a series of the Fermi-liquid relations beyond those of Yamada-Yosida [Prog. Theor. Phys. 54, 316 (1975)]. The ${\ensuremath{\omega}}^{2}$ real part of the self-energy ${\mathrm{\ensuremath{\Sigma}}}_{\ensuremath{\sigma}}^{}(i\ensuremath{\omega})$ is shown to be expressed in terms of the double derivative ${\ensuremath{\partial}}^{2}{\mathrm{\ensuremath{\Sigma}}}_{\ensuremath{\sigma}}^{}(0)/\ensuremath{\partial}{\ensuremath{\epsilon}}_{d\ensuremath{\sigma}}^{2}$ with respect to the impurity energy level ${\ensuremath{\epsilon}}_{d\ensuremath{\sigma}}^{}$, and agrees with the formula obtained recently by Filippone, Moca, von Delft, and Mora (FMvDM) in the Nozi\`eres phenomenological Fermi-liquid theory [Phys. Rev. B 95, 165404 (2017)]. We also calculate the ${T}^{2}$ correction of the self-energy and find that the real part can be expressed in terms of the three-body correlation function $\ensuremath{\partial}{\ensuremath{\chi}}_{\ensuremath{\uparrow}\ensuremath{\downarrow}}/\ensuremath{\partial}{\ensuremath{\epsilon}}_{d,\ensuremath{-}\ensuremath{\sigma}}^{}$, where ${\ensuremath{\chi}}_{\ensuremath{\uparrow}\ensuremath{\downarrow}}$ is the static susceptibility between antiparallel spins. We also provide an alternative derivation of the asymptotic form of the vertex function. Specifically, we calculate the skeleton diagrams for the vertex function ${\mathrm{\ensuremath{\Gamma}}}_{\ensuremath{\sigma}\ensuremath{\sigma};\ensuremath{\sigma}\ensuremath{\sigma}}^{}(i\ensuremath{\omega},0;0,i\ensuremath{\omega})$ for parallel spins up to order ${U}^{4}$ in the Coulomb repulsion $U$. It directly clarifies the fact that the analytic components of order $\ensuremath{\omega}$ vanish as a result of the cancellation of four related Feynman diagrams, which are related to each other through the antisymmetry operation.