A. Neumaier said:
I gave the simplest constructive choice that leads to a valid solution, namely setting to zero enough higher order parameters in a given renormalization fixing scheme
The space of choices of renormalization parameters at each order is not a vector space, but an affine space. There is no invariant meaning of "setting to zero" these parameters, unless one has already
chosen an origin in these affine spaces. The latter may be addressed as a choice of renormalization scheme, but this just gives another name to the choice to be made, it still does not give a canonical choice.
You know this, but here is pointers to the details for those readers who have not see this:
In the original
Epstein-Glaser 73 the choice at order ##\omega## happens on p. 27, where it says "we choose a fixed auxiliary function ##w \in \mathcal{S}(\mathbb{R}^n)## such that...". With the choice of this function they build one solution to the renormalization problem at this order (for them a splitting of distributions) which they call ##(T^+, T^-)##. With this "origin" chosen, every other solution of the renormalization at that order is labeled by a vector space of renormalization constants ##c_\alpha## (on their p. 28, after "The most general solution"). It might superficially seem the as if we could renormalize canonically by declaring "choose all ##c_\alpha## to be zero". But this is an illusion, the choice is now in the "scheme" ##w## relative to which the ##c_\alpha## are given.
In the modern reformulation of Epstein-Glaser's work in terms of extensions of distributions in
Brunetti-Fredenhagen 00 the analogous step happens on p. 24 in or below equation (38), where at order ##\omega## bump functions ##\mathfrak{w}_\alpha## are chosen. The theorem 5.3 below that states then that with this choice, the space of renormalization constants at that order is given by coefficients
relative to these choices ##\mathfrak{w}_\alpha##.
One may succintly summarize this statement by saying that the space of renormalization parameters at each order, while not having a preferred element (in particular not being a vector space with a zero-element singled out) is a torsor over a vector space, meaning that after anyone point is chosen, then the remaining points form a vector space relative to this point. That more succinct formulation of theorem 5.3 in
Brunetti-Fredenhagen 00 is made for instance as corollary 2.6 on p.5 of
Bahns-Wrochna 12.
Hence for a general Lagrangian there is no formula for choosing the renormalization parameters at each order. It is in very special situations only that we may give a formula for choosing the infinitely many renormalization parameters. Three prominent such situations are the following:
1) if the theory is "renormalizable" in that it so happens that after some finite order the space of choices of parameters contain a unique single point. In that case we may make a finite number of choices and then the remaining choices are fixed.
2) If we assume the existence of a "UV-critical hypersurface" (e.g.
Nink-Reuter 12, p. 2), which comes down to postulating a finite dimensional submanifold in the infinite dimensional space of renormalization parameters and postulating/assuming that we make a choice on this submanifold. Major extra assumptions here. If they indeed happen to be met, then the space of choices is drastically shrunk.
3) We assume a UV-completion by a string perturbation series. This only works for field theories which are not in the "swampland" (
Vafa 05). It transforms the space of choices of renormalization parameters into the space of choices of full 2d SCFTS of central charge 15, the latter also known as the "perturbative landscape". Even though this space received a lot of press, it seems that way too little is known about it to say much at the moment. But that's another discussion.
There might be more, but the above three seem to be the most prominent ones.