Computations of N-15 NMR chemical shifts in 93 diverse nitrogen-containing compounds representing almost all known classes are performed at the density functional theory (DFT), second-order Moller-Plesset perturbation theory (MP2), and coupled cluster singles and doubles (CCSD) levels using the composite method approximation (CMA) in comparison with experimental results. It is shown that the CMA-DFT and CMA-CCSD methods provided the best performance characterized by a normalized mean absolute error of 1.1-1.3% as compared to 2.3% for the CMA-MP2 results. Taking into account solvent effects within the conductor-like polarizable continuum model decreased the normalized mean absolute error by 0.4% for the CMA-DFT and by 0.2% for the CMA-CCSD calculations.