Mercurial > octave
comparison scripts/statistics/base/statistics.m @ 3426:f8dde1807dee
[project @ 2000-01-13 08:40:00 by jwe]
author | jwe |
---|---|
date | Thu, 13 Jan 2000 08:40:53 +0000 |
parents | eb27ea9b7ff8 |
children | 71d2e09c15a2 |
comparison
equal
deleted
inserted
replaced
3425:8625164a0a39 | 3426:f8dde1807dee |
---|---|
1 ## Copyright (C) 1995, 1996, 1997 Kurt Hornik | 1 ## Copyright (C) 1995, 1996, 1997 Kurt Hornik |
2 ## | 2 ## |
3 ## This program is free software; you can redistribute it and/or modify | 3 ## This program is free software; you can redistribute it and/or modify |
4 ## it under the terms of the GNU General Public License as published by | 4 ## it under the terms of the GNU General Public License as published by |
5 ## the Free Software Foundation; either version 2, or (at your option) | 5 ## the Free Software Foundation; either version 2, or (at your option) |
6 ## any later version. | 6 ## any later version. |
7 ## | 7 ## |
8 ## This program is distributed in the hope that it will be useful, but | 8 ## This program is distributed in the hope that it will be useful, but |
9 ## WITHOUT ANY WARRANTY; without even the implied warranty of | 9 ## WITHOUT ANY WARRANTY; without even the implied warranty of |
10 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | 10 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
11 ## General Public License for more details. | 11 ## General Public License for more details. |
12 ## | 12 ## |
13 ## You should have received a copy of the GNU General Public License | 13 ## You should have received a copy of the GNU General Public License |
14 ## along with this file. If not, write to the Free Software Foundation, | 14 ## along with this file. If not, write to the Free Software Foundation, |
15 ## 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. | 15 ## 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. |
16 | 16 |
17 ## usage: S = statistics (X) | 17 ## usage: S = statistics (X) |
19 ## If X is a matrix, return a matrix S with the min, first quartile, | 19 ## If X is a matrix, return a matrix S with the min, first quartile, |
20 ## median, third quartile, max, mean, std, skewness and kurtosis of the | 20 ## median, third quartile, max, mean, std, skewness and kurtosis of the |
21 ## columns of X as its rows. | 21 ## columns of X as its rows. |
22 ## | 22 ## |
23 ## If X is a vector, treat it as a column vector. | 23 ## If X is a vector, treat it as a column vector. |
24 | 24 |
25 ## Author: KH <Kurt.Hornik@ci.tuwien.ac.at> | 25 ## Author: KH <Kurt.Hornik@ci.tuwien.ac.at> |
26 ## Description: Compute basic statistics | 26 ## Description: Compute basic statistics |
27 | 27 |
28 function S = statistics (X) | 28 function S = statistics (X) |
29 | 29 |
30 if (nargin != 1) | 30 if (nargin != 1) |
31 usage ("S = statistics (X)"); | 31 usage ("S = statistics (X)"); |
32 endif | 32 endif |
33 | 33 |
34 if (prod (size (X)) > 1) | 34 if (prod (size (X)) > 1) |
35 if (is_vector (X)) | 35 if (is_vector (X)) |
36 X = reshape (X, length (X), 1); | 36 X = reshape (X, length (X), 1); |
37 endif | 37 endif |
38 for k=1:columns(X) | 38 for k=1:columns(X) |
39 S(:,k) = [(min (X(:,k))); | 39 S(:,k) = [(min (X(:,k))); |
40 (empirical_inv ([0.25;0.5;0.75], X(:,k))); | 40 (empirical_inv ([0.25;0.5;0.75], X(:,k))); |
41 (max (X(:,k))); | 41 (max (X(:,k))); |
42 (mean (X(:,k))); | 42 (mean (X(:,k))); |
43 (std (X(:,k))); | 43 (std (X(:,k))); |
44 (skewness (X(:,k))); | 44 (skewness (X(:,k))); |
45 (kurtosis (X(:,k)))]; | 45 (kurtosis (X(:,k)))]; |
46 endfor | 46 endfor |
47 else | 47 else |
48 error ("statistics: invalid argument"); | 48 error ("statistics: invalid argument"); |
49 endif | 49 endif |
50 | 50 |
51 endfunction | 51 endfunction |