Пробы а и b должны быть распределены среди участников случайным образом, причем процедуры рандомизации для а и b должны быть независимы. При этом необходимо, чтобы эксперты-статистики имели точную информацию о том, какие результаты были получены на материале а и какие - на материале b на каждом уровне эксперимента. Однако пробы следует зашифровать так, чтобы скрыть эту информацию от участников эксперимента.
Таблица 1 - Рекомендуемая форма для сравнения данных эксперимента с разделенными уровнями
Номер лаборатории |
Уровень |
|||||||||
|
1 |
2 |
|
j |
|
q |
||||
|
а |
b |
а |
b |
|
а |
b |
|
а |
b |
1 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
|
|
|
4.4 Статистическая модель
4.4.1 Основная модель, используемая в настоящем стандарте, дана равенством (1) в разделе 5 ГОСТ Р ИСО 5725-1. Там установлено, что для оценивания точности (правильности и прецизионности) метода измерений каждый результат измерения полезно представлять как сумму трех составляющих:
yijk = mj + Bij + eijk,(3)
где для определенного испытуемого материала:
mj - общее среднее значение для определенного уровня j = 1, ..., q;
Bij - лабораторная составляющая систематической погрешности в условиях повторяемости в определенной лаборатории i = 1, ..., р на определенном уровне j = 1, ..., q;
eijk - случайная погрешность результата измерений k = 1, ..., п, полученная в лаборатории i на уровне j в условиях повторяемости.
4.4.2 Для эксперимента с разделенными уровнями эта модель принимает вид
yijk = mjk + Bij + eijk.(4)
Это неравенство отличается от равенства (3) только одной деталью: индекс k в mjk означает, что в соответствии с равенством (4) общее среднее значение может теперь зависеть от материала а или b (k = 1 или 2) на уровне j.
Отсутствие индекса k в Вij означает допущение, что систематическая ошибка, связанная с лабораторией i, не зависит от материала а или b на определенном уровне. Вот почему так важно, чтобы эти два материала были бы однородными (одинаковыми).
4.4.3 Определяют среднее значение в базовом элементе (ячейке)
yij = (yija + yijb)/2(5)
и внутриэлементное расхождение (разброс)
Dij = yija - yijb.(6)
4.4.4 Общее среднее значение для уровня j в эксперименте с разделенными уровнями может быть определено как
mj = (mja + mjb)/2.(7)
4.5 Статистический анализ данных эксперимента с разделенными уровнями
4.5.1 Данные эксперимента сводят в таблицу (см. таблицу 1). Каждая комбинация лаборатории и уровня дает базовый элемент (ячейку) в этой таблице, а также содержит два результата yija и yijb.
Рассчитывают Dij - расхождения в элементах и сводят их в таблицу (см. таблицу 2). Метод анализа требует, чтобы все расхождения были рассчитаны с сохранением знака разности
а - b.
Рассчитывают средние значения уij и сводят их в таблицу (см. таблицу 3).
4.5.2 Если элемент в таблице 1 не содержит двух результатов измерений (например потому, что пробы были испорчены или данные исключены в последующем как выбросы), то соответствующие элементы в таблицах 2 и 3 оставляют пустыми.
4.5.3 Для каждого уровня j эксперимента рассчитывают среднее Dj и стандартное sDj отклонения расхождений в графе j таблицы 2 по формулам:
,(8)
,(9)
где Σ - знак суммирования по всем лабораториям i = 1, 2, ..., р.
Если в таблице 2 имеются пустые элементы, то р теперь становится числом элементов в графе j таблицы 2, содержащих данные, и суммирование выполняют без пустых элементов.
4.5.4 Для каждого уровня j в эксперименте рассчитывают среднее и стандартное отклонения средних значений в графе j таблицы 3, используя формулы:
,(10)
,(11)
где Σ - знак суммирования по всем лабораториям i = 1, 2, ..., р.
Если в таблице 3 имеются пустые элементы, то р теперь становится числом элементов в графе j, содержащих данные, и суммирование выполняют без пустых элементов.
4.5.5 Для проверки совместимости данных и наличия выбросов, как описано в 4.6, используют таблицы 2, 3 и статистики, рассчитанные по формулам (8-11). При исключении данных пересчитывают статистики.
4.5.6 Рассчитывают стандартные отклонения повторяемости srj и воспроизводимости sRj по формулам:
,(12)
.(13)
4.5.7 Исследуют, зависят ли srj и sRj от среднего уj и, если это так, находят соответствующие функциональные соотношения, используя методы, описанные в 7.5 ГОСТ Р ИСО 5725-2.
Таблица 2 - Рекомендуемая форма табулирования расхождений в базовых элементах для эксперимента с разделенными уровнями
Номер лаборатории |
Уровень |
|||||
|
1 |
2 |
|
j |
|
q |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
Таблица 3 - Рекомендуемая форма табулирования средних значений в базовых элементах для эксперимента с разделенными уровнями
Номер лаборатории |
Уровень |
|||||
|
1 |
2 |
|
j |
|
q |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
4.6 Исследование данных на совместимость и наличие выбросов
4.6.1 Проверяют данные на совместимость, используя статистику h, описанную в 7.3.1 ГОСТ Р ИСО 5725-2.
Чтобы проконтролировать совместимость расхождений в базовых элементах, рассчитывают серию для статистики h по формуле
hij = (Dij - Dj)/sDj.(14)
Для контроля совместимости средних значений в базовых элементах рассчитывают серию для статистики h по формуле
hij = (yij - yj)/syj.(15)
Для оценки различий лабораторий с точки зрения совместимости полученных данных, наносят на график обе серии в порядке возрастания уровней, но сгруппировав их по лабораториям, как показано на рисунках 2 и 3. Интерпретация этих графиков подробно рассмотрена в 7.3.1 ГОСТ Р ИСО 5725-2. Если лаборатория получила худшую повторяемость по сравнению с другими, это будет видно по необычно большому числу больших значений h на графике, построенном по расхождениям в элементах. Если данные лаборатории, в основном, содержат систематическую погрешность, то это будет видно по значениям h на графике, построенном для средних значений в элементах: большинство из них расположится в одном направлении. В любом случае лаборатория должна изучить причины расхождений и доложить о них организатору эксперимента.
4.6.2 Для контроля данных на наличие квазивыбросов и выбросов используют критерий Граббса, описанный в 7.3.4 ГОСТ Р ИСО 5725-2.
Для контроля наличия квазивыбросов и выбросов во внутриэлементных расхождениях, применяют тестирование по критерию Граббса к значениям в каждой графе таблицы 2 по очереди.
Для контроля наличия квазивыбросов и выбросов в средних значениях элементов применяют тестирование по критерию Граббса к значениям в каждой графе таблицы 3 по очереди.
Интерпретация результатов тестирования полностью рассмотрена в 7.3.2 ГОСТ Р ИСО 5725-2. Их используют для идентификации результатов, которые настолько не соответствуют остальным данным эксперимента, что в случае их включения в расчеты стандартных отклонений повторяемости и воспроизводимости они окажут существенное влияние на значения этих статистик. Обычно данные, идентифицированные как выбросы, исключают из расчетов, а данные, идентифицированные как квазивыбросы, включают в расчеты, если не имеется серьезных оснований для принятия других решений. Если результаты тестирования показывают, что данные в одной из таблиц 2 или 3 должны быть исключены из расчетов стандартных отклонений повторяемости и воспроизводимости, то соответствующие значения в другой таблице также должны быть исключены.
4.7 Представление результатов эксперимента
4.7.1 В 7.7 ГОСТ Р ИСО 5725-2 даны рекомендации по:
- созданию совета экспертов специально для организации эксперимента и рассмотрения его результатов;
- представлению результатов статистического анализа совету экспертов;
- решениям, принимаемым советом экспертов по результатам рассмотрения;
- подготовке полного отчета.
4.7.2 Рекомендации по форме представления установленных стандартных отклонений повторяемости и воспроизводимости стандартного метода измерений даны в 7.1 ГОСТ Р ИСО 5725-1.
4.8 Пример 1. Эксперимент с разделенными уровнями
4.8.1 Таблица 4 содержит данные эксперимента [2] по определению содержания протеина в кормах методом сжигания. Число лабораторий-участниц - девять, эксперимент содержал 14 уровней. В каждом уровне использовались две пробы кормов с одинаковой массовой долей протеина.
Таблица 4 - Пример 1. Определение массовой доли протеина в кормах (в процентах)
Номер лаборатории |
Уровень |
|||||||||||||||||||||||||||
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
||||||||||||||
|
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
а |
b |
1 |
11,11 |
10,34 |
10,91 |
9,81 |
13,74 |
13,48 |
13,79 |
13,00 |
15,89 |
15,26 |
20,14 |
19,78 |
20,33 |
20,06 |
46,45 |
44,42 |
52,05 |
49,40 |
65,84 |
59,14 |
84,16 |
80,86 |
85,38 |
81,71 |
87,64 |
88,23 |
90,24 |
82,10 |
2 |
11,12 |
9,94 |
11,38 |
10,31 |
14,00 |
13,12 |
13,44 |
13,06 |
15,69 |
15,10 |
19,25 |
20,25 |
20,36 |
19,94 |
46,69 |
44,62 |
51,94 |
48,81 |
66,31 |
59,19 |
84,50 |
81,06 |
85,56 |
82,44 |
88,81 |
88,38 |
89,88 |
81,44 |
3 |
11,26 |
10,46 |
10,95 |
10,51 |
13,38 |
12,70 |
13,54 |
13,18 |
15,83 |
15,73 |
20,48 |
19,86 |
20,56 |
20,11 |
46,90 |
44,56 |
52,18 |
48,90 |
66,06 |
58,52 |
82,26 |
79,43 |
85,26 |
82,15 |
88,58 |
88,12 |
89,48 |
81,67 |
4 |
11,07 |
10,41 |
11,66 |
9,95 |
13,01 |
13,16 |
13,58 |
12,88 |
15,08 |
15,63 |
21,54 |
20,06 |
20,64 |
20,46 |
47,13 |
45,29 |
51,73 |
48,56 |
65,93 |
58,93 |
84,39 |
80,08 |
85,20 |
81,76 |
88,47 |
87,98 |
90,04 |
80,73 |
5 |
10,69 |
10,31 |
10,98 |
10,13 |
13,24 |
13,33 |
13,32 |
12,59 |
15,02 |
14,90 |
19,90 |
19,66 |
20,56 |
19,24 |
45,83 |
43,73 |
50,84 |
47,91 |
64,19 |
57,94 |
81,71 |
79,01 |
83,58 |
79,74 |
86,43 |
86,19 |
88,59 |
80,46 |
6 |
11,73 |
11,01 |
12,31 |
10,92 |
14,01 |
13,66 |
14,04 |
13,64 |
16,43 |
15,94 |
20,31 |
20,27 |
20,85 |
20,63 |
46,86 |
43,96 |
52,18 |
49,03 |
65,73 |
58,77 |
82,85 |
81,16 |
84,44 |
80,90 |
87,78 |
86,89 |
89,40 |
80,88 |
7 |
11,13 |
10,36 |
11,38 |
10,44 |
12,94 |
12,44 |
13,63 |
13,06 |
15,75 |
15,56 |
20,00 |
20,56 |
20,25 |
20,19 |
46,25 |
44,31 |
52,25 |
49,44 |
66,06 |
59,19 |
86,25 |
81,00 |
84,88 |
81,44 |
88,06 |
88,00 |
89,31 |
81,38 |
8 |
11,21 |
10,51 |
11,32 |
10,84 |
13,09 |
13,76 |
13,85 |
13,49 |
15,98 |
15,89 |
20,43 |
20,69 |
20,85 |
20,27 |
47,11 |
44,40 |
52,44 |
48,81 |
65,66 |
59,38 |
84,59 |
81,16 |
84,96 |
81,71 |
88,50 |
87,98 |
89,94 |
81,56 |
9 |
11,80 |
11,21 |
11,35 |
9,88 |
13,85 |
14,46 |
13,96 |
13,77 |
16,51 |
15,72 |
20,64 |
21,01 |
20,78 |
20,89 |
47,09 |
45,15 |
52,19 |
48,46 |
66,33 |
59,47 |
83,05 |
80,93 |
84,73 |
81,94 |
88,24 |
88,05 |
89,75 |
81,35 |