### Function Reference: `bisainv`

statistics: x = bisainv (p, beta, gamma)

Inverse of the Birnbaum-Saunders cumulative distribution function (iCDF).

For each element of p, compute the quantile (the inverse of the CDF) of the Birnbaum-Saunders distribution with scale parameter beta and shape parameter gamma. The size of x is the common size of p, beta, and gamma. A scalar input functions as a constant matrix of the same size as the other inputs.

Further information about the Birnbaum-Saunders distribution can be found at https://en.wikipedia.org/wiki/Birnbaum%E2%80%93Saunders_distribution

See also: bisainv, bisapdf, bisarnd, bisafit, bisalike, bisastat

Source Code: bisainv

### Example: 1

     ## Plot various iCDFs from the Birnbaum-Saunders distribution
     p = 0.001:0.001:0.999;
     x1 = bisainv (p, 1, 0.5);
     x2 = bisainv (p, 1, 1);
     x3 = bisainv (p, 1, 2);
     x4 = bisainv (p, 1, 5);
     x5 = bisainv (p, 1, 10);
     plot (p, x1, "-b", p, x2, "-g", p, x3, "-r", p, x4, "-c", p, x5, "-m")
     grid on
     ylim ([0, 10])
     legend ({"β = 1, γ = 0.5", "β = 1, γ = 1", "β = 1, γ = 2", ...
              "β = 1, γ = 5", "β = 1, γ = 10"}, "location", "northwest")
     title ("Birnbaum-Saunders iCDF")
     xlabel ("probability")
     ylabel ("values in x")

### Example: 2

     ## Plot various iCDFs from the Birnbaum-Saunders distribution
     p = 0.001:0.001:0.999;
     x1 = bisainv (p, 1, 0.3);
     x2 = bisainv (p, 2, 0.3);
     x3 = bisainv (p, 1, 0.5);
     x4 = bisainv (p, 3, 0.5);
     x5 = bisainv (p, 5, 0.5);
     plot (p, x1, "-b", p, x2, "-g", p, x3, "-r", p, x4, "-c", p, x5, "-m")
     grid on
     ylim ([0, 10])
     legend ({"β = 1, γ = 0.3", "β = 2, γ = 0.3", "β = 1, γ = 0.5", ...
              "β = 3, γ = 0.5", "β = 5, γ = 0.5"}, "location", "northwest")
     title ("Birnbaum-Saunders iCDF")
     xlabel ("probability")
     ylabel ("values in x")