The function calculates the complete combinatorial of the supplied vectors
using a loop implementation by taking the difference between every
combination of the elements in the vectors. The input vectors must
be numeric, but can be of different lengths.

## Arguments

- x
An input vector

- y
An input vector

## Value

A numeric test value

## References

Poe G. L., Giraud, K. L. & Loomis, J. B., 2005, Computational methods for
measuring the difference of empirical distributions, American Journal of
Agricultural Economics, 87:353-365

## Examples

```
x <- qnorm(runif(100), mean = -0.5, sd = 1)
y <- qnorm(runif(100), mean = 1.5, sd = 2)
poe_test(x, y)
#> [1] 0.8223
```