Mercurial > octave
view scripts/statistics/mode.m @ 24547:fdc9ce839afd
maint: Remove statistics functions which have been shifted to Octave Forge package.
* NEWS: Announce list of functions removed from core.
* scripts/statistics/center.m, scripts/statistics/corr.m,
scripts/statistics/corrcoef.m, scripts/statistics/cov.m,
scripts/statistics/discrete_cdf.m, scripts/statistics/discrete_inv.m,
scripts/statistics/discrete_pdf.m, scripts/statistics/discrete_rnd.m,
scripts/statistics/empirical_cdf.m, scripts/statistics/empirical_inv.m,
scripts/statistics/empirical_pdf.m, scripts/statistics/empirical_rnd.m,
scripts/statistics/histc.m, scripts/statistics/iqr.m,
scripts/statistics/kendall.m, scripts/statistics/kurtosis.m,
scripts/statistics/mean.m, scripts/statistics/meansq.m,
scripts/statistics/median.m, scripts/statistics/mode.m,
scripts/statistics/module.mk, scripts/statistics/moment.m,
scripts/statistics/prctile.m, scripts/statistics/quantile.m,
scripts/statistics/range.m, scripts/statistics/ranks.m,
scripts/statistics/run_count.m, scripts/statistics/runlength.m,
scripts/statistics/skewness.m, scripts/statistics/spearman.m,
scripts/statistics/statistics.m, scripts/statistics/std.m,
scripts/statistics/var.m, scripts/statistics/zscore.m:
Moved from scripts/statistics/base.
* scripts/module.mk: Update build system. Remove @include for directories
base/, tests/, models/, distributions/.
* octave.texi: Update @detailmenu.
* stats.txi: Rewrite statistics chapter with condensed list of functions.
* scripts/statistics/base/module.mk,
scripts/statistics/distributions/module.mk, scripts/statistics/models/module.mk
scripts/statistics/tests/module.mk:
Removed from build system.
* scripts/statistics/base/cloglog.m, scripts/statistics/base/crosstab.m,
scripts/statistics/base/logit.m, scripts/statistics/base/moment.m,
scripts/statistics/base/ppplot.m, scripts/statistics/base/probit.m,
scripts/statistics/base/qqplot.m, scripts/statistics/distributions/betacdf.m,
scripts/statistics/distributions/betainv.m,
scripts/statistics/distributions/betapdf.m,
scripts/statistics/distributions/betarnd.m,
scripts/statistics/distributions/binocdf.m,
scripts/statistics/distributions/binoinv.m,
scripts/statistics/distributions/binopdf.m,
scripts/statistics/distributions/binornd.m,
scripts/statistics/distributions/cauchy_cdf.m,
scripts/statistics/distributions/cauchy_inv.m,
scripts/statistics/distributions/cauchy_pdf.m,
scripts/statistics/distributions/cauchy_rnd.m,
scripts/statistics/distributions/chi2cdf.m,
scripts/statistics/distributions/chi2inv.m,
scripts/statistics/distributions/chi2pdf.m,
scripts/statistics/distributions/chi2rnd.m,
scripts/statistics/distributions/discrete_cdf.m,
scripts/statistics/distributions/discrete_inv.m,
scripts/statistics/distributions/discrete_pdf.m,
scripts/statistics/distributions/discrete_rnd.m,
scripts/statistics/distributions/empirical_cdf.m,
scripts/statistics/distributions/empirical_inv.m,
scripts/statistics/distributions/empirical_pdf.m,
scripts/statistics/distributions/empirical_rnd.m,
scripts/statistics/distributions/expcdf.m,
scripts/statistics/distributions/expinv.m,
scripts/statistics/distributions/exppdf.m,
scripts/statistics/distributions/exprnd.m,
scripts/statistics/distributions/fcdf.m,
scripts/statistics/distributions/finv.m,
scripts/statistics/distributions/fpdf.m,
scripts/statistics/distributions/frnd.m,
scripts/statistics/distributions/gamcdf.m,
scripts/statistics/distributions/gaminv.m,
scripts/statistics/distributions/gampdf.m,
scripts/statistics/distributions/gamrnd.m,
scripts/statistics/distributions/geocdf.m,
scripts/statistics/distributions/geoinv.m,
scripts/statistics/distributions/geopdf.m,
scripts/statistics/distributions/geornd.m,
scripts/statistics/distributions/hygecdf.m,
scripts/statistics/distributions/hygeinv.m,
scripts/statistics/distributions/hygepdf.m,
scripts/statistics/distributions/hygernd.m,
scripts/statistics/distributions/kolmogorov_smirnov_cdf.m,
scripts/statistics/distributions/laplace_cdf.m,
scripts/statistics/distributions/laplace_inv.m,
scripts/statistics/distributions/laplace_pdf.m,
scripts/statistics/distributions/laplace_rnd.m,
scripts/statistics/distributions/logistic_cdf.m,
scripts/statistics/distributions/logistic_inv.m,
scripts/statistics/distributions/logistic_pdf.m,
scripts/statistics/distributions/logistic_rnd.m,
scripts/statistics/distributions/logncdf.m,
scripts/statistics/distributions/logninv.m,
scripts/statistics/distributions/lognpdf.m,
scripts/statistics/distributions/lognrnd.m,
scripts/statistics/distributions/nbincdf.m,
scripts/statistics/distributions/nbininv.m,
scripts/statistics/distributions/nbinpdf.m,
scripts/statistics/distributions/nbinrnd.m,
scripts/statistics/distributions/normcdf.m,
scripts/statistics/distributions/norminv.m,
scripts/statistics/distributions/normpdf.m,
scripts/statistics/distributions/normrnd.m,
scripts/statistics/distributions/poisscdf.m,
scripts/statistics/distributions/poissinv.m,
scripts/statistics/distributions/poisspdf.m,
scripts/statistics/distributions/poissrnd.m,
scripts/statistics/distributions/stdnormal_cdf.m,
scripts/statistics/distributions/stdnormal_inv.m,
scripts/statistics/distributions/stdnormal_pdf.m,
scripts/statistics/distributions/stdnormal_rnd.m,
scripts/statistics/distributions/tcdf.m,
scripts/statistics/distributions/tinv.m,
scripts/statistics/distributions/tpdf.m,
scripts/statistics/distributions/trnd.m,
scripts/statistics/distributions/unidcdf.m,
scripts/statistics/distributions/unidinv.m,
scripts/statistics/distributions/unidpdf.m,
scripts/statistics/distributions/unidrnd.m,
scripts/statistics/distributions/unifcdf.m,
scripts/statistics/distributions/unifinv.m,
scripts/statistics/distributions/unifpdf.m,
scripts/statistics/distributions/unifrnd.m,
scripts/statistics/distributions/wblcdf.m,
scripts/statistics/distributions/wblinv.m,
scripts/statistics/distributions/wblpdf.m,
scripts/statistics/distributions/wblrnd.m,
scripts/statistics/distributions/wienrnd.m,
scripts/statistics/models/logistic_regression.m,
scripts/statistics/models/private/logistic_regression_derivatives.m,
scripts/statistics/models/private/logistic_regression_likelihood.m,
scripts/statistics/tests/anova.m, scripts/statistics/tests/bartlett_test.m,
scripts/statistics/tests/chisquare_test_homogeneity.m,
scripts/statistics/tests/chisquare_test_independence.m,
scripts/statistics/tests/cor_test.m,
scripts/statistics/tests/f_test_regression.m,
scripts/statistics/tests/hotelling_test.m,
scripts/statistics/tests/hotelling_test_2.m,
scripts/statistics/tests/kolmogorov_smirnov_test.m,
scripts/statistics/tests/kolmogorov_smirnov_test_2.m,
scripts/statistics/tests/kruskal_wallis_test.m,
scripts/statistics/tests/manova.m, scripts/statistics/tests/mcnemar_test.m,
scripts/statistics/tests/prop_test_2.m, scripts/statistics/tests/run_test.m,
scripts/statistics/tests/sign_test.m, scripts/statistics/tests/t_test.m,
scripts/statistics/tests/t_test_2.m,
scripts/statistics/tests/t_test_regression.m,
scripts/statistics/tests/u_test.m, scripts/statistics/tests/var_test.m,
scripts/statistics/tests/welch_test.m,
scripts/statistics/tests/wilcoxon_test.m, scripts/statistics/tests/z_test.m,
scripts/statistics/tests/z_test_2.m:
Removed functions.
author | Rik <rik@octave.org> |
---|---|
date | Sun, 07 Jan 2018 17:13:23 -0800 |
parents | scripts/statistics/base/mode.m@194eb4bd202b |
children | 6652d3823428 |
line wrap: on
line source
## Copyright (C) 2007-2017 David Bateman ## ## 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 ## <https://www.gnu.org/licenses/>. ## -*- texinfo -*- ## @deftypefn {} {} mode (@var{x}) ## @deftypefnx {} {} mode (@var{x}, @var{dim}) ## @deftypefnx {} {[@var{m}, @var{f}, @var{c}] =} mode (@dots{}) ## Compute the most frequently occurring value in a dataset (mode). ## ## @code{mode} determines the frequency of values along the first non-singleton ## dimension and returns the value with the highest frequency. If two, or ## more, values have the same frequency @code{mode} returns the smallest. ## ## If the optional argument @var{dim} is given, operate along this dimension. ## ## The return variable @var{f} is the number of occurrences of the mode in ## the dataset. ## ## The cell array @var{c} contains all of the elements with the maximum ## frequency. ## @seealso{mean, median} ## @end deftypefn function [m, f, c] = mode (x, dim) if (nargin < 1 || nargin > 2) print_usage (); endif if (! (isnumeric (x) || islogical (x))) error ("mode: X must be a numeric vector or matrix"); endif nd = ndims (x); sz = size (x); if (nargin < 2) ## Find the first non-singleton dimension. (dim = find (sz > 1, 1)) || (dim = 1); else if (! (isscalar (dim) && dim == fix (dim) && dim > 0)) error ("mode: DIM must be an integer and a valid dimension"); endif endif if (dim > nd) ## Special case of mode over non-existent dimension. m = x; f = ones (size (x)); c = num2cell (x); return; endif sz2 = sz; sz2(dim) = 1; sz3 = ones (1, nd); sz3(dim) = sz(dim); if (issparse (x)) t2 = sparse (sz(1), sz(2)); else t2 = zeros (sz); endif if (dim != 1) perm = [dim, 1:dim-1, dim+1:nd]; t2 = permute (t2, perm); endif xs = sort (x, dim); t = cat (dim, true (sz2), diff (xs, 1, dim) != 0); if (dim != 1) t2(permute (t != 0, perm)) = diff ([find(permute (t, perm))(:); prod(sz)+1]); f = max (ipermute (t2, perm), [], dim); xs = permute (xs, perm); else t2(t) = diff ([find(t)(:); prod(sz)+1]); f = max (t2, [], dim); endif c = cell (sz2); if (issparse (x)) m = sparse (sz2(1), sz2(2)); else m = zeros (sz2, class (x)); endif for i = 1 : prod (sz2) c{i} = xs(t2(:, i) == f(i), i); m(i) = c{i}(1); endfor endfunction %!test %! [m, f, c] = mode (toeplitz (1:5)); %! assert (m, [1,2,2,2,1]); %! assert (f, [1,2,2,2,1]); %! assert (c, {[1;2;3;4;5],[2],[2;3],[2],[1;2;3;4;5]}); %!test %! [m, f, c] = mode (toeplitz (1:5), 2); %! assert (m, [1;2;2;2;1]); %! assert (f, [1;2;2;2;1]); %! assert (c, {[1;2;3;4;5];[2];[2;3];[2];[1;2;3;4;5]}); %!test %! a = sprandn (32, 32, 0.05); %! sp0 = sparse (0); %! [m, f, c] = mode (a); %! [m2, f2, c2] = mode (full (a)); %! assert (m, sparse (m2)); %! assert (f, sparse (f2)); %! c_exp(1:length (a)) = { sp0 }; %! assert (c ,c_exp); %! assert (c2,c_exp); %!assert (mode ([2,3,1,2,3,4],1),[2,3,1,2,3,4]) %!assert (mode ([2,3,1,2,3,4],2),2) %!assert (mode ([2,3,1,2,3,4]),2) %!assert (mode (single ([2,3,1,2,3,4])), single (2)) %!assert (mode (int8 ([2,3,1,2,3,4])), int8 (2)) %!assert (mode ([2;3;1;2;3;4],1),2) %!assert (mode ([2;3;1;2;3;4],2),[2;3;1;2;3;4]) %!assert (mode ([2;3;1;2;3;4]),2) %!test %! x = magic (3); %! [m, f, c] = mode (x, 3); %! assert (m, x); %! assert (f, ones (3,3)); %! assert (c, num2cell (x)); %!shared x %! x(:,:,1) = toeplitz (1:3); %! x(:,:,2) = circshift (toeplitz (1:3), 1); %! x(:,:,3) = circshift (toeplitz (1:3), 2); %!test %! [m, f, c] = mode (x, 1); %! assert (reshape (m, [3, 3]), [1 1 1; 2 2 2; 1 1 1]); %! assert (reshape (f, [3, 3]), [1 1 1; 2 2 2; 1 1 1]); %! c = reshape (c, [3, 3]); %! assert (c{1}, [1; 2; 3]); %! assert (c{2}, 2); %! assert (c{3}, [1; 2; 3]); %!test %! [m, f, c] = mode (x, 2); %! assert (reshape (m, [3, 3]), [1 1 2; 2 1 1; 1 2 1]); %! assert (reshape (f, [3, 3]), [1 1 2; 2 1 1; 1 2 1]); %! c = reshape (c, [3, 3]); %! assert (c{1}, [1; 2; 3]); %! assert (c{2}, 2); %! assert (c{3}, [1; 2; 3]); %!test %! [m, f, c] = mode (x, 3); %! assert (reshape (m, [3, 3]), [1 2 1; 1 2 1; 1 2 1]); %! assert (reshape (f, [3, 3]), [1 2 1; 1 2 1; 1 2 1]); %! c = reshape (c, [3, 3]); %! assert (c{1}, [1; 2; 3]); %! assert (c{2}, [1; 2; 3]); %! assert (c{3}, [1; 2; 3]); ## Test input validation %!error mode () %!error mode (1, 2, 3) %!error <X must be a numeric> mode ({1 2 3}) %!error <DIM must be an integer> mode (1, ones (2,2)) %!error <DIM must be an integer> mode (1, 1.5) %!error <DIM must be .* a valid dimension> mode (1, 0)