Anion order in perovskite oxynitrides is investigated by a combination of first-principles calculations, cluster expansion method and Monte Carlo simulations.