**tabhist
in=**infile [parameter=value]

The median is also optionally computed, for which an alternate sorting routine can be selected (if available).

If multiple columns are used, they are appended to each other and treated
like a single long column. If you need simple statistics for multiple columns,
use **tabstat** instead.

Using the **nsigma** keyword, outliers can be removed
and statistics re-computed.

**in=***infile*- input file, in tabular format [no default].
**xcol=***column*- column(s) from which the (X) values are
taken. Multiple columns can be given in the usual
*nemoinp(3NEMO)*notation for integers. [default**1**]. **xmin=***x_min*- Minimum value along X-axis to include [default: autoscaling on minimum of datapoints].
**xmax=***x_max*- Maximum value along X-axis to include [default: autoscaling on maximum of datapoints].
**bins=***n_bins*- If one number is given, it is the number of (equal size) bins
between min and max. If more than one number is given, they are the bin
edges (thus one more than the number of bins need to be supplied). [default:
**16**]. **maxcount=***count*- Maximum count value per bin plotted along the Y-axis [default: autoscaling].
**nmax=***max_lines*- Maximum number of lines allowed to
read from datafile, if the data comes from a pipe as supposed from file.
A regular file will use the number of lines in the file as default. [Default:
**0**]. **ylog=t|f**- Take log of Y-axis? This option is really redundant, since
*tabmath(1NEMO)*can handle any transformation [default:**f**]. **xlab=***x-label*- Label along the
X-axis [default:
**value**]. **ylab=***y-label*- Label along the Y-axis [default:
**N**or**log(N)**]. **headline=***text*- Random verbiage, will be plotted along right top
of plot for identification [default: none]. The left top will contain the
filename
*infile*. **tab=t|f**- Table output? If false, a plot using your current
*yapp(5NEMO)*device will be given. If true, a simple ascii version of a histogram is shown. The default choice of**bins=16**keeps all the information on a simple 80*24 screen. [Default:**f**] **gauss=t|f**- If true, plot output will contain a
gaussian fit. The "best" gaussian fit is the one that has the same mean,
dispersion and integrated area as the data. [Default:
**t**]. **residual=t|f**- If
true, the residual (data - fit) will be plotted as a dashed line. [Default:
**t**]. **cumul=t|f**- Plot a cumulative histogram instead. If set, the keywords gauss,
residual, and ylog are automatically re-set to false. [Default:
**f**]. **median=t|f**- Compute median too? This can be timeconsuming for large numbers of points.
[Default:
**t**]. **torben=t|f**- Fast median calculator in case N is large. No sort
will be used, but instead an iterative O(N) method will be used. [Default:
**f**].**robust=t|f**Use a robust estimator to remove outliers before taking statistics again. [Default:**f**]. **mad=t|f**- Also compute the Mean Absolute Deviation (MAD).
[Default:
**f**]. **dual=t|f**- Dual pass over the data, subtracting the mean. This can be important if the data have a large average value with a small spread around the mean. Notice that in this mode the mean is subtracted but never added to the min/max/mean/median etc. Default: f
**nsigma=**- Delete points more than nsigma times the sigma from the mean. After each point sigma and mean are recomputed. [Default: none]
**sort=***sort_mode*

http://www.star.bris.ac.uk/~mbt/stilts/

http://arxiv.org/ps/0807.4820 (*choosing the binning for a histogram*)

~/src/kernel/tabsources

xx-mar-88V1.0: createdPJT 15-Apr-88V1.1: higher order moments, Y scalePJT 1-jun-88V2.0: new name, code samePJT 28-oct-88V2.0a: updated doc + labels plotting donePJT 13-nov-93V2.7: added gaussian model + residualsPJT 11-jul-96V2.8: log scale is now 10-based, not ePJT 12-apr-97V3.0: added cumulative optionPJT 24-apr-98V3.0a: fix median calculation for restricted rangePJT 22-dec-99V3.1a: optional median, fix N=1 reporting bugPJT 24-jan-00documentation updated with programPJT 7-jun-013.2: added nsigma, corrected man page optionsPJT 7-may-034.0: multiple columns allowedPJT 28-jan-055.0: separate xmin/xmax=, added sort=, fix median if nsigmaPJT 1-jun-106.0: bins= now allowed to have manual edgesPJT 22-aug-126.2: added torben= option for fast large-N medianPJT 16-jan-146.4: added mad=PJT