File Exchange

image thumbnail

N-Dimensional Histogram Count

version 1.5 (2.46 KB) by

Bins the elements of n matrices into an n-dimensional histogram.

0 Downloads

Updated

View License

HISTCND is similar to HISTC, but creates a histogram with any number of dimensions.

N = HISTCND(X,Y,Z,... XEDGES,YEDGES,ZEDGES,... )

where N is a histogram count with dimensions
         length(XEDGES) x length(YEDGES) x length(ZEDGES) ...

If XEDGES, YEDGES, etc. are monotonically increasing and non-NaN, a data point is assigned to bin N(i,j,k,...) if
  
   XEDGES(i) <= X < XEDGES(i+1)
 
   YEDGES(j) <= Y < YEDGES(j+1)
  
   ZEDGES(k) <= Z < ZEDGES(k+1)

   ....

Note: data outside the ranges of the EDGES vectors are excluded from the histogram, and not placed in the first or last bins.

An example:

x=rand(100,100);
y=rand(100,100);
xedges=0:0.01:1;
yedges=0:0.1:1;

n=histcnd(x,y,xedges,yedges);

pcolor(xedges,yedges,n');
colorbar;

Comments and Ratings (2)

Mathew

Mathew (view profile)

Hi, I've been unable to produce this error. I tried with 4 columns of 1,000,000 data points, with 20 edge points for each dimension.

Using too many edge points (histogram bins) will make the output array too big for matlab to handle.

In the linux program 'top', I saw the MATLAB process was using around 4% of the memory (system has 4GB). Maybe you could try the code below, and see if it works on yours.
Cheers,
Mat

>> x1=rand(1000000,1);
>> x2=rand(1000000,1);
>> x3=rand(1000000,1);
>> x4=rand(1000000,1);
>> e1=linspace(0,1,20);
>> e2=e1;e3=e1;e4=e1;
>> tic;n=histcnd(x1,x2,x3,x4,e1,e2,e3,e4);toc;

Elapsed time is 2.554106 seconds.

>> size(n)

ans =

    20 20 20 20

Justin

Justin (view profile)

Hello,

I'm trying to use your program but it gives me an error that says "out of memory". I've been feeding in 4 columns of 10000 data points into the program.

Regards.

Updates

1.5

Bug Fix: the behaviour is now consistent for data above and below ranges of the edges vectors.

1.3

Added comments to explain the algorithm. Small optimization of code.

1.2

Added an example to the description.

MATLAB Release
MATLAB 7.3 (R2006b)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video