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
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
|
function [resu, idx] = sort(df, varargin)
%# -*- texinfo -*-
%# @deftypefn {Loadable Function} {[@var{s}, @var{i}] =} sort (@var{x})
%# @deftypefnx {Loadable Function} {[@var{s}, @var{i}] =} sort (@var{x}, @var{dim})
%# @deftypefnx {Loadable Function} {[@var{s}, @var{i}] =} sort (@var{x}, @var{mode})
%# @deftypefnx {Loadable Function} {[@var{s}, @var{i}] =} sort (@var{x}, @var{dim}, @var{mode})
%# Return a copy of @var{x} with the elements arranged in increasing
%# order. For matrices, @code{sort} orders the elements in each column.
%#
%# For example:
%#
%# @example
%# @group
%# sort ([1, 2; 2, 3; 3, 1])
%# @result{} 1 1
%# 2 2
%# 3 3
%#
%# @end group
%# @end example
%#
%# The @code{sort} function may also be used to produce a matrix
%# containing the original row indices of the elements in the sorted
%# matrix. For example:
%#
%# @example
%# @group
%# [s, i] = sort ([1, 2; 2, 3; 3, 1])
%# @result{} s = 1 1
%# 2 2
%# 3 3
%# @result{} i = 1 3
%# 2 1
%# 3 2
%# @end group
%# @end example
%#
%# If the optional argument @var{dim} is given, then the matrix is sorted
%# along the dimension defined by @var{dim}. The optional argument @code{mode}
%# defines the order in which the values will be sorted. Valid values of
%# @code{mode} are `ascend' or `descend'.
%#
%# For equal elements, the indices are such that the equal elements are listed
%# in the order that appeared in the original list.
%#
%# The @code{sort} function may also be used to sort strings and cell arrays
%# of strings, in which case the dictionary order of the strings is used.
%#
%# The algorithm used in @code{sort} is optimized for the sorting of partially
%# ordered lists.
%# @end deftypefn
%% Copyright (C) 2009-2017 Pascal Dupuis <cdemills@gmail.com>
%%
%% This file is part of the dataframe package for Octave.
%%
%% This package 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, or (at your option) any later version.
%%
%% This package 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 package; see the file COPYING. If not,
%% see <http://www.gnu.org/licenses/>.
%#
%# $Id$
%#
if (~isa (df, 'dataframe'))
resu = []; return;
end
dim = []; mode = [];
vout= varargin;
indi = 1; while (indi <= length (varargin))
if (isnumeric (varargin{indi}))
if (~isempty (dim))
print_usage ('@dataframe/sort');
resu = [];
return
else
dim = varargin{indi};
if (3 == dim) vout(indi) = 2; end
end
else
if (~isempty (mode))
print_usage ('@dataframe/sort');
resu = [];
return
else
sort = varargin{indi};
end
end
indi = indi + 1;
end;
if (isempty (dim)) dim = 1; end
%# pre-assignation
resu = struct (df);
switch(dim)
case {1}
for indi = (1:resu.x_cnt(2))
[resu.x_data{indi}, idx(:, indi, :)] = sort ...
(resu.x_data{indi}(:, resu.x_rep{indi}), varargin{:});
resu.x_data{indi} = squeeze (resu.x_data{indi});
resu.x_rep{indi} = 1:size(resu.x_data{indi}, 2);
end
if (all ([1 == size(idx, 2) 1 == size(idx, 3)]))
if (size (resu.x_ridx, 1) == resu.x_cnt(1))
resu.x_ridx = resu.x_ridx(idx, :);
end
if (~isempty (resu.x_name{1, 1}))
resu.x_name{1, 1} = resu.x_name{1, 1}(idx);
resu.x_over{1, 1} = resu.x_over{1, 1}(idx);
end
else
%# data where mixed
resu.x_ridx = idx;
resu.x_name{1, 1} = []; resu.x_over{1, 1} = [];
end
case {2}
error ('Operation not implemented');
case {3}
for indi = (1:resu.x_cnt(2))
[resu.x_data{1, indi}, idx(:, indi)] = sort (resu.x_data{1, indi}, vout(:));
end
otherwise
error ('Invalid dimension %d', dim);
end
dummy = dbstack ();
if (any (strmatch ('quantile', {dummy.name})))
resu = df_whole (resu);
else
resu = dataframe (resu);
end
end
|
