Description:
This routine calculates the mean and standard deviation of a double array. If a positive value of CLIP
is given, the data set will be clipped at CLIP sigma repeatedly until all points are within CLIP sigma.
Mean and standard deviation are those values derived after iterative clipping.
Invocation
CALL
SCULIB_STATD(N_POS, SCUDATA, SCUQUAL, BADBIT, NGOOD, MEAN, MEDIAN, SUM, SUMSQ, STDEV,
QSORT, STATUS)
Arguments
N_POS = _INTEGER (Given)
Number of jiggle positions in data set
CLIP = _REAL (Given)
Number of sigma to clip iteratively
SCUDATA(N_POS) = _DOUBLE
(Given)
The data
SCUQUAL(N_POS) = _BYTE (Given)
The data quality
BADBIT
= _BYTE (Given)
Bad bit mask
NGOOD = _INTEGER (Returned)
Number of good
pixels in mean
MEAN = _DOUBLE (Returned)
Mean value
MEDIAN = _DOUBLE
(Returned)
Median of good data
SUM = _DOUBLE (Returned)
Data sum
SUMSQ
= _DOUBLE (Returned)
Sum of squares
STDEV = _DOUBLE (Returned)
Standard
deviation
QSORT = _DOUBLE (Returned)
Sorted data set
STATUS = INTEGER (Given and
Returned)
Global Status value
Copyright
Copyright ©1995,1996,1997,1998,1999 Particle
Physics and Astronomy Research Council. All Rights Reserved.
Implementation Status:
-
Deals with bad pixels
-
Uses PDA sort routine