Peirce's Criterion for Outlier Removal

Versione 1.0.3 (2,04 KB) da cr
Apply Peirce's criterion for outlier detection for arbitrary number of observations 3 or more using Gould's implementation.
57 download
Aggiornato 5 gen 2023

Visualizza la licenza

Apply Peirce's criterion for outlier detection for arbitrary number of observations (≥3). Gould's implementation of Peirce's method is used. This implementation doesn't use lookup tables unlike the other submission and computes the R ratios using Gould's formulae. Hence the application is NOT limited to 60 observations -the size of the limited lookup table from [3] ROSS' paper.
Usage
Input DELTAS is a vector of the observations (or fit errors). Optional second input NUM_UNKNOWS is the number of "unknown quantities". Default value is 1 which is generally the case as argued by the authors. The function returns linear indices of outliers (OUTLIERINDICES), linear indices of observations which are found to be worthy of keeping (INLIERINDICES), and actual observations that are flagged as outliers.
Example
This example is from Peirce's original 1852 publication. Observations 1.01 and -1.4 are flagged as outliers.
z = [-0.3 0.48 0.63 -0.22 0.18 -0.44 -0.24 -0.13 -0.05 0.39 1.01 0.06 -1.4 0.2 0.1];
>> [outlierIdx,inlierIdx,outlierObs] = peirces(z)
outlierIdx =
11 13
inlierIdx =
1 2 3 4 5 6 7 8 9 10 12 14 15
outlierObs =
1.01 -1.4
References
B.A.Gould (1855), ON PEIRCE’S CRITERION FOR THE REJECTION OF DOUBTFUL OBSERVATIONS, WITH TABLES FOR FACILITATING ITS APPLICATION. The Astronomical Journal, No.83, Vol. IV.
B.PEIRCE (1852), CRITERION FOR THE REJECTION OF DOUBTFUL OBSERVATIONS. The Astronomical Journal, No.45, Vol. II.
S.M.ROSS, (2003) Peirce's criterion for the elimination of suspect experimental data. Journal of Engineering Technology.

Cita come

Chandrakanth R. Terupally, 2021. Peirce's Criterion for Outlier Rejection (https://www.mathworks.com/matlabcentral/fileexchange/<...>), MATLAB Central File Exchange. Retrieved November 3, 2021.

Compatibilità della release di MATLAB
Creato con R2017a
Compatibile con R2017a e release successive
Compatibilità della piattaforma
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Versione Pubblicato Note della release
1.0.3

Updated Title and text only.

1.0.2

Uploaded image.

1.0.1

Avoided computing Nth root of Q and raising Q to N.

1.0.0