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
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2003 - INRIA - Carlos Klimann
//
// 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 s = samwr(sizam, numsamp, X)
//
//This function gives a s matrix sizsam x numsamp. It
//contains numsamp random samples (the columns) each of
//sizam ( < size(X, '*')) extractions, without replacement,
//from the vector X.
//
if argn(2) <> 3 then
error(msprintf(gettext("%s: Wrong number of input argument: %d expected.\n"), "samwr", 3)),
end
if X == []|sizam == 0|numsamp == 0 then s = []; return; end
sizx = size(X, "*")
if sizam > sizx then
error(msprintf(gettext("%s: Wrong value of input argument #%d: Lower than or equal to size of input argument #%d expected.\n"), "samwr", 1, 3)),
end
perm = grand(numsamp, "prm", (1:sizx)')
s = matrix(X(perm), -1, numsamp)
s = s(1:sizam, :)
endfunction
|
