Thompson and Fearn (ref 3) - Example Data

HomTest-O-Matic v2.00. May 2002

N in fish food Refs

Best s as RSD 8.6% Based on International Protocol 1. Thompson, M. and Wood, R.,
Target SD as CV 10.0% Best s as RSD 6.9% Based on Fearn and Thompson 2. Thompson, M. and Lowthian, P
uHom 0.005 3. Thompson, M. and Fearn, T., (2

JAR A B A-B A+B D2 CochranCrit

1 2.49 2.39 0.11 4.88 0.0111 -
16 2.54 2.32 0.23 4.86 0.0520 0.9985
23 2.30 2.33 -0.04 4.63 0.0012 0.9669
39 2.30 2.44 -0.14 4.74 0.0197 0.9065
45 2.33 2.42 -0.09 4.75 0.0077 0.8412
51 2.56 2.40 0.16 4.96 0.0249 0.7808 Intermediate calculations

70 2.44 2.53 -0.09 4.96 0.0077 0.7271 s2all 0.005

74 2.51 2.47 0.04 4.98 0.0012 0.6798 s2
an 0.007
88 2.35 2.47 -0.12 4.82 0.0151 0.6385 s2sam 0.000
93 2.49 2.47 0.02 4.96 0.0003 0.6020 F1 1.88
0.00 0.00 0.0000 0.5700 F2 1.01
0.00 0.00 0.0000 0.5410 c 0.02
Replicatates 10 0.1409
GRAND MEAN (CFR) 2.43 Fearn and Thompson Tests
Value Critical Result
Cochran 0.37 0.60 Pass
F-TEST RESULTS Sa/s 0.35 0.5 Pass
F= 1.01 Fcrit = 3.02 s2sam 0.000 0.02 Pass

International Protocol Tests

Sa = 0.08 (Estimate of analytical standard deviation from ANOVA)
Ss = 0.005 (Estimate of sampling standard deviation from ANOVA)
Stotal = 0.08 (Sample standard deviation of all results)
s = 0.24 (Target standard deviation; 10%)
Ss/s = 0.02 (Is the stuff homogeneous?)
Sa/s = 0.35 (Is the method precise enough?) Critical value = 0.4
Stotal/s = NA (Use this test when F<1) Critical value = 0.4
Stotal/root6 = NA (Use this as an estimate of Ss for uncertainty purposes when F<1)

Anova: Single Factor

SUMMARY
Groups Count Sum Average Variance
Row 1 2 4.8755892 2.43779459 0.0055365
Row 2 2 4.8580511 2.42902554 0.0259909
Row 3 2 4.6300559 2.31502795 0.0006152
Row 4 2 4.7352844 2.36764222 0.0098427
Row 5 2 4.7528225 2.37641127 0.0038448
Row 6 2 4.9632796 2.48163981 0.0124572
Row 7 2 4.9632796 2.48163981 0.0038448
Row 8 2 4.9808177 2.49040886 0.0006152
Row 9 2 4.8229749 2.41148745 0.0075358
Row 10 2 4.9632796 2.48163981 0.0001538

ANOVA
Source of Variation SS df MS F P-value F crit
Between Gr0.0639161 9 0.00710179 1.0082484 0.49048 3.0203829471
Within Grou0.0704369 10 0.00704369

Total 0.134353 19

K:/acsl/aqas/pest/01-09/homog [Link] PMc 04/10/2023 [Link]

International Protocol Tests

F>1 Y F>F
F>F Y ss/s>0.4 Y Pass
c c
1.01 3.02 0.02

N N N

Stot/s<0.4 sa/s<0.4 Fail

NA 0.35
Y

N Y N Pass

Fail Pass Fail

Target 10.0 uHom = 0.292

CV %
Thompson and Fearn Tests

Components of Uncertainty
S2total = 0.007071
s2an 0.007
s2sam 0.000
u2hom
u hom
2
0.000025
Cochran
s2sam

s2an
Pass
S2total =
sa/s s2sam

K:/acsl/aqas/pest/01-09/homog [Link] PMc 04/10/2023 [Link]

AMC Kernel density ver 1.0e
Fixed h: .029
Plot Range
Start: -0.047
End: 0.5116
Step: 0.002191

Kernel estimate
X Y
-0.047 0.005328
-0.044809 0.006705
-0.042619 0.008393
-0.040428 0.010448
-0.038238 0.012934
-0.036047 0.015925
-0.033856 0.0195
-0.031666 0.023747
-0.029475 0.028762
-0.027285 0.034646
-0.025094 0.041507
-0.022904 0.049457
-0.020713 0.058609
-0.018522 0.069079
-0.016332 0.080977
-0.014141 0.094412
-0.011951 0.10948
-0.00976 0.126268
-0.007569 0.144844
-0.005379 0.165258
-0.003188 0.187534
-0.000998 0.211669
0.001193 0.237625
0.003384 0.265332
0.005574 0.294681
0.007765 0.325523
0.009955 0.357667
0.012146 0.390883
0.014336 0.424901
0.016527 0.459411
0.018718 0.494074
0.020908 0.528517
0.023099 0.562347
0.025289 0.595156
0.02748 0.626524
0.029671 0.656036
0.031861 0.683284
0.034052 0.707879
0.036242 0.729462
0.038433 0.747709
0.040624 0.762341
0.042814 0.773132
0.045005 0.779912
0.047195 0.782575
0.049386 0.781079
0.051576 0.775449
0.053767 0.765774
0.055958 0.752209
0.058148 0.734965
0.060339 0.714311
0.062529 0.690559
0.06472 0.664065
0.066911 0.635213
0.069101 0.604413
0.071292 0.572084
0.073482 0.538653
0.075673 0.504544
0.077864 0.470165
0.080054 0.43591
0.082245 0.402146
0.084435 0.369213
0.086626 0.337416
0.088816 0.307028
0.091007 0.278285
0.093198 0.25139
0.095388 0.226511
0.097579 0.203784
0.099769 0.183318
0.10196 0.165198
0.104151 0.149487
0.106341 0.136235
0.108532 0.125476
0.110722 0.117242
0.112913 0.111558
0.115104 0.108452
0.117294 0.107954
0.119485 0.110104
0.121675 0.114948
0.123866 0.122544
0.126056 0.132964
0.128247 0.14629
0.130438 0.162619
0.132628 0.182062
0.134819 0.204741
0.137009 0.230792
0.1392 0.260365
0.141391 0.293618
0.143581 0.330725
0.145772 0.371866
0.147962 0.417237
0.150153 0.467041
0.152344 0.521494
0.154534 0.580823
0.156725 0.645266
0.158915 0.715069
0.161106 0.790493
0.163296 0.871806
0.165487 0.959282
0.167678 1.053204
0.169868 1.153853
0.172059 1.26151
0.174249 1.376448
0.17644 1.498925
0.178631 1.629176
0.180821 1.767407
0.183012 1.913783
0.185202 2.068417
0.187393 2.231365
0.189584 2.40261
0.191774 2.582056
0.193965 2.769519
0.196155 2.964717
0.198346 3.167268
0.200536 3.376685
0.202727 3.592375
0.204918 3.813641
0.207108 4.039686
0.209299 4.269618
0.211489 4.502465
0.21368 4.737181
0.215871 4.97266
0.218061 5.207754
0.220252 5.441286
0.222442 5.672066
0.224633 5.898903
0.226824 6.120624
0.229014 6.336084
0.231205 6.544172
0.233395 6.743827
0.235586 6.934032
0.237776 7.113828
0.239967 7.282307
0.242158 7.438612
0.244348 7.581935
0.246539 7.711514
0.248729 7.826628
0.25092 7.926592
0.253111 8.01076
0.255301 8.07852
0.257492 8.129298
0.259682 8.162565
0.261873 8.177839
0.264064 8.174702
0.266254 8.152806
0.268445 8.111894
0.270635 8.051809
0.272826 7.972514
0.275016 7.874106
0.277207 7.756831
0.279398 7.621094
0.281588 7.467471
0.283779 7.29671
0.285969 7.109733
0.28816 6.907633
0.290351 6.691662
0.292541 6.463222
0.294732 6.223839
0.296922 5.97515
0.299113 5.718869
0.301304 5.456765
0.303494 5.19063
0.305685 4.922251
0.307875 4.65338
0.310066 4.385707
0.312256 4.120835
0.314447 3.860256
0.316638 3.605335
0.318828 3.357297
0.321019 3.11721
0.323209 2.885989
0.3254 2.664384
0.327591 2.452992
0.329781 2.252254
0.331972 2.062471
0.334162 1.883808
0.336353 1.716312
0.338544 1.559926
0.340734 1.414498
0.342925 1.279805
0.345115 1.155557
0.347306 1.041418
0.349496 0.937017
0.351687 0.841955
0.353878 0.755821
0.356068 0.678195
0.358259 0.608655
0.360449 0.546786
0.36264 0.492176
0.364831 0.444426
0.367021 0.403145
0.369212 0.367953
0.371402 0.338477
0.373593 0.314354
0.375784 0.295225
0.377974 0.280735
0.380165 0.27053
0.382355 0.264255
0.384546 0.261552
0.386736 0.262059
0.388927 0.265408
0.391118 0.271227
0.393308 0.279139
0.395499 0.288761
0.397689 0.29971
0.39988 0.311602
0.402071 0.324057
0.404261 0.336704
0.406452 0.349181
0.408642 0.361145
0.410833 0.372272
0.413024 0.382264
0.415214 0.390856
0.417405 0.397815
0.419595 0.402948
0.421786 0.406101
0.423976 0.407166
0.426167 0.406079
0.428358 0.402819
0.430548 0.397411
0.432739 0.389921
0.434929 0.380452
0.43712 0.369145
0.439311 0.356168
0.441501 0.341716
0.443692 0.326004
0.445882 0.309257
0.448073 0.291712
0.450264 0.273603
0.452454 0.255164
0.454645 0.236618
0.456835 0.218173
0.459026 0.200024
0.461216 0.182342
0.463407 0.165279
0.465598 0.148961
0.467788 0.13349
0.469979 0.118946
0.472169 0.105384
0.47436 0.092837
0.476551 0.081319
0.478741 0.070824
0.480932 0.061333
0.483122 0.052812
0.485313 0.045216
0.487504 0.038492
0.489694 0.032582
0.491885 0.027422
0.494075 0.022948
0.496266 0.019095
0.498456 0.015798
0.500647 0.012996
0.502838 0.01063
0.505028 0.008646
0.507219 0.006992
0.509409 0.005622
0.5116 0.004495
Calculation Robust value:
C. 2.1 Median, nIQR

No. Atrazine in drinking water (mg L-1)

1 0.0400
2 0.0550
3 0.1780
4 0.2020
5 0.2060
6 0.2270
7 0.2280
8 0.2300
9 0.2300
10 0.2350
11 0.2360
12 0.2370
13 0.2430
14 0.2440
15 0.2450
16 0.2555
17 0.2600
18 0.2640
19 0.2670
20 0.2700
21 0.2730
22 0.2740
23 0.2740
24 0.2780
25 0.2811
26 0.2870
 27 0.2870
28 0.2880
29 0.2890
30 0.2950
31 0.2960
32 0.3110
33 0.3310
34 0.4246
Average
SD
Median xpt
Q1 Quatile 1
Q3 Quatile 3
IQR (Q3-Q1) InterQuartile Range
Normalised IQR spt 0.7413*IQR
Standard uncertainty
p

Kernel Density plot:

spt = 15% From previous rounds of FAPAS PT

Optimum Bandwidth = 0.75 spt =

rtile Range
0.7413*IQR
Robust mean and SD: Algorithm A Atrazine in drinking water (mg L -1)

Iteration 0 |xi-x*| 1 2 3 4 5 6
d = 1.5 s* ---
x* - d ---
x* + d ---
34 0.425
33 0.331
32 0.311
31 0.296
30 0.295
29 0.289
28 0.288
26 0.287
27 0.287
25 0.281
24 0.278
22 0.274
23 0.274
21 0.273
20 0.270
19 0.267
18 0.264
17 0.260
16 0.256
15 0.245
14 0.244
13 0.243
12 0.237
11 0.236
10 0.235
8 0.230
9 0.230
7 0.228
6 0.227
5 0.206
4 0.202
3 0.178
2 0.055
1 0.040
New x*
New s*
p
% Different of Robust SD =

Average Conclusion for assigned value:

SD Robust mean xpt
Original %RSD Robust SD spt
Final %RSD Standard uncertainty
7 8 9 10 11 12 13 14 15
Outlier remove:

No. Atrazine in drinking water (mg L-1)

1 0.0400
2 0.0550
3 0.1780
4 0.2020
5 0.2060
6 0.2270
7 0.2280
8 0.2300
9 0.2300
10 0.2350
11 0.2360
12 0.2370
13 0.2430
14 0.2440
15 0.2450
16 0.2555
17 0.2600
18 0.2640
19 0.2670
20 0.2700
21 0.2730
22 0.2740
23 0.2740
24 0.2780
25 0.2811
26 0.2870
27 0.2870
 28 0.2880
29 0.2890
30 0.2950
31 0.2960
32 0.3110
33 0.3310
34 0.4246
Average xpt
SD spt
Standard uncertainty
n (p)
Q/Hampel method
The method known as Q/Hampel uses the Q method for the calculation of the robust standard deviation s* together

Available online: [Link]

Lab no. Result

1 0.040 Mean (Hampel): 
2 0.055 Standard deviation (Q-Method):
3 0.178 Standard uncertainty
4 0.202
5 0.206
6 0.227
7 0.228
8 0.230
9 0.230
10 0.235
11 0.236
12 0.237
13 0.243
14 0.244
15 0.245
16 0.256
17 0.260
18 0.264
19 0.267
20 0.270
21 0.273
22 0.274
23 0.274
24 0.278
25 0.281
26 0.287
27 0.287
28 0.288
29 0.289
30 0.295
31 0.296
32 0.311
33 0.331
34 0.425
n=
deviation s* together with the Hampel estimator for the calculation of the robust location parameter x*.
Summary of robust assigned value and standard uncertainty

Approach Robust mean/median Robust SD p

Robust: Median, nIQR

Robust: Algorithm A (x*, s*)

Arithmetic, outliers removed

Arithmetic, outliers included

Robust: Algorithm A, SDp from previous round

Target CV (%)

u(xpt) %CV
Horwitz equation
Case I: mg/L
Atrazine in drinking water 0.2570 mg/L

Mass fraction (c) =

= 1

sR = 2
=

SD = mg/L 3
or RSD = %

Case II: mg/kg

Melamine in milk powder 0.0190 mg/kg

Mass fraction (c) =

=

sR =
=
= mg/kg
or = %

Case III: g/100g

Fat in milk powder 18.0 g/100g

Mass fraction (c) =

=

sR =
=
= g/100g
or = %
 Equation used Consider c value

sR = 0.22 x c c < 1.2 x 10-7

sR = 0.02 x c0.8495 1.2 x 10-7 < c < 0.138

sR = 0.01 x c0.5 c > 0.138

Example of c value:

Values Mass fraction (c) Scientific number Select

Fat 15 g/100g

Pb 9 mg/kg

Dioxin 8 ug/kg
 0.8495

0.000001195

1.8610328E-07

0.1861032816