© Owned by the authors, published by EDP Sciences, 2014.Dyson-Schwinger equations are the most common tool for the determination of the correlation functions of Landau gauge Yang-Mills theory in the continuum, in particular in the infrared regime. We shall argue that the use of Callan-Symanzik renormalization group equations has distinctive advantages over the Dyson-Schwinger equations, in particular for the vertex functions. We present a generalization of the infrared safe renormalization scheme proposed by Tissier and Wschebor in 2011. The comparison with the existing lattice data for the gluon and ghost propagators can be used to determine the most appropriate renormalization scheme.