Mercurial > octave-nkf
diff scripts/statistics/tests/kolmogorov_smirnov_test.m @ 3456:434790acb067
[project @ 2000-01-19 06:58:51 by jwe]
author | jwe |
---|---|
date | Wed, 19 Jan 2000 06:59:23 +0000 |
parents | d8b731d3f7a3 |
children | e031284eea27 |
line wrap: on
line diff
--- a/scripts/statistics/tests/kolmogorov_smirnov_test.m Tue Jan 18 19:57:13 2000 +0000 +++ b/scripts/statistics/tests/kolmogorov_smirnov_test.m Wed Jan 19 06:59:23 2000 +0000 @@ -43,19 +43,17 @@ ## If no output argument is given, the p-value is displayed. ## @end deftypefn -## Author: KH <Kurt.Hornik@ci.tuwien.ac.at> -## Description: One-sample Kolmogorov-Smirnov test +## Author: KH <Kurt.Hornik@ci.tuwien.ac.at> +## Description: One-sample Kolmogorov-Smirnov test function [pval, ks] = kolmogorov_smirnov_test (x, dist, ...) if (nargin < 2) - error (sprintf (["usage:\n\t", ... - "[pval, ks] = ", ... - "kolmogorov_smirnov_test (x, dist, [, params] [, alt])"])); + usage ("[pval, ks] = kolmogorov_smirnov_test (x, dist,, params, alt)"); endif if (! is_vector (x)) - error ("kolmogorov_smirnov_test: x must be a vector."); + error ("kolmogorov_smirnov_test: x must be a vector."); endif n = length (x); @@ -80,21 +78,20 @@ endif if (strcmp (alt, "!=") || strcmp (alt, "<>")) - ks = sqrt(n) * max(max([abs(z - (0:(n-1))/n); abs(z - (1:n)/n)])); + ks = sqrt (n) * max (max ([abs(z - (0:(n-1))/n); abs(z - (1:n)/n)])); pval = 1 - kolmogorov_smirnov_cdf (ks); elseif (strcmp (alt, ">")) - ks = sqrt(n) * max (max ([z - (0:(n-1))/n; z - (1:n)/n])); - pval = exp(- 2 * ks^2); + ks = sqrt (n) * max (max ([z - (0:(n-1))/n; z - (1:n)/n])); + pval = exp (- 2 * ks^2); elseif (strcmp (alt, "<")) - ks = - sqrt(n) * min (min ([z - (0:(n-1))/n; z - (1:n)/n])); - pval = exp(- 2 * ks^2); + ks = - sqrt (n) * min (min ([z - (0:(n-1))/n; z - (1:n)/n])); + pval = exp (- 2 * ks^2); else - error (sprintf (["kolmogorov_smirnov_test: ", ... - "alternative %s not recognized"], alt)); + error ("kolmogorov_smirnov_test: alternative %s not recognized", alt); endif if (nargout == 0) - printf ("pval: %g\n", pval); + printf ("pval: %g\n", pval); endif endfunction