It would be great if this library could extend some/all of the functions to work with Decimal data types for cases where we value the additional accuracy over improved performance. Currently I can cast return values from your library to a Decimal with |> Decimal.from_float(), but by that point some accuracy has already been lost due to floating point calculations:
iex(1)> Numerix.Statistics.mean([0.1, 0.1, 0.1])
0.10000000000000002