bzr branch
http://darksoft.org/webbzr/ani/mrses
1
by Suren A. Chilingaryan
Initial import |
1 |
function res=mrses_mtx(A,B,k,Niter,Ncycle,distmod,block) |
2 |
if (nargin<7) |
|
3 |
block=256; |
|
4 |
end
|
|
5 |
if (nargin<6) |
|
6 |
distmod=1; |
|
7 |
end
|
|
8 |
if (nargin<5) |
|
9 |
Ncycle=1000; |
|
10 |
end
|
|
11 |
if (nargin<4) |
|
12 |
Niter=500; |
|
13 |
end
|
|
14 |
if (nargin<3) |
|
15 |
k=5; |
|
16 |
end
|
|
17 |
if (nargin<2) |
|
18 |
error('As minimum two matrixes needed for MRSES'); |
|
19 |
end
|
|
20 |
if (nargin>6) |
|
21 |
error('Too much parameters'); |
|
22 |
end
|
|
23 |
||
24 |
sa=size(A);sb=size(B); |
|
25 |
if (sa(2)==sb(2)) |
|
26 |
genes=sa(2); |
|
27 |
else
|
|
28 |
error('Features dimension mismatch'); |
|
29 |
end
|
|
30 |
||
31 |
ctx = mrses_hw(); |
|
32 |
mrses_hw(ctx, 1, k, block, single(A), single(B), distmod); |
|
33 |
optki=[]; |
|
34 |
for icycle=0:block:(Ncycle-1) |
|
35 |
block_size = min(block, Ncycle - icycle); |
|
36 |
% optki=mrses_hw(ctx, 15, Niter, block_size);
|
|
37 |
optki = [optki, mrses_hw(ctx, 15, Niter, block_size)]; |
|
38 |
end
|
|
39 |
% n = mrses_hw(ctx, 16); |
|
40 |
mrses_hw(ctx); |
|
41 |
||
42 |
%H=[double(n)./Ncycle;1:genes];
|
|
43 |
%res=flipud(sortrows(H'));
|
|
44 |
res=optki; |