We study the rate of convergence w.r.t.~a Wasserstein type distance for random walk approximations of mean field BSDEs. Our method does not use the particle method but instead a freezing technique. We extend results by Briand, Ch. Geiss, S. Geiss, and Labart [Bernoulli, 27(2) 2021] about the rate of convergence of a Donsker-type theorem for BSDEs from the classical setting to the mean field setting. In this connection the mean field setting leads to new phenomena and requires new techniques that should be of independent interest: The Hölder continuous terminal condition causes a singularity in time of the generator when seen as a generator in the non-mean field setting. To handle this singularity we introduce a concept of modified Hölder continuity by which we are able to achieve, up to a logarithmic term, the same polynomial approximation rates as in the classical non-mean field setting (in fact, already when approximating the Brownian motion itself a logarithmic term is necessary). Moreover, the exploited freezing technique of the mean field terms yields to the problem to handle the quantitative behavior of several different generators. Using BMO-techniques we obtain the rate of convergence for the integrated gradient process in the scale of Lorentz (type) spaces of exponential type.