  We show explicitly that four neutrino models of the 2+2 variety still provide
an acceptable global fit to the solar, atmospheric and LSND neutrino data. The
goodness of fit, defined in the usual way, is found to be 0.26 for the simplest
such model. That is, we find that there is a 26% probability of obtaining a
worse global fit to the neutrino data We also make some specific comments on
the paper, ``Ruling out four-neutrino oscillation interpretations of the LSND
anomaly'' [hep-ph/0207157], and explain why they reached drastically different
conclusions.
