/ani/mrses

To get this branch, use:
bzr branch http://suren.me/webbzr/ani/mrses

« back to all changes in this revision

Viewing changes to mrses_mtx.m

  • Committer: Suren A. Chilingaryan
  • Date: 2010-04-28 04:30:08 UTC
  • Revision ID: csa@dside.dyndns.org-20100428043008-vd9z0nso9axezvlp
Initial import

Show diffs side-by-side

added added

removed removed

Lines of Context:
mrses_mtx.m
 
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;