ГОСТ Р ИСО 5725-5-2002 сторінка 10 з 12
uU = 1,
получаем
(w*)2 = 1,0972 ´ 0,2495 + (1,097 ´ 1,645 w*)2/9,
что дает решение (если предположение, что uU = 1, корректно):
w* = 0,69 % креозота.
Можно затем подтвердить, что это значение дает ψ = 1,645 ´ 0,69 = 1,14, как и предполагалось, только w9* превышает ψ, и с последующей заменой w9* на 1,14 получаем снова w* = 0,63 ´ 1,097 = 0,69; значит, представленное решение правильно.
Следовательно, оценка стандартного отклонения повторяемости равна
sr = 0,69/ = 0,49 % креозота.
Это значение лежит между двумя оценками, полученными в 6.5.2 и 6.5.3.
6.5.5 Следующим этапом в анализе является получение робастной оценки стандартного отклонения средних значений в элементах. Применяя алгоритм А к средним значениям, получим результаты, представленные в таблице 26, где средние значения в элементах рассортированы в порядке возрастания. Из четырех итераций, представленных в таблице, ясно, что робастными значениями являются х* = 20,412 и s* ≈ 1,1 и что только два экстремальных средних значения в элементах (х1* = 17,570, х9* = 24,140) отличаются от х* более чем на j. Если бы расчеты выполнялись на компьютере, процесс мог бы быть продолжен, пока изменения в значениях х* и s* от одной итерации до следующей стали бы минимальными.
Таблица 26 - Пример 4. Применение Алгоритма А к средним значениям в элементах (% креозота)
Итерация
01)
1
2
3
4
j
-
1,424
1,478
1,514
1,539
х* - j
-
18,876
18,909
18,893
18,872
x* + j
-
21,724
21,865
21,921
21,950
x1*
17,570
18,876
18,909
18,893
18,872
x2*
19,500
19,500
19,500
19,500
19,500
x3*
20,100
20,100
20,100
20,100
20,100
x4*
20,155
20,155
20,155
20,155
20,155
x5*
20,300
20,300
20,300
20,300
20,300
x6*
20,705
20,705
20,705
20,705
20,705
x7*
20,940
20,940
20,940
20,940
20,940
x8*
21,185
21,185
21,185
21,185
21,185
x9*
24,140
21,724
21,865
21,921
21,950
Среднее
20,511
20,387
20,407
20,411
20,412
Стандартное отклонение
1,727
0,869
0,890
0,905
0,916
Новые x*
20,3002)
20,387
20,407
20,411
20,412
Новые s*
0,9492)
0,985
1,009
1,026
1,039
1) Значения получены из таблицы 24 после перестановки в порядке возрастания.
2) Получены по формулам (56) и (57).
При расчете вручную аналитик должен использовать прямой метод, описанный в 6.2.6, например иL = иU = 1.
Это дает х′ = 20,412 и s′ = 0,573 % креозота.
Отсюда из уравнений (62) и (63)
(s*)2 = 6 ´ (0,573)2/[8/1,1342 - 1,52 (9 + 9 - 4)/7]
получаем s* = 1,070 % креозота и х* = х′ = 20,412 % креозота.
Можно затем подтвердить, что значение s* дает j = 1,605 (тогда как предполагалось, что только x1* и х9* отличаются от х* = 20,412 более чем на j) и что замена x1* на 18,807 и x9* на 22,017 дает новое значение х* = 20,412, а новое значение s* = 0,944 ´ 1,134 = 1,070, так что представленное решение является правильным.
Оценку межлабораторного стандартного отклонения проводят по формуле (72):
% креозота,
а оценку стандартного отклонения воспроизводимости - по формуле (74):
% креозота.
Снова это значение располагается между двумя оценками, полученными в 6.5.2 и 6.5.3.
6.6 Формулы. Робастный анализ для отдельного уровня в эксперименте по модели с разделенными уровнями 6.6.1 Применительно к модели с разделенными уровнями робастная оценка стандартного отклонения повторяемости sr для отдельного уровня может быть получена обработкой данных о расхождениях в элементах на определенном уровне по Алгоритму А с нахождением робастного значения s* из уравнения (61), а затем определением sr по формуле
sr = s* . (75)
6.6.2 Робастная оценка стандартного отклонения средних значений sy в элементах для уровня может быть получена применением Алгоритма А снова к средним значениям в элементах для определенного уровня, нахождением робастного значения s* из уравнения (61), а потом получением sy, с использованием равенства
sy = s*. (76)
Для оценки стандартного отклонения воспроизводимости на определенном уровне модели можно использовать формулы, приведенные в 4.5.6.
6.7 Пример 5. Робастный анализ для отдельного уровня в эксперименте по модели с разделенными уровнями 6.7.1 Данные примера 1 в 4.8 содержали несколько квазивыбросов и один выброс (см. таблицу 8). Кроме того, на рисунке 3 видна отрицательная систематическая погрешность в результатах лаборатории № 5. Если аналитик не может выявить источники этих аномалий, он попадает в трудное положение при принятии решения, какие данные следует исключить из расчетов стандартных отклонений повторяемости и воспроизводимости. Для иллюстрации результатов робастного анализа здесь использованы данные уровня 14 (см. таблицу 4).
6.7.2 Для получения робастной оценки стандартного отклонения повторяемости расхождений в элементах применяют Алгоритм А (см. таблицу 5), что приводит к результатам, показанным в таблице 27, в которой расхождения в элементах рассортированы в порядке возрастания. Из четырех итераций, представленных в таблице, видно, что робастные значения равны х* = 8,29, s* = 0,36, и что только х9* отличается от х* более чем на j.
Применяя метод, описанный в 6.2.6 при uL = 0 и uU = 1, получаем
х′ = 8,219 и s′ = 0,257 % протеина.
Тогда уравнения (62) и (63) в 6.2.6 можно записать в виде
х* = 8,219 + 1,5 ´ s*/8
и
(s*)2 = 7 ´ (0,257)2/[8/1,1342 - 1,52(0 + 9 - 0)/8],
что дает s* = 0,354 % протеина,
а, используя уравнение (75), получим sr = 0,354/ = 0,250 % протеина.
Таблица 27 - Пример 5. Применение Алгоритма А к расхождениям в элементах (% протеина)
Итерация
0
1
2
3
4
j
-
0,53
0,56
0,55
0,54
х* - j
-
7,85
7,74
7,74
7,75
х* + j
-
8,91
8,86
8,84
8,83
x1*
7,81
7,85
7,81
7,81
7,81
x2*
7,93
7,93
7,93
7,93
7,93
x3*
8,13
8,13
8,13
8,13
8,13
x4*
8,14
8,14
8,14
8,14
8,14
x5*
8,38
8,38
8,38
8,38
8,38
x6*
8,40
8,40
8,40
8,40
8,40
x7*
8,44
8,44
8,44
8,44
8,44
x8*
8,52
8,52
8,52
8,52
8,52
x9*
9,31
8,91
8,86
8,84
8,83
Среднее
8,340
8,300
8,290
8,288
8,287
Стандартное отклонение
0,436
0,326
0,322
0,317
0,315
Новые x*
8,3801)
8,300
8,290
8,288
8,287
Новые s*
0,3561)
0,370
0,365
0,359
0,357
1) Получено по формулам (56) и (57).
Робастное среднее значение для расхождений в элементах составляет
х* = 8,219 + 1,5 ´ 0,354/8 = 8,285 % протеина.
При этих значениях х* и s*
j = 1,5 ´ 0,354 = 0,531.
Тогда x* - j = 7,754 и х* + j = 8,816 % протеина.
Это подтверждает, что в расчетах х* и s* только х9* было вне этих пределов. Можно сделать заключение, что это тот самый случай, когда действительно найдено правильное решение.
6.7.3 Применение Алгоритма А к средним значениям в элементах (из таблицы 6) дает результаты, представленные в таблице 28, в которой средние значения в элементах расположены в порядке возрастания. Ситуация подобна представленной в таблице 26, а именно х1* и х9* отличаются более чем на j от х* и х* устремляется к среднему значению от х2* до х8*, равному 85,486. Применяя вновь метод из 6.2.6, но со значениями ul = иU = 1, получим, что среднее значение и стандартное отклонение от х2* до х8* составят:
х΄ = 85,486 и s΄ = 0,209.
Значит, на основе уравнения (63) может быть получено s* из равенства
(s*)2 = 6 ´ (0,209)2/[8/1,1342 - 1,52(9 + 9 - 4)/7],
откуда s* = 0,390 % протеина.
Теперь можно вычислить х* по формуле (62) в 6.2.6, что дает x* = 85,486 % протеина. Для контроля правильности решения, рассчитывают
j = 1,5 ´ 0,390 = 0,585, х* - j = 84,901, х* + j = 86,071 % протеина.
Таблица 28 - Пример 5. Применение Алгоритма А к средним значениям в элементах (% протеина)
Итерация
0
1
2
3
4
j
-
0,446
0,492
0,519
0,537
х* - j
-
85,104
85,009
84,971
84,950
х* + j
-
85,996
85,993
86,009
86,024
x1*
84,525
85,104
85,009
84,971
84,950
x2*
85,140
85,140
85,140
85,140
85,140
x3*
85,345
85,345
85,345
85,345
85,345
x4*
85,385
85,385
85,385
85,385
85,385
x5*
85,550
85,550
85,550
85,550
85,550
x6*
85,575
85,575
85,575
85,575
85,575
x7*
85,660
85,660
85,660
85,660
85,660
x8*
85,750
85,750
85,750
85,750
85,750
x9*
86,170
85,996
85,993
86,009
86,024
Среднее
85,456
85,501
85,490
85,487
85,487
Стандартное отклонение
0,453
0,289
0,305
0,316
0,324
Новые х*
Побачили розбіжність з офіційним текстом?
Дивіться також
Сборник ГОСТ ЕСКД
Збірник із 154 стандартів