Mercurial > octave-antonio

comparison scripts/statistics/tests/anova.m @ 11472:1740012184f9

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
author Rik <octave@nomad.inbox5.com>
date Sun, 09 Jan 2011 21:33:04 -0800
parents 16f53d29049f
children fd0a3ac60b0e
comparison
equal deleted inserted replaced
11471:994e2a93a8e2 11472:1740012184f9
49 49
50 if ((nargin < 1) || (nargin > 2)) 50 if ((nargin < 1) || (nargin > 2))
51 print_usage (); 51 print_usage ();
52 elseif (nargin == 1) 52 elseif (nargin == 1)
53 if (isvector (y)) 53 if (isvector (y))
54 error ("anova: for `anova (y)', y must not be a vector"); 54 error ("anova: for `anova (Y)', Y must not be a vector");
55 endif 55 endif
56 [group_count, k] = size (y); 56 [group_count, k] = size (y);
57 n = group_count * k; 57 n = group_count * k;
58 group_mean = mean (y); 58 group_mean = mean (y);
59 else 59 else
60 if (! isvector (y)) 60 if (! isvector (y))
61 error ("anova: for `anova (y, g)', y must be a vector"); 61 error ("anova: for `anova (Y, G)', Y must be a vector");
62 endif 62 endif
63 n = length (y); 63 n = length (y);
64 if (! isvector (g) || (length (g) != n)) 64 if (! isvector (g) || (length (g) != n))
65 error ("anova: g must be a vector of the same length as y"); 65 error ("anova: G must be a vector of the same length as Y");
66 endif 66 endif
67 s = sort (g); 67 s = sort (g);
68 i = find (s (2 : n) > s(1 : (n-1))); 68 i = find (s (2 : n) > s(1 : (n-1)));
69 k = length (i) + 1; 69 k = length (i) + 1;
70 if (k == 1) 70 if (k == 1)