**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). One rule
of thumb is that the number of bins be the square root of the number of
values. [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**]. **norm=t|f**- Normalize the cumulative histogram to 1. By default the maxcount is used.
[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***pyplot=**- If given, it will be the filename where a template python script that can serve as starting point for more elaborate plotting. Default: none.

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 8-jan-20207.0: added pyplot=PJT 2-mar-20207.1: added norm=PJT