/usr/share/octave/packages/statistics-1.3.0/cdf.m is in octave-statistics 1.3.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 | ## Copyright (C) 2013 Pantxo Diribarne
##
## This program is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 3 of the License, or
## (at your option) any later version.
##
## This program is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
## GNU General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with this program. If not, see <http://www.gnu.org/licenses/>.
## -*- texinfo -*-
## @deftypefn {Function File} {@var{retval} =} cdf (@var{name}, @var{X}, @dots{})
## Return cumulative density function of @var{name} function for value
## @var{x}.
## This is a wrapper around various @var{name}cdf and @var{name}_cdf
## functions. See the individual functions help to learn the signification of
## the arguments after @var{x}. Supported functions and corresponding number of
## additional arguments are:
##
## @multitable @columnfractions 0.02 0.3 0.45 0.2
## @headitem @tab function @tab alternative @tab args
## @item @tab "beta" @tab "beta" @tab 2
## @item @tab "bino" @tab "binomial" @tab 2
## @item @tab "cauchy" @tab @tab 2
## @item @tab "chi2" @tab "chisquare" @tab 1
## @item @tab "discrete" @tab @tab 2
## @item @tab "exp" @tab "exponential" @tab 1
## @item @tab "f" @tab @tab 2
## @item @tab "gam" @tab "gamma" @tab 2
## @item @tab "geo" @tab "geometric" @tab 1
## @item @tab "gev" @tab "generalized extreme value" @tab 3
## @item @tab "hyge" @tab "hypergeometric" @tab 3
## @item @tab "kolmogorov_smirnov" @tab @tab 1
## @item @tab "laplace" @tab @tab 0
## @item @tab "logistic" @tab @tab 0
## @item @tab "logn" @tab "lognormal" @tab 2
## @item @tab "norm" @tab "normal" @tab 2
## @item @tab "poiss" @tab "poisson" @tab 1
## @item @tab "rayl" @tab "rayleigh" @tab 1
## @item @tab "t" @tab @tab 1
## @item @tab "unif" @tab "uniform" @tab 2
## @item @tab "wbl" @tab "weibull" @tab 2
## @end multitable
##
## @seealso{betacdf, binocdf, cauchy_cdf, chi2cdf, discrete_cdf,
## expcdf, fcdf, gamcdf, geocdf, gevcdf, hygecdf,
## kolmogorov_smirnov_cdf, laplace_cdf, logistic_cdf, logncdf,
## normcdf, poisscdf, raylcdf, tcdf, unifcdf, wblcdf}
## @end deftypefn
function [retval] = cdf (varargin)
## implemented functions
persistent allcdf = {{"beta", "beta"}, @betacdf, 2, ...
{"bino", "binomial"}, @binocdf, 2, ...
{"cauchy"}, @cauchy_cdf, 2, ...
{"chi2", "chisquare"}, @chi2cdf, 1, ...
{"discrete"}, @discrete_cdf, 2, ...
{"exp", "exponential"}, @expcdf, 1, ...
{"f"}, @fcdf, 2, ...
{"gam", "gamma"}, @gamcdf, 2, ...
{"geo", "geometric"}, @geocdf, 1, ...
{"gev", "generalized extreme value"}, @gevcdf, 3, ...
{"hyge", "hypergeometric"}, @hygecdf, 3, ...
{"kolmogorov_smirnov"}, @kolmogorov_smirnov_cdf, 1, ...
{"laplace"}, @laplace_cdf, 0, ...
{"logistic"}, @logistic_cdf, 0, ... # ML has 2 args here
{"logn", "lognormal"}, @logncdf, 2, ...
{"norm", "normal"}, @normcdf, 2, ...
{"poiss", "poisson"}, @poisscdf, 1, ...
{"rayl", "rayleigh"}, @raylcdf, 1, ...
{"t"}, @tcdf, 1, ...
{"unif", "uniform"}, @unifcdf, 2, ...
{"wbl", "weibull"}, @wblcdf, 2};
if (numel (varargin) < 2 || ! ischar (varargin{1}))
print_usage ();
endif
name = varargin{1};
x = varargin{2};
varargin(1:2) = [];
nargs = numel (varargin);
cdfnames = allcdf(1:3:end);
cdfhdl = allcdf(2:3:end);
cdfargs = allcdf(3:3:end);
idx = cellfun (@(x) any (strcmpi (name, x)), cdfnames);
if (any (idx))
if (nargs == cdfargs{idx})
retval = feval (cdfhdl{idx}, x, varargin{:});
else
error ("cdf: %s requires %d arguments", name, cdfargs{idx})
endif
else
error ("cdf: %s not implemented", name);
endif
endfunction
%!test
%! assert(cdf ('norm', 1, 0, 1), normcdf (1, 0, 1))
|