Table of Contents

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

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

- The TriMean is reported, which is a robust statistic for the central
- tendency, the weighted average of Median and two quartiles: (Q1 + 2*Q2 + Q3) / 4.

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, see *tabstat(1NEMO)*
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]
**qac=t|f**- If selected, a QAC listing is selected. This gives the mean, rms, min and max. [Default: f]
**sort=***sort_mode*- Default: qsort. If flogger is enabled additional sorting modes are available: bubble, heap, insert, merge, quick, shell.
**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.

% tabgen - 1000 2 | tabmath - - selfie=’range(%1,0.0,0.5)’ | tabhist - 2where column 2 is only used where column 1 is between 0 and 0.5.

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

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

Teuben

~/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 14-nov-20217.4: added qac=PJT