histc_ic.m.svn-base

来自「bayesian network structrue learning mat」· SVN-BASE 代码 · 共 47 行

function [n,xechan] = histc_ic(x,edges)


%HISTC_IC  Histogram count 
%
%   [N XECHAN] = HISTC(X,EDGES), for vector X, counts the number of values
%   in X that fall between the elements in the EDGES vector 
%   (EDGES = cell array returned by HIST_IC function)
%
%   N is cell array containing these counts.
%    (or a vector if X is a column vector)
%   
%   XECHAN = discretized version of X
%
%
%   Example :
% 	X=randn(100,5);
%	Xapp=X(1:50,:);
%	Xtest=X(51:100,:);                        
%
%       % bins are computed with Xapp data
%	[n1,bornes,nbbornes,xechan]=hist_ic(Xapp);
%       % histogram is computed on Xtest data, with previously obtained bins
%	[n2,xtechan]=histc_ic(Xtest,bornes);
%
%
%   05-06-2001 by Ph. Leray
%

[nb_l,nb_c]=size(x);

% Outputs declaration
xechan=zeros(nb_l,nb_c);

if (nb_c==1)
        edges=mat2cell(edges);
end

for j=1:nb_c,
        [n{j} xechan(:,j)]=histc(x(:,j),[-inf edges{j} inf]);
        n{j}=n{j}(1:end-1);
end

if (nb_c==1)
        n=n{1};
end