Atomistic-to-continuum (a/c) coupling methods are a class of computational multiscale schemes that combine the accuracy of atomistic models with the efficiency of continuum elasticity. They are increasingly being utilized in materials science to study the fundamental mechanisms of material failure such as crack propagation and plasticity, which are governed by the interaction between crystal defects and long-range elastic fields. In the construction of a/c coupling methods, various approximation errors are committed. A rigorous numerical analysis approach that classifies and quantifies these errors can give confidence in the simulation results, as well as enable optimization of the numerical methods for accuracy and computational cost. In this article, we present such a numerical analysis framework, which is inspired by recent research activity.