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")
plotted figure

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")
plotted figure