CCG_WMD2x - Estimates the mean of a number of normally distributed data values, some of which may be corrupt

Description:

Not available. Remember to look in EDRS.

Invocation:

CALL CCG_WMD2x( EL, X, W, PBAD, NITER, TOLL, NSIGMA, XMODE, SIGMA, USED, NUSED, STATUS )

Arguments:

EL = INTEGER (Given)

The number of data values.

X( EL ) = ? (Given)

An array of data values.

W( EL ) = ? (Given)

An array of data weights for each data value. The weights are inversely proportional to the square of the relative errors on each data point.

PBAD = REAL (Given)

An estimate of the probability that any one data point will be corrupt. (This value is not critical.)

NITER = INTEGER (Given)

The maximum number of iterations required.

TOLL = REAL (Given)

The absolute accuracy required in the estimate of the mean. Iterations cease when two successive estimates differ by less than this amount.

NSIGMA = REAL (Given)

The sigma level to reject data values at.

XMODE = DOUBLE PRECISION (Returned)

The estimate of the uncorrupted mean.

SIGMA = DOUBLE PRECISION (Returned)

An estimate of the uncorrupted normalised standard deviation of the data points. An estimate of the standard deviation of any one point is: SIGMA / SQRT( W ) where W is its weight.

USED( EL ) = LOGICAL (Returned)

If a value is not rejected then its corresponding used element will be set true.

NUSED = INTEGER (Returned)

Number of the input data values which are actually used.

STATUS = INTEGER (Returned)

The global status.

Notes:

• There is a routine for real and double-precision: replace "x" in the routine name by D or R as appropriate. The X and W arguments supplied to the routine must have the data type specified.