Mercurial > octave-nkf
view scripts/statistics/distributions/unidpdf.m @ 20521:4bb41929286b
Deprecate bitmax.
* NEWS: Announce deprecation and replacement functions.
* scripts/deprecated/bitmax.m: New function. Add deprecation message to
docstring. Add deprecation warning to m-file.
* scripts/deprecated/module.mk: Add deprecated versions to build system.
* libinterp/corefcn/bitfcns.cc: Remove build-in function. Remove function
from docstrings.
* libinterp/corefcn/data.cc: Remove function from docstrings.
* scripts/general/bitcmp.m: Replace bitmax function calls with flintmax.
Remove function from docstrings.
* scripts/general/bitget.m: Replace bitmax function calls with flintmax.
Remove function from docstrings.
* scripts/general/bitset.m: Replace bitmax function calls with flintmax.
Remove function from docstrings.
* scripts/specfun/factor.m: Replace bitmax function calls with flintmax.
Remove function from docstrings.
* scripts/statistics/distributions/unidpdf.m: Remove function from docstrings.
* scripts/general/num2str.m: Replaced function in comment.
* doc/interpreter/doccheck/aspell-octave.en.pws: Remove functions from list.
* doc/interpreter/numbers.txi: Remove function from manual.
author | Kai T. Ohlhus <k.ohlhus@gmail.com> |
---|---|
date | Wed, 09 Sep 2015 18:43:52 +0200 |
parents | 6947b9c82d4b |
children |
line wrap: on
line source
## Copyright (C) 2012 Rik Wehbring ## Copyright (C) 2007-2015 David Bateman ## ## This file is part of Octave. ## ## Octave 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. ## ## Octave 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 Octave; see the file COPYING. If not, see ## <http://www.gnu.org/licenses/>. ## -*- texinfo -*- ## @deftypefn {Function File} {} unidpdf (@var{x}, @var{n}) ## For each element of @var{x}, compute the probability density function (PDF) ## at @var{x} of a discrete uniform distribution which assumes ## the integer values 1--@var{n} with equal probability. ## ## Warning: The underlying implementation uses the double class and will only ## be accurate for @var{n} < @code{flintmax} (@w{@math{2^{53}}}) on ## IEEE 754 compatible systems). ## @end deftypefn function pdf = unidpdf (x, n) if (nargin != 2) print_usage (); endif if (! isscalar (n)) [retval, x, n] = common_size (x, n); if (retval > 0) error ("unidpdf: X and N must be of common size or scalars"); endif endif if (iscomplex (x) || iscomplex (n)) error ("unidpdf: X and N must not be complex"); endif if (isa (x, "single") || isa (n, "single")) pdf = zeros (size (x), "single"); else pdf = zeros (size (x)); endif k = isnan (x) | ! (n > 0 & n == fix (n)); pdf(k) = NaN; k = !k & (x >= 1) & (x <= n) & (x == fix (x)); if (isscalar (n)) pdf(k) = 1 / n; else pdf(k) = 1 ./ n(k); endif endfunction %!shared x,y %! x = [-1 0 1 2 10 11]; %! y = [0 0 0.1 0.1 0.1 0]; %!assert (unidpdf (x, 10*ones (1,6)), y) %!assert (unidpdf (x, 10), y) %!assert (unidpdf (x, 10*[0 NaN 1 1 1 1]), [NaN NaN y(3:6)]) %!assert (unidpdf ([x, NaN], 10), [y, NaN]) ## Test class of input preserved %!assert (unidpdf (single ([x, NaN]), 10), single ([y, NaN])) %!assert (unidpdf ([x, NaN], single (10)), single ([y, NaN])) ## Test input validation %!error unidpdf () %!error unidpdf (1) %!error unidpdf (1,2,3) %!error unidpdf (ones (3), ones (2)) %!error unidpdf (ones (2), ones (3)) %!error unidpdf (i, 2) %!error unidpdf (2, i)