import mir.test: shouldApprox;

logisticInvCDF(0.0).shouldApprox == -double.infinity;
logisticInvCDF(0.25).shouldApprox == -1.098612;
logisticInvCDF(0.5).shouldApprox == 0.0;
logisticInvCDF(0.75).shouldApprox == 1.098612;
logisticInvCDF(1.0).shouldApprox == double.infinity;

// Can also provide location/scale parameters
logisticInvCDF(0.2, 2, 3).shouldApprox == -2.158883;
logisticInvCDF(0.4, 2, 3).shouldApprox == 0.7836047;
logisticInvCDF(0.6, 2, 3).shouldApprox == 3.216395;
logisticInvCDF(0.8, 2, 3).shouldApprox == 6.158883;