Short answer
Not all the axioms of a metric are needed (see below). However I can't think of any published examples using a distance which isn't closely based on a metric.
Metrics recap
The axioms of a metric are:
- $d(y,y') \geq 0$ (non-negativity)
- $d(y,y') = 0$ iff $y=y'$ (identifiability)
- $d(y,y') = d(y',y)$ (symmetry)
- $d(y,y') \leq d(y,z) + d(z,y')$ (triangle inequality)
ABC recap
ABC bases acceptance of simulated (parameter, data) pairs $(\theta, y)$ on the distance $d(y,y_0)$ where $y_0$ is the observed data. Small distance values are more likely to be accepted. This answer concentrates on the simple scheme which accepts if $d$ is below a threshold $\epsilon$.
Necessity of the axioms in ABC
- Non-negativity. This axiom is assumed in all implementations of ABC I've seen.
- Identifiability. ABC produces $\theta$ samples from the exact posterior if acceptance only occurs for $y=y_0$. If there are multiple solutions to $d(y,y_0)=0$ this acceptance rule cannot be attained, which seems like a major deficiency. In particular a limit of increasingly stringent acceptance rules would not produce sampling targets converging on the true posterior. So identifiability seems essential.
- Symmetry. This clearly isn't essential as ABC only uses $d(y,y_0)$. We could quite easily take $d(y_0,y) \neq d(y,y_0)$ without affecting the algorithm.
- Triangle inequality. This is also not needed. For example suppose $y \in \mathbb{R}$ and $d_1(y,y')=|y-y'|$ (a metric). Then $d_2(y,y')=(y-y')^2$ breaks the triangle inequality (e.g. consider $y=0, z=1, y'=2$) but produces the same ABC acceptance regions as $d_1$ for appropriately chosen acceptance thresholds.
Some sufficient conditions for ABC distances
In the appendix of a preprint I make a case for some conditions on ABC distance functions which are sufficient for the ABC target distribution to converge to the posterior.
Other conditions could probably be reached by different arguments though!
(My argument is a by-product of some theory I needed for ABC-SMC algorithms.)
The argument is based on the acceptance regions $A_t=\{ y | d(y,y_0) \leq \epsilon_t \}$ for any sequence $\epsilon_1 > \epsilon_2 > \ldots$ converging on zero.
The conditions are (omitting some technicalities):
- The Lebesgue measure $|A_t|$ converges to 0 as $t \to \infty$.
- The sets $A_t$ have bounded eccentricity. Essentially under projection to any lower dimensional space, the Lebesgue measure of $A_t$ still converges to zero.
Condition 1 ensures that the acceptance region shrinks to a vanishingly small region around the observations. Condition 2 avoids a limiting acceptance region in which some but not all components of $y$ match those of $y_0$.