Mercurial > octave-antonio
annotate scripts/statistics/distributions/wblrnd.m @ 11587:c792872f8942
all script files: untabify and strip trailing whitespace
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Thu, 20 Jan 2011 17:35:29 -0500 |
parents | fd0a3ac60b0e |
children | 19b9f17d22af |
rev | line source |
---|---|
11523 | 1 ## Copyright (C) 1995-2011 Kurt Hornik |
5693 | 2 ## |
3 ## This file is part of Octave. | |
4 ## | |
5 ## Octave is free software; you can redistribute it and/or modify it | |
6 ## under the terms of the GNU General Public License as published by | |
7016 | 7 ## the Free Software Foundation; either version 3 of the License, or (at |
8 ## your option) any later version. | |
5693 | 9 ## |
10 ## Octave is distributed in the hope that it will be useful, but | |
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of | |
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
13 ## General Public License for more details. | |
14 ## | |
15 ## You should have received a copy of the GNU General Public License | |
7016 | 16 ## along with Octave; see the file COPYING. If not, see |
17 ## <http://www.gnu.org/licenses/>. | |
5693 | 18 |
19 ## -*- texinfo -*- | |
10793
be55736a0783
Grammarcheck the documentation from m-files.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
20 ## @deftypefn {Function File} {} wblrnd (@var{scale}, @var{shape}, @var{r}, @var{c}) |
5693 | 21 ## @deftypefnx {Function File} {} wblrnd (@var{scale}, @var{shape}, @var{sz}) |
22 ## Return an @var{r} by @var{c} matrix of random samples from the | |
23 ## Weibull distribution with parameters @var{scale} and @var{shape} | |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
7017
diff
changeset
|
24 ## which must be scalar or of size @var{r} by @var{c}. Or if @var{sz} |
5693 | 25 ## is a vector return a matrix of size @var{sz}. |
26 ## | |
27 ## If @var{r} and @var{c} are omitted, the size of the result matrix is | |
28 ## the common size of @var{alpha} and @var{sigma}. | |
29 ## @end deftypefn | |
30 | |
31 ## Author: KH <Kurt.Hornik@wu-wien.ac.at> | |
32 ## Description: Random deviates from the Weibull distribution | |
33 | |
34 function rnd = wblrnd (scale, shape, r, c) | |
35 | |
36 if (nargin > 1) | |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
37 if (!isscalar(scale) || !isscalar(shape)) |
11128
9cb5c0b7b43b
Fix reversed documentation for shape and scale parameters in Weibull family of functions.
Rik <octave@nomad.inbox5.com>
parents:
11094
diff
changeset
|
38 [retval, scale, shape] = common_size (scale, shape); |
5693 | 39 if (retval > 0) |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
11149
diff
changeset
|
40 error ("wblrnd: SCALE and SHAPE must be of common size or scalar"); |
5693 | 41 endif |
42 endif | |
43 endif | |
44 | |
45 if (nargin == 4) | |
46 if (! (isscalar (r) && (r > 0) && (r == round (r)))) | |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
11149
diff
changeset
|
47 error ("wblrnd: R must be a positive integer"); |
5693 | 48 endif |
49 if (! (isscalar (c) && (c > 0) && (c == round (c)))) | |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
11149
diff
changeset
|
50 error ("wblrnd: C must be a positive integer"); |
5693 | 51 endif |
52 sz = [r, c]; | |
53 | |
11149
fe3c3dfc07eb
style fix: break lines before && and ||, not after
John W. Eaton <jwe@octave.org>
parents:
11128
diff
changeset
|
54 if (any (size (scale) != 1) |
fe3c3dfc07eb
style fix: break lines before && and ||, not after
John W. Eaton <jwe@octave.org>
parents:
11128
diff
changeset
|
55 && ((length (size (scale)) != length (sz)) |
fe3c3dfc07eb
style fix: break lines before && and ||, not after
John W. Eaton <jwe@octave.org>
parents:
11128
diff
changeset
|
56 || any (size (scale) != sz))) |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
11149
diff
changeset
|
57 error ("wblrnd: SCALE and SHAPE must be scalar or of size [R, C]"); |
5693 | 58 endif |
59 elseif (nargin == 3) | |
60 if (isscalar (r) && (r > 0)) | |
61 sz = [r, r]; | |
62 elseif (isvector(r) && all (r > 0)) | |
63 sz = r(:)'; | |
64 else | |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
11149
diff
changeset
|
65 error ("wblrnd: R must be a positive integer or vector"); |
5693 | 66 endif |
67 | |
11149
fe3c3dfc07eb
style fix: break lines before && and ||, not after
John W. Eaton <jwe@octave.org>
parents:
11128
diff
changeset
|
68 if (any (size (scale) != 1) |
fe3c3dfc07eb
style fix: break lines before && and ||, not after
John W. Eaton <jwe@octave.org>
parents:
11128
diff
changeset
|
69 && ((length (size (scale)) != length (sz)) |
fe3c3dfc07eb
style fix: break lines before && and ||, not after
John W. Eaton <jwe@octave.org>
parents:
11128
diff
changeset
|
70 || any (size (scale) != sz))) |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
11149
diff
changeset
|
71 error ("wblrnd: SCALE and SHAPE must be scalar or of size SZ"); |
5693 | 72 endif |
73 elseif (nargin == 2) | |
11128
9cb5c0b7b43b
Fix reversed documentation for shape and scale parameters in Weibull family of functions.
Rik <octave@nomad.inbox5.com>
parents:
11094
diff
changeset
|
74 sz = size(scale); |
5693 | 75 else |
6046 | 76 print_usage (); |
5693 | 77 endif |
78 | |
11128
9cb5c0b7b43b
Fix reversed documentation for shape and scale parameters in Weibull family of functions.
Rik <octave@nomad.inbox5.com>
parents:
11094
diff
changeset
|
79 if (isscalar (scale) && isscalar (shape)) |
9cb5c0b7b43b
Fix reversed documentation for shape and scale parameters in Weibull family of functions.
Rik <octave@nomad.inbox5.com>
parents:
11094
diff
changeset
|
80 if (scale > 0 && scale < Inf && shape > 0 && shape < Inf) |
6350 | 81 rnd = scale .* rande(sz) .^ (1./shape); |
5693 | 82 else |
10525
3306cfcb856e
Replace constructs like "NaN * one()" with "NaN()" and "Inf * ones ()" with "Inf()"
David Bateman <dbateman@free.fr>
parents:
9245
diff
changeset
|
83 rnd = NaN (sz); |
5693 | 84 endif |
85 else | |
6350 | 86 rnd = scale .* rande(sz) .^ (1./shape); |
11128
9cb5c0b7b43b
Fix reversed documentation for shape and scale parameters in Weibull family of functions.
Rik <octave@nomad.inbox5.com>
parents:
11094
diff
changeset
|
87 k = find ((scale <= 0) | (scale == Inf) | ((shape <= 0) & (shape == Inf))); |
6350 | 88 if (any(k)) |
89 rnd(k) = NaN; | |
5693 | 90 endif |
91 endif | |
92 | |
93 endfunction | |
94 |