Category plots, Numeric Values vs Numeric Bins

Okay, here's the short version, and for those interested the longer version:

Short version:
In a category plot of numeric values vs numeric values, how do I center the bars over their x-values?
The top of the bar starts at the value on the bottom axis and is drawn to the next value, BUT the final point's horizontal top is not drawn!
What am I missing?

A plot is attached, and here's some code to plot 3 elucidating examples:
//Graph1: Numeric vs Numeric; Markers
Make /O Vals = {2,3,0,0,1,2,0,2,1}, Bins = {5,10,15,20,25,30,35,40,45};
Display Vals vs Bins; Modifygraph mode=3, msize=4, marker=16, rgb=(65535,0,0);
Label bottom "Bin Center"; Label left "# in each bin";  SetAxis bottom 0, 50; SetAxis left 0, 4
//Graph2: Numeric vs Numeric; Markers and Bars
Display Vals vs Bins; Modifygraph mode=3, msize=4, marker=16, rgb=(65535,0,0);
Label bottom "Bin Center"; Label left "# in each bin";  SetAxis bottom 0, 50; SetAxis left 0, 4
AppendToGraph Vals vs Bins; ModifyGraph mode[1]=5, rgb[1]=(0,0,0);
//Graph3: Text Category
Make /T/N=(dimsize(Bins,0)) Bins_text = Num2Str(Bins[p])
Display Vals vs Bins_text;
Label bottom "Bin Center"; Label left "# in each bin"; SetAxis bottom -1, 10; SetAxis left 0, 4
Execute "TileWindows "+WinName(2,1)+" "+WinName(1,1)+" "+WinName(0,1)



The Longer Version:

I have sets of 3D waves with dimensions (16,8,100). I want to find the histogram of each "beam" in a given wave, store it for later, and then plot it. For example, I want to have the histogram of all points at (0,0,r). Then repeat for all beams (p,q).

Now, histogram won't do multidimensions like this, I think. So I decided to duplicate the 3D input wave, pull out a beam, find the 1D histogram, and then shove the histogram back into the duplicated 3D wave at the appropriate (p,q) position. In this manner the histograms fill up the beams of the new wave, at the appropriate (p,q) coordinates.

However, histogram stores the bin information in the output wave's scaling. So I decided to add a chunk to the output wave, and store the bin values up a beam as well.

So Display OutputWave[0][0][][0] vs OutputWave[0][0][][1]; ModifyGraph mode=5 should spit back out the histogram for the (0,0) pixel.
The coding works okay I think, but I got snagged on the graphing as bars! I can do it as sticks and markers, but kinda like the bars better.
So that's why I'm wondering how to get the "last point" as a bar, or alternatively center the bars on the actual x-values.
I attach the code I wrote to do this, feel free to comment/improve/correct.

Thanks in advance!
HistogramLayers.ipf
How about using an x-offset for the graph?
make /N=11 /O test
setscale /P x,10,2, test
test=exp(-1/10*(x-20)^2)
display test
ModifyGraph manTick(bottom)={0,2,0,0},manMinor(bottom)={0,0}
ModifyGraph offset={-0.5*dimdelta(test,0),0}  // This one is doing the tweak: shift the graph left by half an increment

Maybe the axis range needs to be fixed to a slightly larger value then actual data range as well.
HJ
Concerning the last data point, it seems that an additional dummy 'bin value' is necessary since the width of the bars is defined by the spacing in the bin wave. Try to add another data point in the 'Bin' wave and play with its value (I recommend to try 46, 15, and 50.). The data wave should be NaN at that point.
Looks more like a feature than like a bug to me...
HJ
HJDrescher wrote:
Concerning the last data point, it seems that an additional dummy 'bin value' is necessary since the width of the bars is defined by the spacing in the bin wave. Try to add another data point in the 'Bin' wave and play with its value (I recommend to try 46, 15, and 50.). The data wave should be NaN at that point.


Thanks HJDrescher,

It seems bar plots behave differently when you plot a single wave vs its scaling, and when you plot one wave vs a separate x-wave.
For plotting a single wave vs its own scaling, yes I think your suggestion of an offset would work famously.

However, for a wave vs a separate x-wave IGOR still doesn't draw the last bar, even after changing offset & axis settings.

Attached is a multiplot window created with the commands below.
Top plot is wave test vs its internal scaling
Middle plot is wave test_y vs wave test_x
Bottom plot is wave test_y_addpt vs test_x_addpt where I've artificially added a point with value NaN, 32 to y and x, respectively.

The last method (thanks for the suggestion!) could work but for my case it seems like maybe more trouble than it's worth. I have 128 histograms stored in a single multi-dimensional wave, would like to plot them all on a graph, and then use ModifyGraph hidetrace to show the ones I'm interested in. Figuring out the correct scaling for each histogram's additional point and appending it to the multidimensional wave could (somewhat easily...) be done but I'm hesitant to add "fake" data to an actual data set from our lab... Oh well.

Quote:
Looks more like a feature than like a bug to me...
HJ


Yes clearly a "feature"!
Cheers,
Dave


Display/N=BarHost
Make /N=11/O test; SetScale /P x, 10, 2, test;
test = exp(1/100*(x-20)^2)
Display/HOST=BarHost/N=VsScale/W=(0,0,1,0.333) test, test
ModifyGraph mode[0]=5, mode[1]=3, marker=16, msize=4, rgb[1]=(0,0,65535)
SetAxis bottom 5, 35; SetAxis left 0, 3; Label left "# in each bin"; Label bottom "Bin Center"
TextBox /W=BarHost#VsScale/C/N=Title/A=MT/X=0/Y=1/G=(0,0,65535)/F=0/B=1 "\JC\Zr150y-wave vs scaling"
Duplicate /O test test_x, test_y;
test_x = DimOffset(test,0) + DimDelta(test,0)*p
Display/HOST=BarHost/N=VsXwave/W=(0,0.333,1,0.666) test_y vs test_x;
AppendToGraph /C=(0,0,65535) test_y vs test_x
ModifyGraph mode[0]=5, mode[1]=3, marker=16, msize=4
SetAxis bottom 5, 35; SetAxis left 0, 3; Label left "# in each bin"; Label bottom "Bin Center"
TextBox /W=BarHost#VsXwave/C/N=Title/A=MT/X=0/Y=1/G=(0,0,65535)/F=0/B=1 "\JC\Zr150y-wave vs x-wave"
Duplicate /O test_x test_x_addpt; Duplicate /O test_y test_y_addpt;
InsertPoints 11, 1, test_x_addpt, test_y_addpt
test_y_addpt[11] = NaN; test_x_addpt[11]=32
Display/HOST=BarHost/N=VsXwaveAddpt/W=(0,0.666,1,1) test_y_addpt test_y_addpt vs test_x_addpt
ModifyGraph mode[0]=5, mode[1]=3, marker=16, msize=4, rgb[1]=(0,0,65535)
SetAxis bottom 5, 35; SetAxis left 0, 3; Label left "# in each bin"; Label bottom "Bin Center"
TextBox /W=BarHost#VsXwaveAddpt/C/N=Title/A=MT/X=0/Y=1/G=(0,0,65535)/F=0/B=1 "\JC\Zr150y-wave vs x-wave \radded point"
BarHost.png
Actually, Igor does not behave differently. It plots the bar to the *calculated* position of the next data point. For scaled waves this is obvious. How about two N-point waves: The data part is simple. The position however has a twist. Igor looks up position N+1 in the wave (which does not exist). What happens in this case? I think, Igor produces and catches a run time error (out of range) and returns the last valid entry of the wave (somewhere said in the manual - or try make /N=100 tst=p;print tst2[200]. Hence, it would plot the last bar, e.g., from 42 to 42 but does not for some reason (to avoid a spike probably - actually, this might be a bug...)

The 'feature' part: it provides non-equal widths of the bars - which might be handy at some point.
If I remember correctly, the additional data point is also necessary for images of the 'y vs x' kind (also defining the 'height of the related pixel'). This is quite useful if you want to plot color coded spectra in a temperature series with different steps according the interesting reaction temperatures (similar to http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3920861/figure/fig2/).

I understand your hesitation with the extra data value. However, NaN is somewhat special and might be used here. Good documentation is essential in such a case.

Another approach would leave the original data untouched and use a fixed (N+1) point display wave with the NaN as last data point assignment of only the actual N data points. Something like displaywave[0,N-1]=datawave[p][42] & scaling (Note: Zero based indexing). It might also run faster than hiding/showing a lot of traces (I use this method for 5000 waves with 2000 data points each; the hiding approach froze Igor for a couple of seconds - this method runs smoothly.)

Cheers,
HJ
HJ,
Thanks for clarifying the feature part of it, non-equal bar widths is a useful trick to have in my back pocket for later.
I've settled on going with sticks and markers for now, for ease of implementation.

Also your thoughts about changing a displaywave's definition instead of using hidetrace make sense and sound useful when dealing with lots of waves & points. In my case it's typically <50 waves with <30 points each, and so hidetrace also provides an easy way to keep all the data together when making graph copies as new experiments.

As a side note, I'm not sure why you say Igor does not behave differently. As you say it calculates the N+1 position from scaled data, presumably using something like DimDelta. It does not do this with two N-point waves, understandably, since there's no guarantee the N+1 point will be at the same interval. This is different behavior, no? In any case, this is probably worth it for the capability of unequal bar widths, as you mention.

Thanks again for your help!
Dave