The $\Lambda$CDM simulations of Keller and Wadsley do not account for
  the MOND mass-discrepancy-acceleration relation
  Keller and Wadsley (2016) have smugly suggested, recently, that the end of
MOND may be in view. This is based on their claim that their highly-restricted
sample of $\Lambda$CDM-simulated galaxies are "consistent" with the observed
MOND mass-discrepancy-acceleration relation (MDAR), in particular, with its
recent update by McGaugh et al. (2016), based on the SPARC sample. From this
they extrapolate to "$\Lambda$CDM is fully consistent" with the MDAR. I explain
why these simulated galaxies do not show that $\Lambda$CDM accounts for the
MDAR. a. Their sample of simulated galaxies contains only 18 high-mass
galaxies, within a narrow range of one order of magnitude in baryonic mass, at
the very high end of the observed, SPARC sample, which spans 4.5 orders of
magnitude in mass. More importantly, the simulated sample has none of the
low-mass, low-acceleration galaxies -- abundant in SPARC -- which encapsulate
the crux and the nontrivial aspects of the predicted and observed MDAR. The
low-acceleration part of the simulated MDAR is achieved, rather trivially, from
the flattish-velocity-curve regions of the simulated high-mass galaxies. b.
Half of the simulated galaxies have "wrong" rotation curves that differ greatly
from any observed ones. This, does not prevent these wrong galaxies from lying
on the observed MDAR (for trivial reasons, again). They, in fact, define the
high-acceleration branch of the simulated MDAR. c. To boot, even if
$\Lambda$CDM were made "consistent" with the MDAR through the elaborate
adjustments that the simulations allow, this would not obviate MOND, which
predicts much more than the MDAR.
