  We compute the four-loop contributions to the $\beta$-function and the
anomalous dimension of the field for the $O(N)$-invariant $N$-vector model.
These results are used to compute the second analytic corrections to the
correlation length and the general spin-$n$ susceptibility.
