Формулы, таблицы и рисунки в разделе 6 и приложении B ГОСТ Р ИСО 5725-1 могут быть использованы при выборе числа лабораторий, проб и параллельных определений, но с модификациями, изложенными в 5.3.2 до 5.3.5.
5.3.2. Неопределенность оценки стандартного отклонения повторяемости, полученной из эксперимента на гетерогенном материале, может быть оценена расчетом величины
вместо определенной равенством (9) ГОСТ Р ИСО 5725-1. Однако вышеприведенная формула может быть получена заменой p в уравнении (9) ГОСТ Р ИСО 5725-1 на p' x g. Значит, на рисунке B.1 данные для повторяемости под
Примечание. Формулы (16) для
5.3.3. Неопределенность оценки стандартного отклонения воспроизводимости, полученной из эксперимента на гетерогенном материале, может быть оценена вычислением величины
вместо определенной уравнением (10) ГОСТ Р ИСО 5725-1,
где
; ; ; . (18)Величины Ф и
5.3.4. Детальную организацию эксперимента с гетерогенным материалом осуществляют в соответствии с руководством, изложенным в разделах 5 и 6 ГОСТ Р ИСО 5725-2.
Подпункт 5.1.2 ГОСТ Р ИСО 5725-2 содержит требования для "группы из n испытаний" или "группы из n измерений" (например требование, что группа из n испытаний должна проводиться с соблюдением условий повторяемости). В эксперименте на гетерогенном материале эти требования относят к группе g x n испытаний в элементе, то есть ко всем испытаниям в одной лаборатории на одном уровне.
В эксперименте на гетерогенном материале число проб, которое должно быть приготовлено для каждого уровня, равно p' x g (то есть 2p' в обычном случае, когда g = 2). Важно разместить эти p' x g проб по лабораториям-участницам случайным образом.
5.4.1. Основная модель, использованная в настоящем стандарте, описана в 4.1.1 равенством (3). Для эксперимента с гетерогенным материалом эта модель принимает вид
. (19)Члены m, B и e имеют те же значения, как и в равенстве (3), но равенство (19) содержит особый член
Естественно полагать, что различие между пробами является случайной величиной, не зависящей от лаборатории, но оно может зависеть от уровня в эксперименте. Тогда член
. (20)5.4.2. В обычном случае с двумя пробами для лаборатории и двумя результатами измерений для пробы (g = n = 2), определяют:
a) среднее для пробы и расхождения между результатами испытаний для лаборатории i, уровня j и пробы t (t = 1 или 2)
; (22)b) среднее для элемента и различие между пробами для лаборатории i и уровня j
, (23)c) общее среднее и стандартное отклонение средних для элементов на уровне j
где суммирование осуществляют по всем лабораториям i = 1, 2, ..., p'.
5.5.1. В этом пункте детально рассматривают случай, когда для каждой лаборатории на каждом уровне приготовлены по две пробы и на каждой пробе получены два результата измерений (общий случай рассматривают в 5.9 и 5.10).
Группируют полученные данные в таблицу (см. таблицу 9). Каждая комбинация лаборатории и уровня образует "элемент" в этой таблице, содержащий четыре результата измерений.
Используя уравнения (21) - (26):
a) рассчитывают расхождения между результатами измерений и сводят их в таблицу (см. таблицу 10);
b) рассчитывают расхождения между пробами и сводят их в таблицу (см. таблицу 11);
c) рассчитывают средние для элементов и сводят их в таблицу (см. таблицу 12).
Записывают расхождения как положительные величины (то есть игнорируя знак).
Таблица 10
между результатами измерений в эксперименте
для гетерогенного материала
???????????????????????????????????????????????????????????????????????????
? Номер ? Номер ? Уровень ?
?лаборатории? пробы ???????????????????????????????????????????????????????
? ? ? 1 ? 2 ? ? j ? ? q ?
???????????????????????????????????????--------??????????--------??????????
? 1 ? 1 ? ? ? ? ? ? ?
? ? 2 ? ? ? ? ? ? ?
???????????????????????????????????????--------??????????--------??????????
? 2 ? 1 ? ? ? ? ? ? ?
? ? 2 ? ? ? ? ? ? ?
???????????????????????????????????????--------??????????--------??????????
| | | | | | | | |
???????????????????????????????????????--------??????????--------??????????
? i ? 1 ? ? ? ? ? ? ?
? ? 2 ? ? ? ? ? ? ?
???????????????????????????????????????--------??????????--------??????????
| | | | | | | | |
???????????????????????????????????????--------??????????--------??????????
? p' ? 1 ? ? ? ? ? ? ?
? ? 2 ? ? ? ? ? ? ?
???????????????????????????????????????--------??????????--------??????????
Таблица 11
между пробами в эксперименте для гетерогенного материала
???????????????????????????????????????????????????????????????????????????
? Номер ? Уровень ?
?лаборатории ??????????????????????????????????????????????????????????????
? ? 1 ? 2 ? ? j ? ? q ?
???????????????????????????????????---------???????????---------???????????
? 1 ? ? ? ? ? ? ?
???????????????????????????????????---------???????????---------???????????
? 2 ? ? ? ? ? ? ?
???????????????????????????????????---------???????????---------???????????
| | | | | | | |
???????????????????????????????????---------???????????---------???????????
? i ? ? ? ? ? ? ?
???????????????????????????????????---------???????????---------???????????
| | | | | | | |
???????????????????????????????????---------???????????---------???????????
? p' ? ? ? ? ? ? ?
???????????????????????????????????---------???????????---------???????????
Таблица 12
по элементам в эксперименте для гетерогенного материала
???????????????????????????????????????????????????????????????????????????
? Номер ? Уровень ?
?лаборатории ??????????????????????????????????????????????????????????????
? ? 1 ? 2 ? ? j ? ? q ?
???????????????????????????????????---------???????????---------???????????
? 1 ? ? ? ? ? ? ?
???????????????????????????????????---------???????????---------???????????
? 2 ? ? ? ? ? ? ?
???????????????????????????????????---------???????????---------???????????
| | | | | | | |
???????????????????????????????????---------???????????---------???????????
? i ? ? ? ? ? ? ?
???????????????????????????????????---------???????????---------???????????
| | | | | | | |
???????????????????????????????????---------???????????---------???????????
? p' ? ? ? ? ? ? ?
???????????????????????????????????---------???????????---------???????????
5.5.2. Если элемент в таблице 9 содержит менее четырех результатов измерений (например, по причине порчи проб или исключения данных после применения методов контроля наличия выбросов, описанных ниже), тогда:
a) либо используют формулы для общего случая, приведенные ниже;
b) либо игнорируют все данные в элементе.
Альтернатива a) является предпочтительной. Выбор b) - бросовые данные, допускает применение простых формул.
a) сумму квадратов расхождений между результатами измерений в графе j таблицы 10 (суммируют по p' лабораториям и двум пробам)
b) сумму квадратов расхождений между пробами в графе j таблицы 11 (суммируют все p' лабораторий)
c) среднее значение и стандартное отклонение средних для элементов в графе j таблицы 12 с использованием уравнений (25) и (26).
5.5.4. Используют таблицы 10 - 12 и статистические результаты, рассчитанные по 5.5.3, чтобы оценить данные на однородность и наличие выбросов, как описано в 5.6. Если какие-то данные исключают, пересчитывают статистические результаты.
5.5.5. Рассчитывают стандартные отклонения повторяемости
. (30)Если это дает
, (31)тогда устанавливают
. (32)Рассчитывают оценку стандартного отклонения
Примечания. 1. Может показаться интересным выполнить испытание на значимость, чтобы определить, является ли расхождение между пробами статистически значимым, однако это не является необходимой частью анализа. Некорректно использовать такое испытание, чтобы решить, можно ли пренебречь расхождением между пробами в анализе (так как результаты измерений в каждом элементе обрабатывают так, как если бы они все были получены на одной и той же пробе). Это внесло бы систематическую погрешность в оценку стандартного отклонения повторяемости, поскольку утверждение о том, что расхождение между пробами не является статистически значимым, не доказывает, что этим расхождением можно пренебречь.
2. В случае, описанном в 5.1.5 (когда имеются три фактора: "лаборатории", "испытания внутри лабораторий" и "параллельные определения при выполнении испытаний"), стандартные отклонения повторяемости и воспроизводимости должны рассчитываться по формулам:
, .Эти формулы применяют, когда результаты испытаний рассчитывают как среднее результатов двух определений.
5.5.6. Исследуют зависимость
5.6.1. При проверке данных на совместимость используют статистики h и k, описанные в 7.3.1 ГОСТ Р ИСО 5725-2.
Для контроля совместимости средних значений в элементах рассчитывают статистику h по формуле
. (34)Наносят статистические данные на график, чтобы показать, в каких лабораториях имеет место несовместимость, выстраивают данные по уровням, а также группируют их по лабораториям.
Для контроля совместимости расхождений между пробами рассчитывают статистику k по формуле
. (35)Наносят статистические данные на график, чтобы показать, в каких лабораториях имеет место несовместимость, выстраивают данные по уровням, а также группируют их по лабораториям.
Для контроля совместимости расхождений между результатами измерений рассчитывают статистику k по формуле
. (36)Наносят эти статистические данные на график, чтобы показать, в каких лабораториях имеет место несовместимость, выстраивают данные по уровням, а также группируют их по лабораториям.
Интерпретация графиков полностью описана в 7.3.1 ГОСТ Р ИСО 5725-2. Если лаборатория сообщает результаты с систематическими погрешностями, то для нее большинство данных по статистике h для средних значений в элементах на соответствующем графике будет большим и иметь одно направление. Если лаборатория не провела измерение внутри уровней с соблюдением условий повторяемости (и допустила наличие посторонних факторов, увеличивших расхождение между пробами), то на соответствующем графике для статистики k будут видны необычно большие статистические данные для расхождений между пробами. Если лаборатория имеет плохую повторяемость, это проявится на графике в виде необычно больших значений статистики k для расхождений между результатами измерений.
5.6.2. Исследуют данные эксперимента с точки зрения наличия квазивыбросов и выбросов по критериям Кохрена и Граббса, как это описано в 7.3.3 и 7.3.4 ГОСТ Р ИСО 5725-2.
Для проверки наличия квазивыбросов и выбросов в расхождениях между результатами измерений рассчитывают значения статистики Кохрена для каждого уровня j по формуле
, (37)где
Для пользования таблицей критических значений подпункта 8.1 ГОСТ Р ИСО 5725-2 следует в таблице найти ряд, соответствующий p = 2p', на краю слева и в графе n = 2.
Чтобы проверить наличие квазивыбросов и выбросов в расхождениях между пробами, рассчитывают значения статистики Кохрена для каждого уровня j по формуле
, (38)где
Для пользования таблицей критических значений в ГОСТ Р ИСО 5725-2 следует в таблице найти ряд, соответствующий p = p', на краю слева и в графе n = 2.
Чтобы проверить на наличие квазивыбросов и выбросов средние значения в элементах для каждого уровня j, рассчитывают по этим значениям статистику Граббса, как это показано в 7.3.4 ГОСТ Р ИСО 5725-2 (где s - это
Интерпретация этих проверок полностью описана в 7.3.2 ГОСТ Р ИСО 5725-2. В эксперименте на гетерогенном материале результаты этих проверок должны быть использованы в следующем порядке. Сначала нужно применить тест Кохрена к расхождениям между результатами измерений. Если на основе этого анализа решено, что расхождение между результатами измерений является выбросом и должно быть исключено, тогда оба результата измерений, которые дали выброс, должны быть исключены при расчетах стандартных отклонений повторяемости и воспроизводимости (но при этом другие результаты измерений в элементе должны быть оставлены). Далее применяют тест Кохрена к расхождениям между пробами и, наконец, - тесты Граббса к средним значениям в элементах. Если решено, что расхождение между пробами или среднее значение в элементе является выбросом и что результаты, которые стали источником таких выбросов, подлежат исключению, тогда все экспериментальные данные для соответствующих элементов исключают из расчетов стандартных отклонений повторяемости и воспроизводимости.
Рекомендации, предложенные в 4.7, в равной степени применимы к эксперименту на гетерогенном материале.
5.8.1. Агрегатированные частицы материалов (связанный цемент или битум), служащие для покрытия аэродромов и дорог, должны обладать определенной влаго- и морозостойкостью. Метод, который применяют для измерения этих их возможностей, - это испытание на прочность с использованием сульфата магния согласно BS 812-12 [6], при котором испытуемую навеску материала подвергают пропитке (в несколько циклов) в насыщенном растворе сульфата магния с последующей сушкой. Изначально навеску готовят из остатка на сите с отверстиями 10 мм после отсева. В процессе испытаний частицы измельчают, и результатом измерения является массовая доля от испытуемой навески, которая проходит через сито с отверстиями 10 мм. Высокий результат (свыше 10% до 20% по массе) означает агрегатное состояние с плохой прочностью.
5.8.2. Данные, представленные в таблице 13, были получены в эксперименте, в котором пары проб, отобранные от восьми образцов материала, были направлены в 11 лабораторий, и на каждой пробе были получены два результата измерений на прочность с применением сульфата магния. Пробы были массой около 100 кг (они использовались в ряде других испытаний), а испытуемые навески были массой около 350 г.
Таблица 13
с помощью сульфата магния, %
???????????????????????????????????????????????????????????????????????????
? Номер ? Номер пробы ? Уровень ?
? лаборатории ? ?????????????????????????????????????????
? ? ? 1 ? 2 ? 3 ? 4 ?
? ? ?????????????????????????????????????????
? ? ? Номер результата измерений ?
? ? ?????????????????????????????????????????
? ? ? 1 2 ? 1 2 ? 1 2 ? 1 2 ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 1 ?69,2 67,0?7,4 8,0 ?4,1 3,5 ?10,4 10,1?
? ? 2 ?69,7 71,7?6,6 5,7 ?10,5 13,1?13,9 13,8?
???????????????????????????????????????????????????????????????????????????
? 2 ? 1 ?66,5 64,1?1,9 2,1 ?3,0 3,2 ?8,7 6,7 ?
? ? 2 ?65,7 65,8?4,2 3,3 ?1,9 1,1 ?8,3 4,8 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 1 ?68,7 69,5?6,3 5,8 ?2,4 2,9 ?11,7 7,0 ?
? ? 2 ?67,7 77,7?9,7 5,3 ?2,1 3,3 ?7,9 12,0?
???????????????????????????????????????????????????????????????????????????
? 4 ? 1 ?77,5 75,3?2,0 3,6 ?2,4 1,4 ?9,4 7,1 ?
? ? 2 ?76,3 77,2?4,7 3,8 ?6,4 2,3 ?10,7 7,7 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 1 ?55,4 63,2?3,8 4,1 ?1,3 0,8 ?3,7 6,3 ?
? ? 2 ?65,9 54,7?2,1 3,1 ?0,7 1,7 ?3,3 3,7 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 1 ?64,8 70,9?8,4 6,1 ?6,0 9,7 ?16,5 12,3?
? ? 2 ?78,2 73,4?8,3 10,6?12,4 9,8 ?13,2 16,8?
???????????????????????????????????????????????????????????????????????????
? 7 ? 1 ?64,8 63,4?4,3 5,7 ?2,9 3,0 ?7,5 9,3 ?
? ? 2 ?67,0 63,4?7,7 3,9 ?4,3 6,4 ?11,1 8,3 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 1 ?64,9 68,4?4,4 2,8 ?1,3 2,8 ?5,7 6,8 ?
? ? 2 ?65,4 65,5?5,4 6,7 ?2,7 2,8 ?4,8 5,5 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 1 ? - - ? - - ?1,1 0,0 ?6,6 7,0 ?
? ? 2 ? - - ? - - ?0,7 3,7 ?4,9 6,3 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 1 ?57,0 57,7?3,3 0,4 ?2,1 2,4 ?5,5 5,8 ?
? ? 2 ?57,1 52,7?4,2 2,3 ?3,6 3,5 ?3,9 5,7 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 1 ?70,6 75,2?5,3 6,4 ?5,7 1,9 ?9,5 7,2 ?
? ? 2 ?77,9 68,2?3,5 7,1 ?1,4 3,0 ?8,1 7,4 ?
???????????????????????????????????????????????????????????????????????????
Окончание таблицы 13
???????????????????????????????????????????????????????????????????????????
? Номер ? Номер пробы ? Уровень ?
? лаборатории ? ?????????????????????????????????????????
? ? ? 1 ? 2 ? 3 ? 4 ?
? ? ?????????????????????????????????????????
? ? ? Номер результата измерений ?
? ? ?????????????????????????????????????????
? ? ? 1 2 ? 1 2 ? 1 2 ? 1 2 ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 1 ?8,9 7,4 ?31,1 28,5?38,7 41,7?4,2 4,1 ?
? ? 2 ?7,6 9,1 ?23,0 23,1?44,2 41,1?7,3 4,4 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 1 ?3,2 3,5 ?16,5 15,4?36,6 45,2?3,2 5,4 ?
? ? 2 ?2,8 4,0 ?10,3 12,8?43,2 40,5?1,7 2,5 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 1 ?4,4 6,1 ?24,3 16,7?38,9 43,1?3,7 7,7 ?
? ? 2 ?6,0 6,0 ?20,8 22,2?46,1 47,4?3,5 5,6 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 1 ?2,7 3,1 ?20,2 16,2?32,0 35,5?2,9 2,2 ?
? ? 2 ?2,3 2,9 ?20,0 11,9?26,5 35,7?3,2 2,3 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 1 ?1,3 1,4 ?13,8 15,1?36,7 39,5?1,1 1,2 ?
? ? 2 ?1,5 1,3 ?11,5 13,3?37,6 34,1?0,6 1,7 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 1 ?8,2 4,2 ?20,3 24,7?49,4 50,6?11,9 18,5?
? ? 2 ?3,7 4,6 ?21,0 18,9?48,2 52,4?14,9 8,1 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 1 ?3,1 5,5 ?27,2 23,3?38,9 29,9? - 1,7 ?
? ? 2 ?5,6 5,5 ?21,5 22,7?34,4 38,3?2,2 5,0 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 1 ?1,8 2,2 ?13,6 12,0?27,0 37,0?0,3 2,2 ?
? ? 2 ?4,0 4,0 ?15,6 16,7?39,7 34,6?3,6 3,7 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 1 ?3,8 3,8 ?17,7 17,1?33,4 33,1?1,8 2,0 ?
? ? 2 ?3,5 2,8 ?21,4 16,8?26,5 25,2?2,5 1,6 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 1 ?3,5 3,0 ?21,7 23,9?35,3 26,5?0,5 4,3 ?
? ? 2 ?3,2 3,5 ?27,0 32,5?18,0 18,2?2,0 2,1 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 1 ?3,5 2,5 ?11,0 18,4?27,0 33,5?5,1 3,9 ?
? ? 2 ?2,0 2,8 ?16,4 8,1 ?35,4 29,3?2,1 5,0 ?
???????????????????????????????????????????????????????????????????????????
5.8.3. Таблицы 14 - 16 показывают расхождения между результатами измерений, различия между пробами и средние значения в элементах, рассчитанные с использованием равенств (21) - (24), только для уровня 6 эксперимента.
Таблица 14
для уровня 6
???????????????????????????????????????????????????????????????????????????
?Номер лаборатории ? Номер пробы ?Расхождение между ? Статистика k ?
? ? ? результатами ? ?
? ? ? измерений, % ? ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 1 ? 2,6 ? 0,624 ?
? ? 2 ? 0,1 ? 0,024 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 1 ? 1,1 ? 0,264 ?
? ? 2 ? 2,5 ? 0,600 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 1 ? 7,6 ? 1,825 ?
? ? 2 ? 1,4 ? 0,336 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 1 ? 4,0 ? 0,960 ?
? ? 2 ? 8,1 ? 1,945 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 1 ? 1,3 ? 0,312 ?
? ? 2 ? 1,8 ? 0,432 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 1 ? 4,4 ? 1,056 ?
? ? 2 ? 2,1 ? 0,504 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 1 ? 3,9 ? 0,936 ?
? ? 2 ? 1,2 ? 0,288 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 1 ? 1,6 ? 0,384 ?
? ? 2 ? 1,1 ? 0,264 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 1 ? 0,6 ? 0,144 ?
? ? 2 ? 4,6 ? 1,104 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 1 ? 2,2 ? 0,528 ?
? ? 2 ? 5,5 ? 1,320 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 1 ? 7,4 ? 1,777 ?
? ? 2 ? 8,1 ? 1,945 ?
???????????????????????????????????????????????????????????????????????????
Таблица 15
???????????????????????????????????????????????????????????????????????????
? Номер лаборатории ? Расхождение между ? Статистика k ?
? ? пробами, % ? ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 6,75 ? 1,767 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 4,40 ? 1,152 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 1,00 ? 0,262 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 2,25 ? 0,589 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 2,05 ? 0,537 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 2,55 ? 0,668 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 3,15 ? 0,825 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 3,35 ? 0,877 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 1,70 ? 0,445 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 6,95 ? 1,819 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 2,55 ? 0,668 ?
???????????????????????????????????????????????????????????????????????????
Таблица 16
???????????????????????????????????????????????????????????????????????????
? Номер лаборатории ? Среднее элемента, % ? Статистика h ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 26,425 ? 1,475 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 13,750 ? -1,043 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 21,000 ? 0,397 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 17,075 ? -0,382 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 13,425 ? -1,108 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 21,225 ? 0,442 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 23,675 ? 0,929 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 14,475 ? -0,899 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 18,250 ? -0,149 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 26,275 ? 1,445 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 13,425 ? -1,108 ?
???????????????????????????????????????????????????????????????????????????
Подставляя в равенства (27) и (28) расхождения между результатами измерений из таблицы 14 и между пробами из таблицы 15, получаем
, . (общее среднее), .Так что, используя уравнения (29) - (33), для стандартных отклонений повторяемости и воспроизводимости и стандартного отклонения, которое измеряет расхождение между пробами, получим:
, , .Таблица 17 дает результаты расчетов по другим уровням.
Таблица 17
и стандартные отклонения, рассчитанные по данным всех
восьми уровней в таблице 13 (исключая элементы
с опущенными данными)
???????????????????????????????????????????????????????????????????????????
?Уровень j? Число ? Общее ?Суммы квадратов? Стандартные отклонения, % ?
? ?лабораторий?среднее? расхождений, ? ?
? ? p' ? y , % ? 2 ? ?
? ? ? j ? % ? ?
? ? ? ?????????????????????????????????????????????
? ? ? ? SS ? SS ? s ? s ? s ? s ?
? ? ? ? rj? Hj ? yj ? rj ? Rj ? Hj ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 11 ? 3,7 ?82,99 ?96,3725 ? 2,62 ? 1,37 ? 2,56 ? 1,85 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 11 ? 4,0 ?34,70 ?11,2550 ? 1,88 ? 0,89 ? 2,01 ? 0,34 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 10 ? 4,1 ?155,39?29,4225 ? 3,49 ? 1,97 ? 3,92 ? 0,00 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 10 ? 5,0 ?83,51 ?25,2375 ? 1,95 ? 1,44 ? 2,29 ? 0,47 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 11 ? 8,2 ?131,07?23,5775 ? 3,10 ? 1,73 ? 3,47 ? 0,00 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 11 ? 19,0 ?381,66?160,5300? 5,03 ? 2,95 ? 5,51 ? 1,72 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 11 ? 36,5 ?636,19?305,4775? 7,28 ? 3,80 ? 7,78 ? 2,58 ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 10 ? 67,4 ?529,71?92,9225 ? 6,23 ? 3,64 ? 7,05 ? 0,00 ?
???????????????????????????????????????????????????????????????????????????
5.8.4. Рисунок 4 представляет гистограммы расхождений между результатами измерений, расхождений между пробами и средних значений в элементах для уровня 6. Графики такого типа позволяют легко определить расхождения, возникающие от различных источников (между результатами измерений, пробами и лабораториями). Рисунок 4 показывает, что в этом эксперименте на уровне 6 имеется широкая вариация в средних значениях по элементам, так что, если метод испытаний будет соответствовать спецификации, то, вероятно, будут возникать разногласия между продавцом и покупателем из-за расхождений в результатах. Расхождения между пробами, которые меньше расхождений между результатами измерений (испытаний), означают, что разница между пробами на уровне 6 не является значительной.
%
? 1 A
?
25 ??
? 7
?
? 3 6 % %
?
20 ?? 10 ?? 10 ??
? 9 ? ? 4 B
? 4 ? 1 A ? 3 B
? ? ?
15 ?? 5 ?? 5 ?? A
? 8 ? 2 ? 4 6 9
? 2 5 B ? 7 8 ? 7
? ? 4 5 6 B ? 1 2 6 A
? ? 3 9 ? 2 3 5 5 7 8 8
10 ?? 0 ?? 0 ?? 1 9
Средние значения Расхождения Расхождения между
в элементах между пробами результатами измерений
(испытаний)
A, B - лаборатории N 10 и 11 соответственно
5.8.5. Значения статистик h и k, рассчитанные согласно 5.6.1, также представлены в таблицах 14 - 16 лишь для уровня 6. Для всех уровней эти значения изображены графически на рисунках 5 - 7; порядок уровней изменен так, чтобы общие средние по уровню располагались в порядке их возрастания, как показано в таблице 17. Рисунок 5 показывает, что лаборатория N 6 получила несколько высоких значений статистики k для расхождений между результатами измерений, что свидетельствует о ее худшей повторяемости по сравнению с остальными лабораториями. Рисунок 6 показывает, что три лаборатории (N 1, 6 и 10) получили высокие значения статистики k для расхождений между пробами, что могло произойти из-за отклонений от рекомендованной процедуры подготовки испытуемых навесок из проб. Рисунок 7 показывает устойчивые положительные или отрицательные значения статистики h в большинстве лабораторий (в лабораториях N 1, 6 и 10 вновь достигнуты наибольшие значения). Это прямое доказательство того, что в большинстве лабораторий имеется систематическая погрешность, свидетельствующая, что метод измерений (испытаний) неадекватно реализуется.
![]() между результатами измерений (сгруппированных
по лабораториям)
![]() между пробами (сгруппированных по лабораториям)
![]() средних значений в элементах (сгруппированных
по лабораториям)
5.8.6. Применение анализа данных по критериям Кохрена и Граббса, как описано в 5.6.2, дает результаты, представленные в таблице 18. Установлены два выброса. В отсутствие другой информации, данные, ответственные за это, должны быть исключены, а расчеты повторены. Анализ может быть затем продолжен в направлении исследования функциональных связей таким же способом, как в эксперименте по модели с однородными уровнями, рассмотренном в ГОСТ Р ИСО 5725-2.
Таблица 18
Пример 2. Значения статистик Кохрена и Граббса
??????????????????????????????????????????????????????????????????????????????????
? Уровень j ? Число ? Статистика Кохрена ? Статистика Кохрена для ?
? ?лабораторий ? для расхождений ? расхождений между пробами ?
? ? p' ? между результатами ? ?
? ? ? измерений ? ?
??????????????????????????????????????????????????????????????????????????????????
? 3 ? 11 ? 0,203 ? 0,664 <*> (1) ?
??????????????????????????????????????????????????????????????????????????????????
? 5 ? 11 ? 0,461 <**> (6) ? 0,374 ?
??????????????????????????????????????????????????????????????????????????????????
? 8 ? 10 ? 0,298 ? 0,465 ?
??????????????????????????????????????????????????????????????????????????????????
? 2 ? 10 ? 0,232 ? 0,238 ?
??????????????????????????????????????????????????????????????????????????????????
? 4 ? 11 ? 0,169 ? 0,550 ?
??????????????????????????????????????????????????????????????????????????????????
? 6 ? 11 ? 0,172 ? 0,301 ?
??????????????????????????????????????????????????????????????????????????????????
? 7 ? 11 ? 0,157 ? 0,536 ?
??????????????????????????????????????????????????????????????????????????????????
? 1 ? 10 ? 0,237 ? 0,680 <*> (6) ?
??????????????????????????????????????????????????????????????????????????????????
? Статистика Граббса для средних значений в элементах ?
??????????????????????????????????????????????????????????????????????????????????
? Уровень j ? Число ? Одно ? Два ?Два наибольших ? Одно ?
? ?лабораторий ?наименьшее ?наименьших? ? наибольшее ?
? ? p' ? ? ? ? ?
??????????????????????????????????????????????????????????????????????????????????
? 3 ? 11 ? 0,970 ? 0,791 ? 0,098 <**> ? 2,219 ?
? ? ? ? ? (1; 6) ? ?
??????????????????????????????????????????????????????????????????????????????????
? 5 ? 11 ? 1,396 ? 0,709 ? 0,302 ? 2,266 ?
??????????????????????????????????????????????????????????????????????????????????
? 8 ? 10 ? 0,849 ? - ? - ? 2,643 <**> ?
? ? ? ? ? ? (6) ?
??????????????????????????????????????????????????????????????????????????????????
? 2 ? 10 ? 1,259 ? 0,614 ? 0,466 ? 1,713 ?
??????????????????????????????????????????????????????????????????????????????????
? 4 ? 11 ? 1,290 ? 0,681 ? 0,294 ? 2,082 ?
??????????????????????????????????????????????????????????????????????????????????
? 6 ? 11 ? 1,108 ? 0,700 ? 0,479 ? 1,475 ?
??????????????????????????????????????????????????????????????????????????????????
? 7 ? 11 ? 1,649 ? 0,562 ? 0,453 ? 1,875 ?
??????????????????????????????????????????????????????????????????????????????????
? 1 ? 10 ? 1,808 ? 0,345 ? 0,590 ? 1,476 ?
??????????????????????????????????????????????????????????????????????????????????
? Примечание. Числа в скобках указывают лаборатории, которые обусловили ?
?квазивыбросы или выбросы. ?
? ?
? Критические значения следующие?
? ?
? критерий лабораторий в таблице ?
? ГОСТ Р ИСО ?
? 5725-2 ?
? ?
?Тест Кохрена Расхождения p' p ?
? между 10 20 0,389 0,480 ?
? результатами 11 22 0,365 0,450 ?
? измерений ?
? ?
?Тест Кохрена Расхождения 10 10 0,602 0,718 ?
? между 11 11 0,570 0,684 ?
? пробами ?
? ?
?Тест Граббса 10 10 2,290 2,482 ?
?- для Средние 11 11 2,355 2,564 ?
?одиночного значения ?
?выброса в элементах ?
? ?
?- для пары Средние 10 10 0,1864 0,1150 ?
?выбросов значения 11 11 0,2213 0,1448 ?
? в элементах ?
??????????????????????????????????????????????????????????????????????????????????
Для каждого уровня j вычисляют следующие статистики.
a) Общее среднее (суммирование по i, t и k)
, (39)где
b) "Вклады" лабораторий для каждой i (суммирование по t и k)
, (40)то есть равно среднему лаборатории минус общее среднее, где
c) "Вклады" проб для каждых i и t (суммирование по k)
, (41)то есть равно среднему пробы минус среднее лаборатории, где
d) Остатки для каждых i, t и k
, (42)то есть равно результату измерения минус среднее пробы.
e) Сумма квадратов для проб (суммирование по i)
. (43)f) Сумма квадратов для проб (суммирование по i и t)
. (44)g) Сумма квадратов для повторяемости (суммирование по i, t и k)
. (45)h) Степени свободы:
, , , (46)где
i) Факторы для каждого i (суммирование по t):
, (47) . (48)j) Факторы (суммирование по i):
, (49) , (50) . (51)k) Стандартные отклонения повторяемости
, (53) , (54)Примечание. Формулы (52 - 55) были получены с использованием статистической теории, разработанной Шеффе [7].
5.10.1. В качестве примера применения общих формул, необходимость которого возникает в связи с исключением некоторых результатов измерений, использованы данные примера 2 - уровень 4 (см. таблицу 19). Формулы, представленные в 5.9, дают общее среднее, указанное в таблице 19, а также суммы квадратов, степени свободы и факторы, приведенные в таблицах 20 - 22.
Таблица 19
Пример 3. Определение прочности с использованием
сульфата магния для уровня 4
???????????????????????????????????????????????????????????????????????????
? Номер ? Номер пробы t ? Результат измерения, % ?
? лаборатории i ? ??????????????????????????????????????
? ? ? k = 1 ? k = 2 ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 1 ? - ? 10,1 ?
? ? 2 ? 13,9 ? 13,8 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 1 ? - ? - ?
? ? 2 ? 8,3 ? 4,8 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 1 ? - ? 7,0 ?
? ? 2 ? - ? 12,0 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 1 ? 9,4 ? - ?
? ? 2 ? - ? - ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 1 ? 3,7 ? 6,3 ?
? ? 2 ? 3,3 ? 3,7 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 1 ? 16,5 ? 12,3 ?
? ? 2 ? 13,2 ? 16,8 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 1 ? 7,5 ? 9,3 ?
? ? 2 ? 11,1 ? 8,3 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 1 ? 5,7 ? 6,8 ?
? ? 2 ? 4,8 ? 5,5 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 1 ? 6,6 ? 7,0 ?
? ? 2 ? 4,9 ? 6,3 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 1 ? 5,5 ? 5,8 ?
? ? 2 ? 3,9 ? 5,7 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 1 ? 9,5 ? 7,2 ?
? ? 2 ? 8,1 ? 7,4 ?
???????????????????????????????????????????????????????????????????????????
? Общее среднее m = 8,1111%. ?
? j ?
? ?
? Число результатов измерений n = 36. ?
? j ?
???????????????????????????????????????????????????????????????????????????
Таблица 20
???????????????????????????????????????????????????????????????????????????
? Номер ? Среднее ? Число ? "Вклад" ? Фактор K ?
?лаборатории i?лаборатории, % ? результатов ? лабораторий ? ij ?
? ? ?измерений n ? B , % ? ?
? ? ? ij ? ij ? ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 12,600 ? 3 ? 4,4889 ? 5 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 6,550 ? 2 ? -1,5611 ? 4 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 9,500 ? 2 ? 1,3889 ? 2 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 9,400 ? 1 ? 1,2889 ? 1 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 4,250 ? 4 ? -3,8611 ? 8 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 14,700 ? 4 ? 6,5889 ? 8 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 9,050 ? 4 ? 0,9389 ? 8 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 5,700 ? 4 ? -2,4111 ? 8 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 6,200 ? 4 ? -1,9111 ? 8 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 5,225 ? 4 ? -2,8861 ? 8 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 8,050 ? 4 ? -0,0611 ? 8 ?
???????????????????????????????????????????????????????????????????????????
? 2 ?
? Сумма квадратов для лабораторий SS = 378,8531% . ?
? Lj ?
? ?
? Степени свободы для лабораторий ню = 11 - 1 = 10. ?
? Lj ?
? ?
? Факторы K = 130, К' = 68, K" = 19,6667. ?
? j j j ?
???????????????????????????????????????????????????????????????????????????
Таблица 21
Пример 3. Расчет суммы квадратов для проб
???????????????????????????????????????????????????????????????????????????
? Номер ? Номер пробы t ?Среднее пробы,? Число ?"Вклад" пробы?
?лаборатории i? ? % ? результатов ? H ,% ?
? ? ? ?измерений n ? ijt ?
? ? ? ? ijt? ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 1 ? 10,10 ? 1 ? -2,500 ?
? ? 2 ? 13,85 ? 2 ? 1,250 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 1 ? - ? 0 ? - ?
? ? 2 ? 6,55 ? 2 ? 0,000 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 1 ? 7,00 ? 1 ? -2,500 ?
? ? 2 ? 12,00 ? 1 ? 2,500 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 1 ? 9,40 ? 1 ? 0,000 ?
? ? 2 ? - ? 0 ? - ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 1 ? 5,00 ? 2 ? 0,750 ?
? ? 2 ? 3,50 ? 2 ? -0,750 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 1 ? 14,40 ? 2 ? -0,300 ?
? ? 2 ? 15,00 ? 2 ? 0,300 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 1 ? 8,40 ? 2 ? -0,650 ?
? ? 2 ? 9,70 ? 2 ? 0,650 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 1 ? 6,25 ? 2 ? 0,550 ?
? ? 2 ? 5,15 ? 2 ? -0,550 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 1 ? 6,80 ? 2 ? 0,600 ?
? ? 2 ? 5,60 ? 2 ? -0,600 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 1 ? 5,65 ? 2 ? 0,425 ?
? ? 2 ? 4,80 ? 2 ? -0,425 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 1 ? 8,35 ? 2 ? 0,300 ?
? ? 2 ? 7,75 ? 2 ? -0,300 ?
???????????????????????????????????????????????????????????????????????????
? 2 ?
? Сумма квадратов для проб SS = 29,9075% . ?
? Hj ?
? ?
? Степени свободы для проб ню = 20 - 11 = 9. ?
? Hj ?
???????????????????????????????????????????????????????????????????????????
Таблица 22
???????????????????????????????????????????????????????????????????????????
? Номер ? Номер пробы t ? Результат измерений, % ?
? лаборатории i ? ??????????????????????????????????????
? ? ? k = 1 ? k = 2 ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 1 ? - ? 0,00 ?
? ? 2 ? 0,05 ? -0,05 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 1 ? - ? - ?
? ? 2 ? 1,75 ? -1,75 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 1 ? - ? 0,00 ?
? ? 2 ? - ? 0,00 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 1 ? 0,00 ? - ?
? ? 2 ? - ? - ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 1 ? -1,30 ? 1,30 ?
? ? 2 ? -0,20 ? 0,20 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 1 ? 2,10 ? -2,10 ?
? ? 2 ? -1,80 ? 1,80 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 1 ? -0,90 ? 0,90 ?
? ? 2 ? 1,40 ? -1,40 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 1 ? -0,55 ? 0,55 ?
? ? 2 ? -0,35 ? 0,35 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 1 ? -0,20 ? 0,20 ?
? ? 2 ? -0,70 ? 0,70 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 1 ? -0,15 ? 0,15 ?
? ? 2 ? -0,90 ? 0,90 ?
???????????????????????????????????????????????????????????????????????????
? 11 ? 1 ? 1,15 ? -1,15 ?
? ? 2 ? 0,35 ? -0,35 ?
???????????????????????????????????????????????????????????????????????????
? 2 ?
? Сумма квадратов для повторяемости SS = 36,895% . ?
? rj ?
? ?
? Степени свободы для повторяемости ню = 36 - 20 = 16. ?
? rj ?
???????????????????????????????????????????????????????????????????????????
5.10.2. Применяя уравнения (52) - (55), получим:
,тогда
,и
,тогда
,и
,тогда
,и
.6.1.1. В ГОСТ Р ИСО 5725-2 данные, полученные в эксперименте по оценке прецизионности, рекомендуют проверять двумя тестами на наличие выбросов: тестами Кохрена и Граббса; при этом любые данные, которые увеличивают тестовую статистику в том или ином из этих тестов до значений, превышающих критические на уровне 1% значимости, должны быть отброшены (если у статистика нет обоснованного повода оставить эти данные). На практике применить эту процедуру часто нелегко. Рассмотрим результаты теста на выбросы в примере 1 в 4.8, представленные в таблице 8. Лаборатория N 5 дает только одно среднее значение в элементе (на уровне 10), достаточно экстремальное, чтобы по критерию Граббса квалифицировать его как выброс, но также дает три других квазивыброса, а данные на рисунке 3 прямо указывают, что в этой лаборатории что-то не в порядке. В этой ситуации специалист по статистике должен принять одно из решений:
a) сохранить все данные по лаборатории N 5;
b) отбросить только данные из уровня 10 по лаборатории N 5;
c) отбросить все данные лаборатории N 5.
Решение специалиста будет иметь существенное влияние на рассчитываемые значения стандартных отклонений повторяемости и воспроизводимости. В обычной практике анализа результатов экспериментов по оценке прецизионности данные, лежащие на линии, разделяющей квазивыбросы и выбросы, обнаруживаются достаточно часто, что может повлиять на результаты расчетов, что нежелательно. Робастные методы, описываемые в этом пункте, позволяют проанализировать полученные данные таким способом, при котором не требуется принимать решения, влияющие на результаты расчетов. Таким образом, если имеется основание ожидать, что результаты эксперимента по оценке прецизионности могут содержать выбросы, робастные методы могут быть предпочтительнее.
6.1.2. Основная модель, рассмотренная в разделе 5 ГОСТ Р ИСО 5725-1, содержит допущение по обоснованности установления общего значения для стандартного отклонения повторяемости для всех лабораторий, применяющих подтвержденный метод измерений. На практике часто оказывается, что некоторые лаборатории имеют худшую повторяемость, чем другие. Посмотрим, например, рисунок 5 для примера 2 в 5.8. Очевидно, что лаборатория N 6 имеет намного худшую повторяемость, чем лаборатория N 9 в этом эксперименте, так что допущение, что они достигли одинаковой повторяемости не кажется достоверным в этом случае. Некоторые участники эксперимента по оценке прецизионности могут получать плохую повторяемость, когда метод измерений подвергается такому эксперименту впервые или когда они имеют небольшой опыт в реализации этого метода измерений. В этих ситуациях использование робастных методов будет особенно предпочтительным.
6.1.3. Примером применения робастных методов [8] является случай, когда при анализе данных эксперимента по оценке прецизионности значения стандартных отклонений повторяемости и воспроизводимости рассчитывают таким образом, что на них не влияет наличие выбросов. Если всех участников эксперимента можно разделить на два класса: производящих данные хорошего и плохого качества, то робастные методы дадут значения стандартных отклонений повторяемости и воспроизводимости, которые действительны для класса с хорошим качеством данных, и не окажут воздействия на данные плохого качества (при условии, что класс данных плохого качества не слишком велик).
6.1.4. Использование робастных методов для анализа данных не влияет на планирование, организацию или выполнение эксперимента по оценке прецизионности. Решение об использовании робастных методов или методов выявления и удаления выбросов должно приниматься экспертом по статистике и представляться в совет экспертов. При использовании робастных методов в ходе обработки данных необходимо, как и в других случаях, проводить тесты на наличие выбросов, проверку совместимости (однородности), как это описано в ГОСТ Р ИСО 5725-2 или ГОСТ Р ИСО 5725-5, а также исследовать причины отдельных выбросов или графики по статистикам h и k. Однако сами исходные данные не должны исключаться как результаты этих измерений и проверок.
6.1.5. Знаменатели в формулах для статистик h и k являются стандартными отклонениями, которые в соответствии с методами расчета этих статистик, описанными в ГОСТ Р ИСО 5725-2, рассчитывают на основе представленных данных. Присутствие выбросов в этих данных будет изменять знаменатели, что приведет к искажениям в графиках этих статистик. Например, если на каком-то уровне эксперимента одна лаборатория выдает, что среднее значение в элементе является необычно большим выбросом, так что его величина намного больше, чем у любых других выбросов на том же уровне, то на графике статистики h это будет выглядеть в виде непомерно большого значения h для этого уровня. Однако значение статистики и для всех других лабораторий на этом же уровне будет малым, даже если несколько других лабораторий имеют выбросы. К подобному эффекту в расчетах статистики h может привести и использование общего среднего. В то же время использование робастных оценок стандартных отклонений как знаменателей в статистиках h и k и робастных оценок общих средних в расчете статистики h позволяет избежать этого искажения. Поэтому их и рекомендуется использовать для этих целей.
6.1.6. Данные эксперимента по оценке прецизионности позволяют рассчитать статистики двух типов:
a) средние значения в элементах, по которым рассчитывают стандартное отклонение, определяющее оценку межлабораторного расхождения;
b) стандартные отклонения или расхождения в пределах элементов (в том числе расхождения в эксперименте с распределенными уровнями), которые объединяют, чтобы получить оценку внутрилабораторного расхождения (вариации).
Робастные методы, описанные здесь, не подменяют эти средние значения в элементах, стандартные отклонения или расхождения (или вариации), различия, а обеспечивают альтернативные способы их сочетания для получения статистик, используемых для расчетов стандартных отклонений повторяемости и воспроизводимости.
Например, для значений одного уровня в эксперименте по модели с однородными уровнями, рассмотренном в ГОСТ Р ИСО 5725-2, первым этапом анализа является расчет среднего и стандартного отклонений результатов измерений в каждом элементе. Средние значения в элементах затем используют для расчетов стандартного отклонения, которое является оценкой межлабораторного расхождения. Когда используют робастные методы, изложенные в этом пункте, расчет выполняют с использованием Алгоритма A и средние значения в элементах не исключают из расчетов в результате применения к ним критерия Граббса. Также по этой модели эксперимента стандартные отклонения в элементах объединяют, чтобы оценить стандартное отклонение повторяемости. Если при этом использовать робастный анализ, то применяют Алгоритм S, который позволяет не исключать стандартные отклонения в элементах в результате использования критерия Кохрена. С любым подходом (описанным либо в ГОСТ Р ИСО 5725-2, либо здесь) обе эти оценки затем одинаковым образом используют для расчетов оценок стандартных отклонений повторяемости и воспроизводимости.
Более сложный пример шестифакторного ступенчато вложенного эксперимента приведен в приложении С ГОСТ Р ИСО 5725-3. Согласно этой модели первым этапом анализа является расчет средних значений по данным для каждой лаборатории (на каждом уровне), обозначаемых
6.1.7. Робастные методы, включенные в эту часть ГОСТ Р ИСО 5725, были выбраны потому, что они могут быть применимы ко всем экспериментальным моделям, приведенным в частях 2 - 5 ГОСТ Р ИСО 5725, а также потому, что предлагаемые в них расчеты относительно просты. Необходимо заметить, однако, что при этом обеспечиваются робастные способы объединения лишь средних значений, стандартных отклонений и расхождений в элементах. Описанные робастные методы не объединяют индивидуальные результаты измерений (испытаний), то есть они начинают с арифметических средних и стандартных отклонений в элементах. Имеются, однако, методы, которые объединяют результаты измерений (испытаний) в пределах элементов робастным способом, но они могут быть более сложными при применении на практике.
6.2.1. Этот алгоритм дает робастные величины среднего и стандартного отклонений данных, к которым он применяется, а именно:
a) средним значениям в элементах для любой модели;
b) расхождениям в элементах для модели с распределенными уровнями.
6.2.2. Обозначим индексом p общее число данных, расположенных в порядке возрастания:
Обозначим робастные среднее и стандартное отклонения этих данных x* и s*.
6.2.3. Рассчитаем первоначальные значения для x* и s* в виде:
Рассчитаем
. (58)Для каждого значения
рассчитывают: (59)Рассчитывают новые значения x* и s* по формулам:
6.2.5. Робастные оценки x* и s* могут быть получены итеративным расчетом, то есть повторением расчетов по 6.2.4 несколько раз, до тех пор, пока изменения в оценках x* и s* от одного расчета до следующего станут минимальными. Этот метод прост для программирования на компьютере.
6.2.6. Альтернативный метод без итерации легко применим для расчетов вручную с использованием уравнений (60), (61), которые можно представить в виде:
где
; ;x' и s' - средние значения и стандартные отклонения
значений элементов .Эти данные можно использовать, чтобы прямо рассчитать x* и s*, если известны
Еще одна возможность состоит в том, чтобы использовать итеративный метод для нахождения приближенного, а затем точного решения с помощью уравнений (62) и (63). Этот подход использован в примерах, приведенных ниже.
6.3.1. Этот алгоритм применяют для внутрилабораторного стандартного отклонения (или внутрилабораторных расхождений) в любой модели эксперимента. Он дает робастное среднеквадратичное значение для стандартных отклонений или расхождений, к которым применен.
6.3.2. Обозначим индексом p общее число данных, расположенных в порядке возрастания:
.(Это могут быть расхождения или стандартные отклонения).
Обозначим робастные среднеквадратичные значения w*, а число степеней свободы, связанных с каждым
). В таблице 23 находим соответствующие значения Таблица 23
???????????????????????????????????????????????????????????????????????????
? Степень свободы ? Ограничительный фактор ? Согласующий фактор ?
? ню ? эта ? кси ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 1,645 ? 1,097 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 1,517 ? 1,054 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 1,444 ? 1,039 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 1,395 ? 1,032 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 1,359 ? 1,027 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 1,332 ? 1,024 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 1,310 ? 1,021 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 1,292 ? 1,019 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 1,277 ? 1,018 ?
???????????????????????????????????????????????????????????????????????????
? 10 ? 1,264 ? 1,017 ?
???????????????????????????????????????????????????????????????????????????
? Примечание. Значения эта и кси выведены согласно Приложению B. ?
???????????????????????????????????????????????????????????????????????????
6.3.3. Найдем первоначальное значение для w* в виде
Рассчитывают
. (65)Для каждого
рассчитывают (66)Рассчитывают новое значение w* по формуле
6.3.5. Робастная оценка w* может быть получена итеративным способом повторением расчетов по 6.3.4 несколько раз, пока изменение оценки w* от первого расчета до последующего станет минимальным. Это простой метод для программирования на компьютере.
6.3.6. Альтернативный метод без использования итерации легко применим для расчетов вручную, аналогично описанному в 6.2.6. Уравнение (67) может быть представлено в виде
где
Это можно решить подбором, положив
Подход, который используют в примерах, приведенных ниже, состоит в использовании итеративного метода для приближенного решения, а затем в вычислении уравнения (68) для нахождения точного решения.
6.4. Формулы. Робастный анализ для отдельного уровня в эксперименте по модели с однородными уровнями
6.4.1. Робастная оценка стандартного отклонения повторяемости
Если в элементе имеются два результата измерений и Алгоритм S применяют к расхождениям в элементах, то
. (70)6.4.2. Робастная оценка стандартного отклонения средних значений в элементах
6.4.3. Затем может быть получено межлабораторное стандартное отклонение
где n - число результатов измерений в элементе.
Если выражение под корнем отрицательное, тогда принимают
Стандартное отклонение воспроизводимости для определенного уровня равно
6.5. Пример 4. Робастный анализ для отдельного уровня в эксперименте по модели с однородными уровнями
6.5.1. Пример 3 в ГОСТ Р ИСО 5725-2 иллюстрирует модель с однородными уровнями, в котором данные содержат квазивыбросы и выбросы. Уровень 5 в этом примере представляет определенный интерес, поскольку лаборатория N 1 дала среднее значение в элементах, близкое к квазивыбросу по критерию Граббса, а лаборатория N 6 - по критерию Кохрена. Эти данные представлены в таблице 24.
Таблица 24
(% креозота)
???????????????????????????????????????????????????????????????????????????
? Номер ? Данные ? Среднее ? Расхождение ?
?лаборатории i ? ? значение ? в элементе ?
? ? ? в элементе ? ?
???????????????????????????????????????????????????????????????????????????
? 1 ? 24,28 ? 24,00 ? 24,140 ? 0,28 ?
???????????????????????????????????????????????????????????????????????????
? 2 ? 20,40 ? 19,91 ? 20,155 ? 0,49 ?
???????????????????????????????????????????????????????????????????????????
? 3 ? 19,30 ? 19,70 ? 19,500 ? 0,40 ?
???????????????????????????????????????????????????????????????????????????
? 4 ? 20,30 ? 20,30 ? 20,300 ? 0,00 ?
???????????????????????????????????????????????????????????????????????????
? 5 ? 20,53 ? 20,88 ? 20,705 ? 0,35 ?
???????????????????????????????????????????????????????????????????????????
? 6 ? 18,56 ? 16,58 ? 17,570 ? 1,98 ?
???????????????????????????????????????????????????????????????????????????
? 7 ? 19,70 ? 20,50 ? 20,100 ? 0,80 ?
???????????????????????????????????????????????????????????????????????????
? 8 ? 21,10 ? 20,78 ? 20,940 ? 0,32 ?
???????????????????????????????????????????????????????????????????????????
? 9 ? 20,71 ? 21,66 ? 21,185 ? 0,95 ?
???????????????????????????????????????????????????????????????????????????
6.5.2. Если сохранить данные всех лабораторий, то стандартные отклонения повторяемости и воспроизводимости могут быть оценены с использованием формул по 7.4 ГОСТ Р ИСО 5725-2, которые дают:
p = 9;
m = 20,511;
; ; ; .6.5.3. Однако, согласно ГОСТ Р ИСО 5725-2, аналитик использовал информацию по другим уровням в эксперименте и сомневается в идентичности проб, испытанных лабораторией N 6, чтобы оправдать исключение обеих лабораторий N 1 и N 6 из расчетов, получая:
p = 7;
m = 20,412;
; ; ; .Ясно, что решение об исключении данных двух лабораторий оказало существенное влияние на оценки стандартных отклонений повторяемости и воспроизводимости.
6.5.4. Первым этапом в анализе является получение робастной оценки стандартного отклонения повторяемости. Расчеты могут быть представлены согласно таблице 25, в которой расхождения в элементах рассортированы в порядке возрастания. Применяя Алгоритм S, использующий итерацию, получим результаты, представленные в этой таблице. В этом примере число степеней свободы каждого расхождения в элементах составляет
и . Из четырех итераций, приведенных в таблице, ясно, что робастное значение w* = 0,7, и только одно расхождение в элементе превышает Таблица 25
(% креозота)
![]() ???????????????????????????????????????????????????????????????????????????
? Итерация ? 0 <1> ? 1 ? 2 ? 3 ? 4 ?
???????????????????????????????????????????????????????????????????????????
? пси ? - ? 0,66 ? 0,86 ? 1,00 ? 1,09 ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,00 ? 0,00 ? 0,00 ? 0,00 ? 0,00 ?
? 1 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,28 ? 0,28 ? 0,28 ? 0,28 ? 0,28 ?
? 2 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,32 ? 0,32 ? 0,32 ? 0,32 ? 0,32 ?
? 3 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,35 ? 0,35 ? 0,35 ? 0,35 ? 0,35 ?
? 4 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,40 ? 0,40 ? 0,40 ? 0,40 ? 0,40 ?
? 5 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,49 ? 0,49 ? 0,49 ? 0,49 ? 0,49 ?
? 6 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,80 ? 0,66 ? 0,80 ? 0,80 ? 0,80 ?
? 7 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,95 ? 0,66 ? 0,86 ? 0,95 ? 0,95 ?
? 8 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,98 ? 0,66 ? 0,86 ? 1,00 ? 1,09 ?
? 9 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Среднеквадратичные w? 0,83 ? 0,47 ? 0,56 ? 0,60 ? 0,62 ?
?Новые w* ?0,40 <2> ? 0,52 ? 0,61 ? 0,66 ? 0,68 ?
???????????????????????????????????????????????????????????????????????????
? <1> Значения получены из таблицы 24 после перестановки в порядке?
?возрастания. ?
? <2> Медианное (срединное) расхождение [см. уравнение (64)]. ?
???????????????????????????????????????????????????????????????????????????
Решение может быть также прямо получено следующим образом. Используя уравнение (68), в котором:
,получаем
,что дает решение (если предположение, что
w* = 0,69% креозота.
Можно затем подтвердить, что это значение дает
, как и предполагалось, только Следовательно, оценка стандартного отклонения повторяемости равна
креозота.6.5.5. Следующим этапом в анализе является получение робастной оценки стандартного отклонения средних значений в элементах. Применяя Алгоритм A к средним значениям, получим результаты, представленные в таблице 26, где средние значения в элементах рассортированы в порядке возрастания. Из четырех итераций, представленных в таблице, ясно, что робастными значениями являются x* = 20,412 и
отличаются от x* более чем на Таблица 26
в элементах (% креозота)
???????????????????????????????????????????????????????????????????????????
? Итерация ? 0 <1> ? 1 ? 2 ? 3 ? 4 ?
???????????????????????????????????????????????????????????????????????????
? фи ? - ? 1,424 ? 1,478 ? 1,514 ? 1,539 ?
???????????????????????????????????????????????????????????????????????????
? x* - фи ? - ? 18,876 ? 18,909 ? 18,893 ? 18,872 ?
? x* + фи ? - ? 21,724 ? 21,865 ? 21,921 ? 21,950 ?
???????????????????????????????????????????????????????????????????????????
? x* ? 17,570 ? 18,876 ? 18,909 ? 18,893 ? 18,872 ?
? 1 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 19,500 ? 19,500 ? 19,500 ? 19,500 ? 19,500 ?
? 2 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 20,100 ? 20,100 ? 20,100 ? 20,100 ? 20,100 ?
? 3 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 20,155 ? 20,155 ? 20,155 ? 20,155 ? 20,155 ?
? 4 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 20,300 ? 20,300 ? 20,300 ? 20,300 ? 20,300 ?
? 5 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 20,705 ? 20,705 ? 20,705 ? 20,705 ? 20,705 ?
? 6 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 20,940 ? 20,940 ? 20,940 ? 20,940 ? 20,940 ?
? 7 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 21,185 ? 21,185 ? 21,185 ? 21,185 ? 21,185 ?
? 8 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 24,140 ? 21,724 ? 21,865 ? 21,921 ? 21,950 ?
? 9 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Среднее ? 20,511 ? 20,387 ? 20,407 ? 20,411 ? 20,412 ?
???????????????????????????????????????????????????????????????????????????
?Стандартное ? 1,727 ? 0,869 ? 0,890 ? 0,905 ? 0,916 ?
?отклонение ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Новые x* ?20,300 <2> ? 20,387 ? 20,407 ? 20,411 ? 20,412 ?
?Новые s* ? 0,949 <2> ? 0,985 ? 1,009 ? 1,026 ? 1,039 ?
???????????????????????????????????????????????????????????????????????????
? <1> Значения получены из таблицы 24 после перестановки в порядке?
?возрастания. ?
???????????????????????????????????????????????????????????????????????????
При расчете вручную аналитик должен использовать прямой метод, описанный в 6.2.6, например
.Это дает x' = 20,412 и s' = 0,573% креозота.
Отсюда из уравнений (62) и (63)
![]() получаем s* = 1,070% креозота и x* = x' = 20,412% креозота.
Можно затем подтвердить, что значение s* дает
(тогда как предполагалось, что только Оценку межлабораторного стандартного отклонения проводят по формуле (72):
,а оценку стандартного отклонения воспроизводимости - по формуле (74):
.6.6. Формулы. Робастный анализ для отдельного уровня в эксперименте по модели с разделенными уровнями
6.6.1. Применительно к модели с разделенными уровнями робастная оценка стандартного отклонения повторяемости
6.6.2. Робастная оценка стандартного отклонения средних значений
Для оценки стандартного отклонения воспроизводимости на определенном уровне модели можно использовать формулы, приведенные в 4.5.6.
6.7. Пример 5. Робастный анализ для отдельного уровня в эксперименте по модели с разделенными уровнями
6.7.1. Данные примера 1 в 4.8 содержали несколько квазивыбросов и один выброс (см. таблицу 8). Кроме того, на рисунке 3 видна отрицательная систематическая погрешность в результатах лаборатории N 5. Если аналитик не может выявить источники этих аномалий, он попадает в трудное положение при принятии решения, какие данные следует исключить из расчетов стандартных отклонений повторяемости и воспроизводимости. Для иллюстрации результатов робастного анализа здесь использованы данные уровня 14 (см. таблицу 4).
6.7.2. Для получения робастной оценки стандартного отклонения повторяемости расхождений в элементах применяют Алгоритм A (см. таблицу 5), что приводит к результатам, показанным в таблице 27, в которой расхождения в элементах рассортированы в порядке возрастания. Из четырех итераций, представленных в таблице, видно, что робастные значения равны x* = 8,29, s* = 0,36 и что только
Таблица 27
Пример 5. Применение Алгоритма A к расхождениям в элементах
(% протеина)
???????????????????????????????????????????????????????????????????????????
? Итерация ? 0 ? 1 ? 2 ? 3 ? 4 ?
???????????????????????????????????????????????????????????????????????????
? фи ? - ? 0,53 ? 0,56 ? 0,55 ? 0,54 ?
???????????????????????????????????????????????????????????????????????????
? x* - фи ? - ? 7,85 ? 7,74 ? 7,74 ? 7,75 ?
? x* + фи ? - ? 8,91 ? 8,86 ? 8,84 ? 8,83 ?
???????????????????????????????????????????????????????????????????????????
? x* ? 7,81 ? 7,85 ? 7,81 ? 7,81 ? 7,81 ?
? 1 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 7,93 ? 7,93 ? 7,93 ? 7,93 ? 7,93 ?
? 2 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 8,13 ? 8,13 ? 8,13 ? 8,13 ? 8,13 ?
? 3 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 8,14 ? 8,14 ? 8,14 ? 8,14 ? 8,14 ?
? 4 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 8,38 ? 8,38 ? 8,38 ? 8,38 ? 8,38 ?
? 5 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 8,40 ? 8,40 ? 8,40 ? 8,40 ? 8,40 ?
? 6 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 8,44 ? 8,44 ? 8,44 ? 8,44 ? 8,44 ?
? 7 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 8,52 ? 8,52 ? 8,52 ? 8,52 ? 8,52 ?
? 8 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 9,31 ? 8,91 ? 8,86 ? 8,84 ? 8,83 ?
? 9 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Среднее ? 8,340 ? 8,300 ? 8,290 ? 8,288 ? 8,287 ?
???????????????????????????????????????????????????????????????????????????
?Стандартное ? 0,436 ? 0,326 ? 0,322 ? 0,317 ? 0,315 ?
?отклонение ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Новые x* ?8,380 <1>? 8,300 ? 8,290 ? 8,288 ? 8,287 ?
?Новые s* ?0,356 <1>? 0,370 ? 0,365 ? 0,359 ? 0,357 ?
???????????????????????????????????????????????????????????????????????????
? <1> Получено по формулам (56) и (57). ?
???????????????????????????????????????????????????????????????????????????
x' = 8,219 и s' = 0,257% протеина.
Тогда уравнения (62) и (63) в 6.2.6 можно записать в виде
![]() и
,что дает s* = 0,354% протеина, а, используя уравнение (75), получим
протеина.Робастное среднее значение для расхождений в элементах составляет
.При этих значениях x* и s*
.Тогда
и протеина.Это подтверждает, что в расчетах x* и s* только
6.7.3. Применение Алгоритма A к средним значениям в элементах (из таблицы 6) дает результаты, представленные в таблице 28, в которой средние значения в элементах расположены в порядке возрастания. Ситуация подобна представленной в таблице 26, а именно
, получим, что среднее значение и стандартное отклонение от x' = 85,486 и s' = 0,209.
Таблица 28
в элементах (% протеина)
???????????????????????????????????????????????????????????????????????????
? Итерация ? 0 ? 1 ? 2 ? 3 ? 4 ?
???????????????????????????????????????????????????????????????????????????
? фи ? - ? 0,446 ? 0,492 ? 0,519 ? 0,537 ?
???????????????????????????????????????????????????????????????????????????
? x* - фи ? - ? 85,104 ? 85,009 ? 84,971 ? 84,950 ?
? x* + фи ? - ? 85,996 ? 85,993 ? 86,009 ? 86,024 ?
???????????????????????????????????????????????????????????????????????????
? x* ? 84,525 ? 85,104 ? 85,009 ? 84,971 ? 84,950 ?
? 1 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 85,140 ? 85,140 ? 85,140 ? 85,140 ? 85,140 ?
? 2 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 85,345 ? 85,345 ? 85,345 ? 85,345 ? 85,345 ?
? 3 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 85,385 ? 85,385 ? 85,385 ? 85,385 ? 85,385 ?
? 4 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 85,550 ? 85,550 ? 85,550 ? 85,550 ? 85,550 ?
? 5 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 85,575 ? 85,575 ? 85,575 ? 85,575 ? 85,575 ?
? 6 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 85,660 ? 85,660 ? 85,660 ? 85,660 ? 85,660 ?
? 7 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 85,750 ? 85,750 ? 85,750 ? 85,750 ? 85,750 ?
? 8 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 86,170 ? 85,996 ? 85,993 ? 86,009 ? 86,024 ?
? 9 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Среднее ? 85,456 ? 85,501 ? 85,490 ? 85,487 ? 85,487 ?
???????????????????????????????????????????????????????????????????????????
?Стандартное ? 0,453 ? 0,289 ? 0,305 ? 0,316 ? 0,324 ?
?отклонение ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Новые x* ?85,550 <1>? 85,501 ? 85,490 ? 85,487 ? 85,487 ?
?Новые s* ?0,297 <1> ? 0,328 ? 0,346 ? 0,358 ? 0,367 ?
???????????????????????????????????????????????????????????????????????????
? <1> Получено по формулам (56) и (57). ?
???????????????????????????????????????????????????????????????????????????
Значит, на основе уравнения (63) может быть получено s* из равенства
,откуда s* = 0,390% протеина.
Теперь можно вычислить x* по формуле (62) в 6.2.6, что дает x* = 85,486% протеина.
Для контроля правильности решения рассчитывают
, , .Очевидно, что, как и предполагалось, только
Для получения стандартного отклонения воспроизводимости используют уравнение (76) в 6.6.2, которое дает
протеина, а затем уравнение (13) в 4.5.6, чтобы получить протеина.Из этого примера следует, что робастный метод дает несколько меньшие оценки
6.8.1. В модели для гетерогенного материала в обычном случае, когда две пробы подготовлены для каждой из p' лабораторий на каждом уровне и при этом получают по два результата испытаний на каждой пробе, робастные оценки стандартных отклонений повторяемости и воспроизводимости могут быть получены троекратным применением Алгоритмов A и S, как описано ниже.
a) Применяют Алгоритм S к расхождениям между результатами испытаний, чтобы получить робастное значение w* из уравнения (67), и определяют
b) Применяют Алгоритм S к расхождениям между пробами, чтобы получить другое робастное значение w* из уравнения (67), и определяют
c) Применяют Алгоритм A к средним значениям в элементах, чтобы получить робастное значение s* из уравнения (61), и находят
Эти расчеты могут быть удобно представлены в табличной форме с размещением в первой графе значений расхождений или средних значений в порядке возрастания, как это показано на примерах, описанных ниже.
6.8.2. Затем для расчетов оценок стандартных отклонений повторяемости и воспроизводимости и стандартного отклонения
6.9.1. Данные для уровня 6 в примере 2 (см. 5.8) не содержат выбросов или квазивыбросов, и они были использованы для иллюстрации результатов, которые получают робастным анализом в таком случае.
6.9.2. Применение Алгоритма S к расхождениям между результатами измерений (из таблицы 14) иллюстрируется данными, представленными в таблице 29. Здесь число степеней свободы
, и число индексированных значений p = 2p' = 22. Из четырех итераций, представленных в таблице, видно, что робастным значением является и что значения от ,так что уравнение (68) принимает вид
.Отсюда получаем w* = 4,30%.
Таблица 29
между результатами измерений
![]() ???????????????????????????????????????????????????????????????????????????
? Итерация ? 0 ? 1 ? 2 ? 3 ? 4 ?
???????????????????????????????????????????????????????????????????????????
? пси ? - ? 3,9 ? 5,1 ? 5,9 ? 6,4 ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,1 ? 0,1 ? 0,1 ? 0,1 ? 0,1 ?
? 1 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 0,6 ? 0,6 ? 0,6 ? 0,6 ? 0,6 ?
? 2 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,1 ? 1,1 ? 1,1 ? 1,1 ? 1,1 ?
? 3 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,1 ? 1,1 ? 1,1 ? 1,1 ? 1,1 ?
? 4 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,2 ? 1,2 ? 1,2 ? 1,2 ? 1,2 ?
? 5 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,3 ? 1,3 ? 1,3 ? 1,3 ? 1,3 ?
? 6 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,4 ? 1,4 ? 1,4 ? 1,4 ? 1,4 ?
? 7 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,6 ? 1,6 ? 1,6 ? 1,6 ? 1,6 ?
? 8 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,8 ? 1,8 ? 1,8 ? 1,8 ? 1,8 ?
? 9 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 2,1 ? 2,1 ? 2,1 ? 2,1 ? 2,1 ?
? 10 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 2,2 ? 2,2 ? 2,2 ? 2,2 ? 2,2 ?
? 11 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 2,5 ? 2,5 ? 2,5 ? 2,5 ? 2,5 ?
? 12 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 2,6 ? 2,6 ? 2,6 ? 2,6 ? 2,6 ?
? 13 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 3,9 ? 3,9 ? 3,9 ? 3,9 ? 3,9 ?
? 14 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 4,0 ? 3,9 ? 4,0 ? 4,0 ? 4,0 ?
? 15 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 4,4 ? 3,9 ? 4,4 ? 4,4 ? 4,4 ?
? 16 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 4,6 ? 3,9 ? 4,6 ? 4,6 ? 4,6 ?
? 17 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 5,5 ? 3,9 ? 5,1 ? 5,5 ? 5,5 ?
? 18 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 7,4 ? 3,9 ? 5,1 ? 5,9 ? 6,4 ?
? 19 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 7,6 ? 3,9 ? 5,1 ? 5,9 ? 6,4 ?
? 20 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 8,1 ? 3,9 ? 5,1 ? 5,9 ? 6,4 ?
? 21 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 8,1 ? 3,9 ? 5,1 ? 5,9 ? 6,4 ?
? 22 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Новые w ? 4,17 ? 2,80 ? 3,29 ? 3,55 ? 3,70 ?
?Новые w* ?2,35 <1> ? 3,07 ? 3,61 ? 3,89 ? 4,06 ?
???????????????????????????????????????????????????????????????????????????
? <1> Получено по формуле (64). ?
???????????????????????????????????????????????????????????????????????????
Подтверждением правильности решения по определению w* является то, что если w* = 4,3%, то
и четыре значения от Используя уравнение (77), получим
.6.9.3. Применяя второй раз Алгоритм S к расхождениям между пробами (из таблицы 15), получаем результаты, приведенные в таблице 30. Из четырех итераций, представленных в этой таблице, видно, что робастным значением является w* = 40 и что
,так что уравнение (68) принимает вид
.Отсюда получаем w* = 4,23%.
Таблица 30
между пробами
![]() ???????????????????????????????????????????????????????????????????????????
? Итерация ? 0 ? 1 ? 2 ? 3 ? 4 ?
???????????????????????????????????????????????????????????????????????????
? пси ? - ? 4,19 ? 5,43 ? 6,10 ? 6,45 ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,00 ? 1,00 ? 1,00 ? 1,00 ? 1,00 ?
? 1 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 1,70 ? 1,70 ? 1,70 ? 1,70 ? 1,70 ?
? 2 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 2,05 ? 2,05 ? 2,05 ? 2,05 ? 2,05 ?
? 3 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 2,25 ? 2,25 ? 2,25 ? 2,25 ? 2,25 ?
? 4 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 2,55 ? 2,55 ? 2,55 ? 2,55 ? 2,55 ?
? 5 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 2,55 ? 2,55 ? 2,55 ? 2,55 ? 2,55 ?
? 6 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 3,15 ? 3,15 ? 3,15 ? 3,15 ? 3,15 ?
? 7 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 3,35 ? 3,35 ? 3,35 ? 3,35 ? 3,35 ?
? 8 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 4,40 ? 4,19 ? 4,40 ? 4,40 ? 4,40 ?
? 9 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 6,75 ? 4,19 ? 5,43 ? 6,10 ? 6,45 ?
? 10 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? w* ? 6,95 ? 4,19 ? 5,43 ? 6,10 ? 6,45 ?
? 11 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Новые w ? 3,82 ? 3,01 ? 3,38 ? 3,58 ? 3,69 ?
?Новые w* ?2,55 <1> ? 3,30 ? 3,71 ? 3,92 ? 4,05 ?
???????????????????????????????????????????????????????????????????????????
? <1> Получено по формуле (64). ?
???????????????????????????????????????????????????????????????????????????
К сожалению, этому соответствует значение
и это - недействительное решение, так как Приняв
.Тогда уравнение (68) примет вид
.Отсюда получаем w* = 4,18%.
Теперь
, и можно убедиться в действительности решения, так как только Используя уравнение (78) в 6.8.1, получим
.6.9.4. Применяя Алгоритм A к средним значениям в элементах (из таблицы 16), получим результаты, представленные в таблице 31. Расчеты сошлись после двух итераций и дали s* = 5,70 (при этом ни одно
Таблица 31
Пример 6. Применение Алгоритма A к средним значениям
в элементах (%)
???????????????????????????????????????????????????????????????????????????
? Итерация ? 0 ? 1 ? 2 ? 3 ? 4 ?
???????????????????????????????????????????????????????????????????????????
? фи ? - ? 10,005 ? 8,550 ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* - фи ? - ? 8,245 ? 10,450 ? ? ?
? x* + фи ? - ? 28,255 ? 27,550 ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 13,425 ? 13,425 ? 13,425 ? ? ?
? 1 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 13,425 ? 13,425 ? 13,425 ? ? ?
? 2 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 13,750 ? 13,750 ? 13,750 ? ? ?
? 3 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 14,475 ? 14,475 ? 14,475 ? ? ?
? 4 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 17,075 ? 17,075 ? 17,075 ? ? ?
? 5 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 18,250 ? 18,250 ? 18,250 ? ? ?
? 6 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 21,000 ? 21,000 ? 21,000 ? ? ?
? 7 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 21,225 ? 21,225 ? 21,225 ? ? ?
? 8 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 23,675 ? 23,675 ? 23,675 ? ? ?
? 9 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 26,275 ? 26,275 ? 26,275 ? ? ?
? 10 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
? x* ? 26,425 ? 26,425 ? 26,425 ? ? ?
? 11 ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Среднее ? 19,00 ? 19,00 ? 19,00 ? ? ?
???????????????????????????????????????????????????????????????????????????
?Стандартное ? 5,03 ? 5,03 ? 5,03 ? ? ?
?отклонение ? ? ? ? ? ?
???????????????????????????????????????????????????????????????????????????
?Новые x* ?18,25 <1>? 19,00 ? 19,00 ? ? ?
?Новые s* ?6,67 <1> ? 5,70 ? 5,70 ? ? ?
???????????????????????????????????????????????????????????????????????????
? <1> Получено по формулам (56) и (57). ?
???????????????????????????????????????????????????????????????????????????
Используя уравнение (79), получим
.6.9.5. С использованием результатов, полученных в 6.9.2 - 6.9.4, уравнения (29) - (33) в 5.5.5 теперь дают:
, , .Тогда
, , .Следовательно, в этом примере робастный метод дает оценки
ИСПОЛЬЗУЕМЫЕ В ГОСТ Р ИСО 5725
a - отсекаемый на оси ординат отрезок в соотношении s = a + bm
A - показатель, используемый для расчета неопределенности оценки
b - угловой коэффициент прямой в соотношении s = a + bm
B - лабораторная составляющая систематической погрешности измерений при реализации конкретного метода - разность между систематической погрешностью лаборатории при реализации конкретного метода измерений (конкретной МВИ) и систематической погрешностью метода измерений
и т.д. - составляющие величины B, представляющие факторы, которые изменяются в условиях промежуточной прецизионностиc - отсекаемый на оси ординат отрезок в соотношении
![]() C, C', C" - тестовые статистики
- критические значения для статистикd - угловой коэффициент прямой в соотношении
![]() e - составляющая результата измерений, представляющая случайную погрешность каждого результата измерений
f - коэффициент критического диапазона
- p-квантиль F-распределения с G - статистика Граббса
h - статистика Манделя для межлабораторной совместимости
k - статистика Манделя для внутрилабораторной совместимости
LCL - нижний предел контроля (действия либо предупреждения)
m - общее среднее значение измеряемой характеристики; уровень
M - количество факторов, рассматриваемых в условиях промежуточной прецизионности
N - количество повторений (итераций)
n - количество результатов измерений, полученных в одной лаборатории на одном уровне (т.е. в пределах ячейки - базового элемента)
p - количество лабораторий, участвующих в межлабораторном эксперименте
P - вероятность
q - количество уровней измеряемой характеристики в межлабораторном эксперименте
r - предел повторяемости (сходимости)
R - предел воспроизводимости
RM - стандартный образец
s - оценка стандартного (среднеквадратического) отклонения
T - итог или сумма какого-либо выражения
t - количество объектов испытаний или групп объектов
UCL - верхний предел контроля (действия либо предупреждения)
W - весовой коэффициент, используемый при расчете взвешенной регрессии
w - диапазон изменения выборки результатов измерений
x - заданная величина, используемая для критерия Граббса
y - результат измерений (или результат испытаний)
![]() Символы, используемые в качестве подстрочных индексов
C - различие, определяемое калибровкой
E - различие, определяемое оборудованием
i - идентификатор для конкретной лаборатории
I() - идентификатор для промежуточных мер прецизионности; в скобках - идентификация типа промежуточной ситуации
j - идентификатор для уровня (ГОСТ Р ИСО 5725-2)
Идентификатор для группы испытаний или для фактора (ГОСТ Р ИСО 5725-3)
k - идентификатор для конкретного результата испытаний в лаборатории i на уровне j
L - межлабораторный
m - идентификатор для поддающейся обнаружению систематической погрешности
M - различие, обусловленное неидентичностью проб (образцов)
O - различие, определяемое сменой оператора
P - вероятность
r - повторяемость
R - воспроизводимость
T - различие, обусловленное периодом (временем), в течение которого проводят измерения или оценочный эксперимент
W - внутрилабораторный
1, 2, 3... для результатов измерений, нумеруемых в порядке их получения
(1), (2), (3)... для результатов измерений (или результатов испытаний), нумеруемых в порядке возрастания измеряемой величины
Дополнительные условные обозначения и сокращения,
использованные в ГОСТ Р ИСО 5725-5
D - расхождение в пределах базового элемента в эксперименте с разделенными уровнями
g - число проб (образцов), испытанных в лаборатории на одном уровне
H - составляющая погрешности результата измерения, представляющая случайную погрешность, вызванную неоднородностью (различиями) проб (пробы)
K - функция числа результатов измерений в элементах
p' - число лабораторий, участвующих в межлабораторном эксперименте
SS - сумма квадратов
z - остаток
Дополнительные условные обозначения, используемые
в качестве подстрочных индексов в ГОСТ Р ИСО 5725-5
a, b - идентификаторы для проб в эксперименте по модели с разделенными уровнями
t - идентификатор для пробы в лаборатории i на уровне j
H - между пробами
Дополнительное условное обозначение,
используемое в качестве надстрочного индекса
в ГОСТ Р ИСО 5725-5
* - робастная оценка
(справочное)
B.1. Введение
Использование робастных методов анализа данных экспериментов по оценке прецизионности предложено Комитетом аналитических методов Королевского Химического Общества Соединенного Королевства [8]. Алгоритм A в настоящем стандарте взят из публикации так же, как и коэффициент 1,134, использованный для расчета s* в Алгоритме A (это значение
Алгоритм S схож с процедурой, приведенной в [8] для специального случая, в котором каждая лаборатория представляет n = 2 измерений для каждого уровня. Это делает удобным применение метода робастного анализа в экспериментах по оценке прецизионности с более чем двумя факторами (такими, как эксперимент для гетерогенного материала в разделе 5 настоящего стандарта и модели в ГОСТ Р ИСО 5725-3). Вывод факторов, используемых в алгоритме S, приведен ниже.
B.2. Условные обозначения, применяемые в настоящем приложении
s - стандартное отклонение, дающее оценку
.B.3. Вывод ограничительного фактора
Согласующий фактор
. (B.1)Это требование может быть записано в виде
где случайная величина в фигурных скобках, тесно связанная с
, имеет распределение Плотность вероятности распределения
, (B.3)так что
потому что предел
.Второй член в правой части (B.4) равен
Для Алгоритма S ограничительный фактор
Биометрические таблицы для распределения
Первый член в правой части уравнения (B.4) можно представить в виде
.При
хорошо известное свойство гамма-функции имеет вид .Тогда этот первый член можно переписать в виде
Следовательно, для данного числа степеней свободы
![]() или
. (B.8)Это равенство может быть использовано для получения значений согласующего фактора
(справочное)
ВЫВОД РАВЕНСТВ, ИСПОЛЬЗУЕМЫХ ДЛЯ РОБАСТНОГО АНАЛИЗА
Равенства (62) и (63), используемые для расчета робастных величин среднего значения и стандартного отклонения методом, описанным в 6.2.6, могут быть получены из соотношений (60) и (61) Алгоритма A следующим образом.
(C.2)и
где
значений по пунктам .Значит, уравнение (C.1) может быть записано в виде
.Тогда
![]() или
, (C.4)что является равенством (62).
Для получения уравнения (63) из уравнения (61) заметим, что сумма в уравнении (61) может быть представлена следующим образом
. (C.5)Подставляя x* в сумму в правой части выражения для x*, после некоторых алгебраических преобразований получим
. (C.6)Используя определение s' в уравнении (C.3), можно записать
. (C.7)Подставив уравнение (C.7) в уравнение (61), получим уравнение (63).
(справочное)
БИБЛИОГРАФИЯ
[1] ISO 3534-1:1993. Statistics - Vocabulary and symbols - Part 1: Statistical methods. Terms and definitions
[3] Mandel, J. and Lashof, T.W. Interpretation and Generalization of Youden's Two-Sample Diagram Journal of Quality Technology
[5] BS 812-103:1985, Testing aggregates - Part 103: Methods for determination of particle size distribution British Standards Institution
[6] BS 812-121:1989, Testing aggregates - Part 121: Methods for determination of soundness British Standards Institution
[8] Analytical Methods Committee Robust statistics - How not to reject outliers - Part 1: Basic concepts - Part 2: Inter-laboratory trials The Analyst
[9] SWEENEY, An inter-laboratory study of the determination of protein by combustion in feeds Journal of the Association of Official Analytical Chemists.
Вернуться в "Каталог нормативных документов"
Источник информации: https://internet-law.ru/documents/prod/gost-r_gosudarstvennyj-standart/46/gost_45875.html
На правах рекламы:
|