Kelly Kearney

I've been having some difficulty with this function when trying to animate figures that include a large mix of object types, particularly ones with the same "type" value but different properties. For example, lines produced by line vs plot, hggroup objects of different types, etc. Anymate fails when trying to gather properties of these due to its reliance on a single sample handle from each type. An easy fix (or at least, easier than trying to handle every property gracefully) would be to have a switch to exclude certain properties from being collected (right now, excluding from anymation doesn't exclude from collection).

Overall, this is a great function, and I use it quite often to analyze model ensemble output. A few enhancements that could be nice:

- Add the option to display in a horizontal orientation.

- Add the option to filter outliers when calculating bin widths and kernal densities. Could also be nice to display these as points, as in boxplot, rather than connecting them via long lines to the main histogram.

- This is an edge case, but the function will error under the addSpread option if a column/group contains only NaNs and/or Infs.

@Mukhtar Ullah,

I don't actually manually set the interpreter in this function, so it relies on the default. You can actually change this after the fact using set, as long as you return the object_h handles (see last syntax option); the first n indices of this array will correspond to the n text objects in the legend. This slightly obscure way of returning handles was chosen to be consistent with the original legend command.

great function! This saved me many hours of working around all the bugs in contourf.

The only (minor) short coming is that this function doesn't deal with saturated low values. Eg. if you enter:
contourfcmap(x,y,z,clev,cmap,lo,hi)
and there are values in z < min(clev), they will be plotted as white spaces. This is an issue with contourf really, which I won't get into here. An easy fix is set all values in z outside the clev limits to min(clev) or max(clev). Eg.
toohigh = z>max(clev); z(toohigh) = max(clev);
loolow = z<min(clev); z(toolow) = min(clev);