comparison toolboxes/FullBNT-1.0.7/netlab3.3/gtmlmode.m @ 0:e9a9cd732c1e tip

first hg version after svn

+function modes = gtmlmode(net, data)
+%GTMLMODE Mode responsibility for data in a GTM.
+%
+%	Description
+%	 MODES = GTMLMODE(NET, DATA) takes a GTM structure NET, and computes
+%	the modes of the responsibility  distributions for each data point in
+%	DATA.  These will always lie at one of the latent space sample points
+%	NET.X.
+%
+%	See also
+%	GTM, GTMPOST, GTMLMEAN
+%
+
+%	Copyright (c) Ian T Nabney (1996-2001)
+
+% Check for consistency
+errstring = consist(net, 'gtm', data);
+if ~isempty(errstring)
+  error(errstring);
+end
+
+R = gtmpost(net, data);
+% Mode is maximum responsibility
+[max_resp, max_index] = max(R, [], 2);
+modes = net.X(max_index, :);