The vector autoregressive (VAR) model has been widely used for modeling temporal dependence in a multivariate time series. For large (and even moderate) dimensions, the number of the AR coefficients can be prohibitively large, resulting in noisy estimates, unstable predictions, and difficult-to-interpret temporal dependence. To overcome such drawbacks, we propose a two-stage approach for fitting sparse VAR (sVAR) models in which many of the AR coefficients are zero. The first stage selects nonzero AR coefficients based on an estimate of the partial spectral coherence (PSC) together with the use of BIC. The PSC is useful for quantifying the conditional relationship between marginal series in a multivariate process. A refinement second stage is then applied to further reduce the number of parameters. The performance of this two-stage approach is illustrated with simulation and real data examples. Supplementary materials for this article are available online.