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
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 1999 - INRIA - Carlos Klimann
// Copyright (C) 2003 - Jean-Sebastien Giet & Bruno Pincon
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
//
function [m] = tabul(X, order)
//
// PURPOSE
// This function computes the frequency of values of
// the components of a vector or matrix X of numbers or
// string characters.
//
// If X is a numerical vector or matrix then
// m is a two column matrix who contains in
// the first column the distinct values of X
// and in the other column the number of occurrences
// of those values (m(i,2) is the number of occurrences
// of m(i,1)).
//
// If X is a vector or matrix of strings, m is a list
// whose first member is a string vector composed with
// the distinct values of X and the second member is a
// vector whose components are the number of occurrences
// of those values ( m(i)(2) is the number of occurrences
// of the string m(i)(1) ).
//
// The optional parameter order must be "d" or "i"
// (by default order = "d") and it gives the order
// of the distinct vector values of X (first column
// or first part of m) :
// order = "d" means that these values are sorted
// in decreasing order
// = "i" means by increasing order
//
// AUTHORS
// Original version by Carlos Klimann
// This version by Jean-Sebastien Giet & Bruno Pincon
//
// date: 1999-04-09 (original version)
// 2003-Mars-26 (new version)
//
// NOTES
// The new version :
// is faster (by using no more loop)
// performs a complete check of the input arguments
// add the order option
// used gsort in place of sort : the sort function
// sorts strings vector not the usual way !
rhs = argn(2)
if rhs<1 | 2<rhs then
error(msprintf(gettext("%s: Wrong number of input arguments: %d to %d expected.\n"),"tabul",1,2)),
elseif rhs == 1 then
order = "d"
end
typeX = type(X)
if typeX ~= 1 & typeX ~= 10 then
error(msprintf(gettext("%s: Wrong type for input argument #%d: Vector, matrix of numbers or strings expected.\n"),"tabul",1))
end
if type(order) ~= 10 then
error(msprintf(gettext("%s: Wrong value for input argument #%d: ''%s'' or ''%s'' expected.\n"),"tabul",2,"i","d"))
end
if order~="i" & order~="d" then
error(msprintf(gettext("%s: Wrong value for input argument #%d: ''%s'' or ''%s'' expected.\n"),"tabul",2,"i","d"))
end
if ( X == [] ) then
m = %nan
return
end
X = X(:)
X = gsort(X,"g",order)
n = size(X,"*")
ind = [find(X(1:$-1)~=X(2:$)) n]
val = X(ind)
occ = diff([0 ind])'
if typeX == 1 then
m = [val occ]
else // X (and so val) is a vector of strings
m = list(val,occ)
end
endfunction
|
