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. It is returned as the bad value should there be only one or no valid data values.
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: