view toolboxes/FullBNT-1.0.7/KPMstats/test_dir.m @ 0:cc4b1211e677 tip

initial commit to HG from
Changeset:
646 (e263d8a21543) added further path and more save "camirversion.m"

% # of sample points
n_samples = 1000;

p = ones(3,1)/3;

% Low Entropy
alpha = 0.5*p;

% High Entropy
%alpha = 10*p;

% draw n_samples random points from the 3-d dirichlet(alpha),
% and plot the results
points = zeros(3,n_samples);
for i = 1:n_samples
    points(:,i) = dirichletrnd(alpha);
end

scatter3(points(1,:)', points(2,:)', points(3,:)', 'r', '.', 'filled');