/usr/share/octave/packages/tsa-4.3.3/histo2.m is in octave-tsa 4.3.3-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 | function R = histo2(Y, W)
% HISTO2 calculates histogram for multiple columns with separate bin values
% for each data column.
%
% R = HISTO2(Y)
% R = HISTO2(Y, W)
% Y data
% W weight vector containing weights of each sample,
% number of rows of Y and W must match.
% default W=[] indicates that each sample is weighted with 1.
%
% R = HISTO(...)
% R is a struct with th fields
% R.X the bin-values, bin-values are computed separately for each
% data column, thus R.X is a matrix, each column contains the
% the bin values of for each data column, unused elements are indicated with NaN.
% In order to have common bin values, use HISTO3.
% R.H is the frequency of occurence of value X
% R.N are the number of valid (not NaN) samples (i.e. sum of weights)
%
% more histogram-based results can be obtained by HIST2RES2
%
% see also: HISTO, HISTO2, HISTO3, HISTO4
%
% REFERENCE(S):
% C.E. Shannon and W. Weaver "The mathematical theory of communication" University of Illinois Press, Urbana 1949 (reprint 1963).
% $Id: histo2.m 12710 2014-09-23 00:39:17Z schloegl $
% Copyright (C) 1996-2002,2008,2011,2014 by Alois Schloegl <alois.schloegl@ist.ac.at>
% This is part of the TSA-toolbox
% http://pub.ist.ac.at/~schloegl/matlab/tsa/
%
% This program 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.
%
% This program 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 this program. If not, see <http://www.gnu.org/licenses/>.
%%%%% check input arguments %%%%%
[yr,yc] = size(Y);
if nargin < 2,
W = [];
end;
if ~isempty(W) && (yr ~= numel(W)),
error('number of rows of Y does not match number of elements in W');
end;
%%%%% identify all possible X's and generate overall Histogram %%%%%
N = sum(~isnan(Y), 1);
if all(N==0)
H=[];
X=[];
end;
NN = N;
if isempty(W)
sY = sort(Y,1);
else
[sY, idx] = sort(Y,1);
W = cumsum(W(idx)); %% W becomes cumulative sum
end;
[ix,iy] = find( diff(sY, [], 1) > 0);
nn0 = 0;
for k = 1:yc,
tmp = [ix(iy==k); N(k)];
nn1 = length(tmp);
if isempty(W)
H(1:nn1,k) = [tmp(1); diff(tmp)];
else
%%% Note that W is the cumulative sum
H(1:nn1,k) = [W(tmp(1),k); diff(W(tmp,k))];
NN(k) = W(N(k), k);
end;
X(1:nn1, k) = sY(tmp, k);
if k==1;
nn0 = nn1;
elseif nn1 < nn0,
H (1+nn1:nn0, k) = NaN;
X (1+nn1:nn0, k) = NaN;
elseif nn1 > nn0,
H (1+nn0:nn1, 1:k-1) = NaN;
X (1+nn0:nn1, 1:k-1) = NaN;
nn0 = nn1;
end;
end;
R.datatype = 'HISTOGRAM';
R.H = H;
R.X = X;
R.N = NN;
%!assert(getfield(histo2([]),'N'), 0)
%!assert(getfield(histo2(1),'N'), 1)
%!assert(getfield(histo2([1;1]),'H'), 2)
|