Mercurial > octave-libgccjit
view scripts/statistics/distributions/poissrnd.m @ 11523:fd0a3ac60b0e
update copyright notices
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 14 Jan 2011 05:47:45 -0500 |
parents | 1740012184f9 |
children | c792872f8942 |
line wrap: on
line source
## Copyright (C) 1995-2011 Kurt Hornik ## ## 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} {} poissrnd (@var{lambda}, @var{r}, @var{c}) ## Return an @var{r} by @var{c} matrix of random samples from the ## Poisson distribution with parameter @var{lambda}, which must be a ## scalar or of size @var{r} by @var{c}. ## ## If @var{r} and @var{c} are omitted, the size of the result matrix is ## the size of @var{lambda}. ## @end deftypefn ## Author: KH <Kurt.Hornik@wu-wien.ac.at> ## Description: Random deviates from the Poisson distribution function rnd = poissrnd (lambda, r, c) if (nargin == 3) if (! (isscalar (r) && (r > 0) && (r == round (r)))) error ("poissrnd: R must be a positive integer"); endif if (! (isscalar (c) && (c > 0) && (c == round (c)))) error ("poissrnd: C must be a positive integer"); endif sz = [r, c]; if (any (size (lambda) != 1) && ((length (size (lambda)) != length (sz)) || any (size (lambda) != sz))) error ("poissrnd: LAMBDA must be scalar or of size [R, C]"); endif elseif (nargin == 2) if (isscalar (r) && (r > 0)) sz = [r, r]; elseif (isvector(r) && all (r > 0)) sz = r(:)'; else error ("poissrnd: R must be a positive integer or vector"); endif if (any (size (lambda) != 1) && ((length (size (lambda)) != length (sz)) || any (size (lambda) != sz))) error ("poissrnd: LAMBDA must be scalar or of size sz"); endif elseif (nargin == 1) sz = size (lambda); else print_usage (); endif if (isscalar (lambda)) if (!(lambda >= 0) || !(lambda < Inf)) rnd = NaN (sz); elseif (lambda > 0 && lambda < Inf) rnd = randp(lambda, sz); else rnd = zeros (sz); endif else rnd = zeros (sz); k = find (!(lambda >= 0) | !(lambda < Inf)); if (any (k)) rnd(k) = NaN; endif k = find ((lambda > 0) & (lambda < Inf)); if (any (k)) rnd(k) = randp(lambda(k), size(k)); endif endif endfunction