ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ТЕХНИЧЕСКОМУ РЕГУЛИРОВАНИЮ И МЕТРОЛОГИИ

НАЦИОНАЛЬНЫМ

СТАНДАРТ

РОССИЙСКОЙ

ФЕДЕРАЦИИ

Статистические методы СПОСОБНОСТЬ ОБНАРУЖЕНИЯ

Часть 6

Методология определения критического значения и минимального обнаруживаемого значения с применением аппроксимации распределения Пуассона нормальным распределением

(ISO 11843-6:2013, ЮТ)

Издание официальное

Москва

Стандартинформ

2016

Предисловие

1    ПОДГОТОВЛЕН Открытым акционерным обществом «Научно-исследовательский центр контроля и диагностики технических систем» (АО «НИЦ КД») на основе собственного аутентичного перевода на русский язык англоязычной версии международного стандарта, указанного в пункте 4

2    ВНЕСЕН Техническим комитетом по стандартизации ТК 125 «Применение статистических методов»

3    УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 6 октября 2015 г. № 1471-ст

4    Настоящий стандарт идентичен международному стандарту ИСО 11843-6:2013 «Способность обнаружения. Часть 6. Методология определения критического значения и минимального обнаруживаемого значения с применением аппроксимации распределения Пуассона нормальным распределением» (ISO 11843-6:2013 «Capability of detection — Part 6: Methodology for the determination of the critical value and the minimum detectable value in Poisson distributed measurements by normal approximations», IDT).

Наименование настоящего стандарта изменено относительно наименования указанного международного стандарта для приведения в соответствие с ГОСТ Р 1.5 (подраздел 3.5).

При применении настоящего стандарта рекомендуется использовать вместо ссылочных международных стандартов соответствующие им национальные стандарты Российской Федерации и межгосударственные стандарты, сведения о которых приведены в дополнительном приложении ДА

5    ВВЕДЕН ВПЕРВЫЕ

Правила применения настоящего стандарта установлены в ГОСТ Р 1.0-2012 (раздел 8). Информация об изменениях к настоящему стандарту публикуется в ежегодном (по состоянию на 1 января текущего года) информационном указателе «Национальные стандарты», а официальный текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ближайшем выпуске ежемесячного информационного указателя «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет (www.gost.ru)

© Стандартинформ, 2016

Настоящий стандарт не может быть полностью или частично воспроизведен, тиражирован и распространен в качестве официального издания без разрешения Федерального агентства по техническому регулированию и метрологии

ГОСТ Р ИСО 11843-6-2015

Приложение А (справочное)

X — У — J —

к —

N —

х — У — Ус -

Обозначения, используемые в ИСО 11843-6

параметр состояния; переменная отклика;

количество повторных измерений на образце сравнения, представляющих значения в базовом состоянии (пустая проба);

количество повторных измерений в реальном состоянии (контрольная выборка);

количество повторных измерений каждого образца сравнения при определении оценки способности обнаружения;

значение переменной состояния; значение переменной отклика;

критическое значение переменной отклика в соответствии с ИСО 11843-1 и ИСО 11843-3;

заданное значение, которое проверяют для определения, больше оно минимального обнаруживаемого значения или нет;

°ь -°д -

Ль — Лд -

Уь -

?д -

Уь -х —

а —

р -

1—а —

1-Р -

^z1^ — ^z1-|3

То -

минимальное обнаруживаемое значение переменной состояния; стандартное отклонение отклика в фактических условиях в базовом состоянии;

стандартное отклонение отклика в фактических условиях для выборки с переменной состояния, равной

^х9;

математическое ожидание отклика в фактических условиях в базовом состоянии;

математическое ожидание отклика в фактических условиях для выборки с переменной состояния, равной

среднее арифметическое результатов измерений отклика в базовом состоянии по выборке с переменной состояния, равной х;

среднее арифметическое результатов измерений отклика в базовом состоянии по выборке с переменной состояния, равной х_д;

минимальное обнаруживаемое значение отклика с переменной состояния, равной x_d;

среднее, соответствующее среднему количеству событий в распределении Пуассона;

вероятность ошибки первого рода;

вероятность ошибки второго рода;

доверительная вероятность;

доверительная вероятность;

квантиль нормированного нормального распределения уровня (1-а); квантиль нормированного нормального распределения уровня (1-р); нижняя доверительная граница.

7

Приложение В (справочное)

Определение оценок среднего и дисперсии при аппроксимации распределения Пуассона

нормальным распределением

Функция вероятности распределения Пуассона р(у X) имеет вид:

Р(уЛ) = ^-е-\    (В.1)

У!

где X — среднее, соответствующее среднему количеству событий за установленный период времени; у — фактическое количество событий, зарегистрированных за то же время.

Так как случайная величина У подчиняется распределению Пуассона с параметром X, математическое ожидание и дисперсия этой случайной величины равны X, т. е. £(У) = X и Var(Y) =Х. Таким образом, необходимо определить оценку только параметра X. Эта оценка на основе J независимых результатов измерений имеет вид:

j

При аппроксимации распределения Пуассона нормальным распределением случайную величину У заменяют случайной величиной Z, которая подчиняется нормальному распределению Л/(Х, X).

8

ГОСТ Р ИСО 11843-6-2015

Приложение С (справочное)

Точность аппроксимации

В данном приложении минимальные обнаруживаемые значения отклика, вычисленные с использованием аппроксимации распределения Пуассона нормальным распределением, сопоставлены с минимальными обнаруживаемыми значениями, полученными на основе распределения Пуассона. Это позволяет определить зависимость точности аппроксимации от количества данных.

Минимальное обнаруживаемое значение отклика для распределения Пуассона вычисляют в соответствии со следующей процедурой.

Сумма случайных величин, подчиняющихся распределению Пуассона, также подчиняется распределению Пуассона. Однако это несправедливо для разности случайных величин, подчиняющихся распределению Пуассона. Если эта разность четко описана, используется следующая функция вероятности. Отклик в базовом состоянии соответствует фоновому шуму при измерениях у_ь и переменная отклика в реальном состоянии y_d отражает каждую из двух выборок в условиях нулевой гипотезы.

Это означает, что распределение можно описать формулой (С.1), где у = |y_b-y_d|.

оо

(С.1)

Рг[у] = е-^2е Е %^Mhy) [j!(/- у)!]"¹ = е-²⁰ /_у( 20).

/_к(-) — модифицированная функция Бесселя первого вида.

Распределение можно описать формулой (С.2) в условиях альтернативной гипотезы.

/ 01 \ у/2

_Рг[у] = _е-(⁰1 ^{+ 0}2) L 0/ 0₂^J-^y [7! (J - у)! Г¹ = е-<⁰1⁺ е2> (-^-) 1_у ( 2    (С.2)

Минимальный обнаруживаемый отклик в фактическом состоянии может быть получен из этих двух уравнений. Альтернативно минимальный обнаруживаемый отклик при использовании аппроксимации может быть получен с помощью формул (7) и (11), если количество повторных измерений N заменить на бесконечность.

В таблице С.1 приведено минимальное обнаруживаемое значение, когда параметр у_ь соответствует значению в базовом состоянии от 1 до 200, а также разности значений y_d, вычисленных с помощью точного и приближенного метода [11].

Вычисления для распределения Пуассона и при использовании нормальной аппроксимации дают довольно близкие результаты по всему диапазону (см. рисунок С.1).

Если минимальное обнаруживаемое значение отклика должно быть определено с погрешностью не более 5 %, условия измерений должны быть отрегулированы так, чтобы значение у_ь составляло не менее 18.

Таблица С.1 — Сопоставление распределения Пуассона и нормальной аппроксимации

Значение отклика в базовом состоянии Распреде ление Пуассона Нормальное прибли жение Разность Значение отклика в базовом состоянии Распределение Пуассона Нормальное прибли жение Разность Уь Уь Уь Уь Уь Уь 1 8,2 8,4 -0,1 8 24,7 23,9 0,9 2 11,3 11,3 0,0 9 26,1 25,7 0,4 3 14,1 13,8 0,3 10 27,4 27,4 0,0 4 17,1 16,0 1,0 11 29,9 29,1 0,7 5 18,9 18,1 0,8 12 31,2 30,8 0,3 6 20,8 20,1 0,7 13 32,5 32,5 0,0 7 22,2 22,0 0,2 14 34,9 34,1 0,7

9

Значение отклика в базовом состоянии	Распреде ление Пуассона	Нормальное прибли жение	Разность	Значение отклика в базовом состоянии	Распределение Пуассона	Нормальное прибли жение	Разность
Уь	Уь	Уь		Уь	yd	yd
15	36,1	35,7	0,4	50	85,6	85,6	0,0
16	37,4	37,3	0,1	51	87,8	86,9	0,9
17	39,8	38,9	0,9	52	88,9	88,3	0,7
18	41,0	40,4	0,6	53	90,1	89,6	0,5
19	42,3	42,0	0,3	54	91,2	90,9	0,3
20	43,5	43,5	0,0	55	92,4	92,2	0,2
21	45,8	45,0	0,8	56	93,5	93,5	0,0
22	47,1	46,5	0,5	57	95,7	94,8	0,9
23	48,3	48,0	0,3	58	96,9	96,1	0,8
24	49,5	49,5	0,0	59	98,0	97,4	0,6
25	51,8	51,0	0,8	60	99,2	98,7	0,4
26	53,0	52,4	0,6	61	100,3	100,0	0,3
27	54,2	53,9	0,3	62	101,5	101,3	0,1
28	55,4	55,3	0,1	63	102,6	102,6	0,0
29	57,7	56,8	1,0	64	104,8	103,9	0,9
30	58,9	58,2	0,7	65	105,9	105,2	0,7
31	60,1	59,6	0,5	66	107,1	106,5	0,6
32	61,3	61,0	0,3	67	108,2	107,8	0,4
33	62,5	62,4	0,0	68	109,3	109,1	0,3
34	64,7	63,8	0,9	69	110,5	110,4	0,1
35	65,9	65,2	0,7	70	111,6	111,6	0,0
36	67,1	66,6	0,5	71	113,8	112,9	0,9
37	68,3	68,0	0,3	72	114,9	114,2	0,7
38	69,5	69,4	0,1	73	116,0	115,5	0,6
39	71,7	70,8	1,0	74	117,2	116,7	0,4
40	72,9	72,1	0,8	75	118,3	118,0	0,3
41	74,1	73,5	0,6	76	119,4	119,3	0,2
42	75,2	74,9	0,4	77	120,5	120,5	0,0
43	76,4	76,2	0,2	78	122,7	121,8	0,9
44	77,5	77,6	0,0	79	123,9	123,1	0,8
45	79,8	78,9	0,9	80	125,0	124,3	0,7
46	80,9	80,3	0,7	81	126,1	125,6	0,5
47	82,1	81,6	0,5	82	127,2	126,8	0,4
48	83,3	82,9	0,3	83	128,3	128,1	0,2
49	84,4	84,3	0,1	84	129,5	129,3	0,1

Значение отклика в базовом состоянии	Распреде ление Пуассона	Нормальное прибли жение	Разность	Значение отклика в базовом состоянии	Распределение Пуассона	Нормальное прибли жение	Разность
Уь	Уь	Уь		Уь	yd	yd
85	130,6	130,6	0,0	120	173,6	173,7	0,0
86	132,8	131,9	0,9	121	175,8	174,9	0,9
87	133,9	133,1	0,8	122	176,9	176,1	0,8
88	135,0	134,3	0,6	123	178,0	177,3	0,7
89	136,1	135,6	0,5	124	179,1	178,5	0,6
90	137,2	136,8	0,4	125	180,2	179,7	0,4
91	138,3	138,1	0,3	126	181,3	180,9	0,3
92	139,5	139,3	0,1	127	182,4	182,1	0,2
93	140,6	140,6	0,0	128	183,5	183,3	0,1
94	142,7	141,8	0,9	129	184,6	184,5	0,0
95	143,9	143,1	0,8	130	186,7	185,8	1,0
96	145,0	144,3	0,7	131	187,8	187,0	0,9
97	146,1	145,5	0,6	132	188,9	188,2	0,8
98	147,2	146,8	0,4	133	190,0	189,4	0,6
99	148,3	148,0	0,3	134	191,1	190,6	0,5
100	149,4	149,2	0,2	135	192,2	191,8	0,4
101	150,5	150,5	0,1	136	193,3	193,0	0,3
102	151,6	151,7	-0,1	137	194,4	194,2	0,2
103	153,8	152,9	0,9	138	195,5	195,4	0,1
104	154,9	154,2	0,7	139	196,6	196,6	0,0
105	156,0	155,4	0,6	140	198,7	197,8	1,0
106	157,1	156,6	0,5	141	199,8	198,9	0,9
107	158,2	157,8	0,4	142	200,9	200,1	0,8
108	159,3	159,1	0,3	143	202,0	201,3	0,6
109	160,4	160,3	0,2	144	203,1	202,5	0,6
110	161,5	161,5	0,0	145	204,2	203,7	0,5
111	163,7	162,7	1,0	146	205,3	204,9	0,3
112	164,8	163,9	0,9	147	206,4	206,1	0,2
113	165,9	165,2	0,7	148	207,5	207,3	0,1
114	167,0	166,4	0,6	149	208,6	208,5	0,1
115	168,1	167,6	0,5	150	209,6	209,7	0,0
116	169,2	168,8	0,4	151	211,8	210,9	0,9
117	170,3	170,0	0,3	152	212,9	212,1	0,8
118	171,4	171,2	0,2	153	214,0	213,3	0,7
119	172,5	172,5	0,1	154	215,0	214,4	0,6

Значение отклика в базовом состоянии	Распреде ление Пуассона	Нормальное прибли жение	Разность	Значение отклика в базовом состоянии	Распределение Пуассона	Нормальное прибли жение	Разность
Уь	Уь	Уь		Уь	yd	yd
155	216,1	215,6	0,5	178	243,2	242,8	0,4
156	217,2	216,8	0,4	179	244,3	244,0	0,3
157	218,3	218,0	0,3	180	245,4	245,1	0,2
158	219,4	219,2	0,2	181	246,5	246,3	0,2
159	220,5	220,4	0,1	182	247,5	247,5	0,1
160	221,6	221,6	0,0	183	248,6	248,6	0,0
161	223,7	222,7	1,0	184	250,7	249,8	0,9
162	224,8	223,9	0,9	185	251,8	251,0	0,8
163	225,9	225,1	0,8	186	252,9	252,2	0,8
164	227,0	226,3	0,7	187	254,0	253,3	0,7
165	228,1	227,5	0,6	188	255,1	254,5	0,6
166	229,1	228,6	0,5	189	256,2	255,7	0,5
167	230,2	229,8	0,4	190	257,2	256,8	0,4
168	231,3	231,0	0,3	191	258,3	258,0	0,3
169	232,4	232,2	0,2	192	259,4	259,2	0,2
170	233,5	233,4	0,1	193	260,5	260,3	0,1
171	234,6	234,5	0,0	194	261,6	261,5	0,1
172	236,7	235,7	1,0	195	262,6	262,7	0,0
173	237,8	236,9	0,9	196	264,8	263,8	0,9
174	238,9	238,1	0,8	197	265,8	265,0	0,8
175	240,0	239,3	0,7	198	266,9	266,2	0,7
176	241,0	240,4	0,6	199	268,0	267,3	0,7
177	242,1	241,6	0,5	200	269,1	268,5	0,6

12

ГОСТ РИСО 11843-6—2015

х

Рисунок С.1 — Разность значений y_d> полученных для распределения Пуассона и с использованием его аппроксимации нормальным распределением в процентах

Приложение D (справочное)

Выбор количества каналов для датчика

Количество каналов датчика, выбранных для измерений, определяет диапазон, который может быть охвачен в процессе измерений. Важно обеспечить отсутствие наложения соседних пиков и одинаковое количество каналов при измерении фонового шума и выборки.

Рисунок D.1 —Требования к областям сигналов

13

Приложение Е (справочное)

Примеры расчетов

Е.1 Пример 1. Определение содержания опасных веществ с помощью рентгеновской

дифрактометрии

Асбест хризотиловый (опасное вещество) может быть обнаружен при использовании рентгеновской дифрактометрии. Образец асбеста хризотилового массой 0,10 мг был точно взвешен, смешан с чистой водой, осажден на фильтровальную бумагу и проанализирован. Результаты показали содержание асбеста 0,10 % или 0,10 мг/100 мг (максимально допустимое значение в строительных материалах). Проводилось пять повторных измерений содержания асбеста в холостой пробе с х_ь = 0 и в фактической пробе с х_д. Результаты исследования порошка асбеста хризотилового, представленные в виде зависимости интенсивности дифрагированного излучения от Брэгговского угла 29, приведены на рисунке Е.1. Способность обнаружения вычислена для K=J=1 и а = р = 0,05.

Условия измерений:

Средство измерений: рентгеновский дифрактометр;

Источник рентгеновского излучения: Си, обеспечивающий К-альфа монохроматическое рентгеновское излучение с 0,154 нм;

Выходная мощность: 40 кВ, 40 мА;

Ширина канала: 0,02 градуса Брэгговского угла 20;

Общее количество каналов: 23 (для фонового шума и для пиковой области);

Количество измерений: 5;

Время аккумуляции в каждом канале: 2 с;

Полная ширина на уровне половины максимума (FWHM); 0,46 градуса Брэгговского угла 20.

X — Брэгговский угол 20; Y — относительная интенсивность сигналов Рисунок Е.1 — Образец дифракции рентгеновских лучей асбеста хризотилового Е.1.1 Статистический анализ

В результате эксперимента у_ь = 174, у_д = 261. Для этих значений нижняя граница доверительного интервала, вычисленная по формуле (11) для уровня доверия 95 % имеет вид;

<У,~9ь) ^{+ г}(1-)^    (261-174)-1,645

14

ГОСТ Р ИСО 11843-6-2015

Эта нижняя граница больше

74 + 261 = 65,0

(см. формулы (6) и (7)).

Таким образом, с уровнем доверия 95 % можно принять решение о том, что в исследованном порошке присутствует 0,1 % асбеста хризотилового.

Е.1.2 Определение оценки минимальной обнаруживаемой концентрации асбеста

Минимальная обнаруживаемая концентрация асбеста в строительных материалах может быть определена с использованием тех же самых условий. Минимальная обнаруживаемая концентрация может быть улучшена при увеличении количества измерений. Увеличение количества измерений до бесконечности дает предельную минимальную обнаруживаемую концентрацию, как показано в следующих вычислениях:

yd - 174 = 1,645- (V2 ■ 174 + V174 + yd)

Было получено 238 сигналов для минимального обнаруживаемого значения переменной отклика, по отношению к которому содержание асбеста может быть оценено в процентах следующим образом:

0,1 % / (261 - 174) (количество сигналов) = 1,15 • 10“³ % / (количество сигналов).

Предельное минимальное обнаруживаемое содержание асбеста хризотилового получено следующим образом:

1,15 • 10^-3 % • (238 - 174) (количество сигналов) = 0,074 %.

Таким образом, оценка минимального обнаруживаемого значения содержания асбеста хризотилового составляет 0,074 %.

Е.2 Пример 2. Определение загрязнения на поверхности кремниевых пластин с помощью

рентгеновской фотоэлектронной спектроскопии

Загрязнение поверхности кремниевой пластины веществами органического происхождения определено на основе измерения энергии связи на уровне Is орбитали углерода, полученной с помощью фотоэлектронной спектроскопии.

Полученные значения для трех повторных измерений в двух областях одного и того же образца, представляющих х_ь = 0 и неизвестное содержание х_д, приведены в таблице Е.1 вместе с Is-спектром углерода (рисунок Е.2), полученным на основе фотоэлектронной спектроскопии.

При использовании трех повторных измерений способность обнаружения вычислена для К = J = 1 и а = р = 0,05.

Условия измерений:

Средство измерений: рентгеновский фотоэлектронный спектрометр;

Источник рентгеновского излучения: AI, обеспечивающий К-альфа монохроматическое рентгеновское излучение;

Выходная мощность: 0,9 Вт;

Ширина канала 0,025 эВ;

Количество каналов: 11 (для фонового шума и пиковой области);

Количество измерений: 3;

Время аккумуляции в каждом канале: 100 мс.

15

Таблица Е.1 — Наблюдаемое количество сигналов для фактической и холостой пробы

Фоновый шум Количество сигналов Пиковая область Количество сигналов эВ п = 1 п = 2 п = 3 Среднее эВ п = 1 п = 2 /7=3 Среднее 291,85 102 78 81 283,98 111 98 113 291,83 99 77 87 283,95 98 103 111 291,80 96 64 88 283,93 107 103 128 291,78 112 85 72 283,90 99 104 91 291,75 100 86 90 283,88 111 95 109 291,73 99 88 63 283,85 108 131 110 291,70 109 81 76 283,83 103 95 87 291,68 95 76 98 283,80 121 115 118 291,65 91 79 75 283,78 97 116 91 291,63 107 90 78 283,75 102 103 102 291,60 92 90 72 283,73 118 95 105 у, (суммарное количество сигналов/ канал) 1102 894 880 959 у, (суммарное количество сигналов/канал) 1175 1158 1165 1166

У

Рисунок Е.2 — Is-спектр углерода, полученный на основе фотоэлектронной спектроскопии

16

ГОСТ Р ИСО 11843-6-2015

Содержание

1    Область применения.................................................................1

2    Нормативные ссылки.................................................................1

3    Термины и определения ..............................................................2

4    Измерительная система и обработка данных.............................................2

5    Аппроксимация распределения Пуассона................................................2

6    Отчет о полученной оценке способности обнаружения.....................................6

7    Отчет о результатах применения метода ................................................6

Приложение А (справочное) Обозначения, используемые в ИСО 11843-6 .......................7

Приложение В (справочное) Определение оценок среднего и дисперсии при аппроксимации

распределения Пуассона нормальным распределением ........................8

Приложение С (справочное) Точность аппроксимации.......................................9

Приложение D (справочное) Выбор количества каналов для датчика .........................13

Приложение Е (справочное) Примеры расчетов...........................................14

Приложение ДА (справочное) Сведения о соответствии ссылочных международных стандартов национальным стандартам Российской Федерации

и действующим в этом качестве межгосударственным стандартам..............18

Библиография.......................................................................19

ГОСТ Р ИСО 11843-6-2015

Е.2.1 Статистический анализ

В результате эксперимента у_ь = 959, у_д = 1166. Используя эти значения по формуле (10), определена нижняя граница доверительного интервала для уровня доверия 95 %:

Эта нижняя граница больше

(см. неравенство (6)). Поэтому с уровнем доверия 95 % можно утверждать, что поверхность кремниевой пластины загрязнена углеводородом.

Е.2.2 Определение минимального обнаруживаемого содержания углеводорода с использованием аппроксимации

Минимальное обнаруживаемое содержание загрязнения углеводородом поверхности кремниевой пластины получено на основе наблюдаемого количества сигналов всех составляющих элементов и факторов чувствительности установленного средства измерений.

17

Введение

В приборах различного типа для обнаружения сигналов часто используют метод подсчета импульсов (датчик рентгеновской, электронной и ионной спектроскопии, такие как рентгеновские дифрактометры, рентгеновские флуоресцентные спектрометры, рентгеновские фотоэлектронные спектрометры, оже-электронные спектрометры, масс-спектрометры вторичной ионизации, хроматомасспектрометры и др.). В этом случае сигналы состоят из серии импульсов, появляющихся случайным образом через нерегулярные интервалы времени. Для их описания может быть применено распределение Пуассона, а методология определения минимального обнаруживаемого значения может быть выведена с применением статистической теории.

Определение минимального обнаруживаемого значения иногда важно в практической работе. Это значение обеспечивает критерий принятия решения, когда «сигнал не обнаружен или когда сигнал значимо отличается от фонового шума» [1—8]. Например, это важно, когда измеряемая величина представляет опасные вещества или поверхностное загрязнение полупроводниковых материалов. Ограничения на опасные вещества устанавливают пределы использования шести опасных материалов (шестивалентный хром, свинец, ртуть, кадмий и огнезащитные агенты, пербромодифенил) в производстве электронных компонентов и соответствующих товаров, продаваемых в ЕС. Для проверки используют специальные приборы. Рентгеновские флуоресцентные спектрометры и газохроматографические масс-спектрометры являются такими приборами. Рентгеновские дифрактометры используют для измерения уровня опасности асбеста и кристаллического кремния в окружающей среде или в строительных материалах.

Методы, используемые для определения минимального обнаруживаемого значения, должны в течение некоторого времени широко использоваться в области аналитической химии, но не там, где используют метод подсчета импульсов. Необходимо установить методологию определения минимального обнаруживаемого значения в этой области.

В настоящем стандарте распределение Пуассона аппроксимировано нормальным распределением, обеспечивающим согласованность с подходом IUPAC, лежащем в основе серии стандартов ISO 11843. Обычную аппроксимацию используют для генерации дисперсии, критического значения отклика, критерия способности обнаружения и минимального обнаруживаемого значения [1].

В настоящем стандарте:

а — вероятность ошибочного обнаружения того, что система не находится в базовом состоянии, если в действительности она находится в базовом состоянии;

(3 — вероятность ошибочного не обнаружения того, что система не находится в базовом состоянии, если значение переменной состояния равно минимальному обнаруживаемому значению (x_d).

Настоящий стандарт полностью согласован с ИСО 11843-1, ИСО 11843-3, ИСО 11843-4.

IV

НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ

Статистические методы СПОСОБНОСТЬ ОБНАРУЖЕНИЯ Часть 6

Методология определения критического значения и минимального обнаруживаемого значения с применением аппроксимации распределения Пуассона нормальным распределением

Statistical methods. Capability of detection. Part 6. Methodology for the determination of the critical value and the minimum detectable value in Poisson distributed measurements by normal approximations

Дата введения — 2016—12—01

1    Область применения

В настоящем стандарте установлены методики определения критического значения переменной отклика и минимального обнаруживаемого значения по результатам измерений, подчиняющихся распределению Пуассона. Стандарт применим в тех случаях, когда изменчивость фонового шума и сигнала могут быть описаны распределением Пуассона. Обычно для распределения Пуассона используют аппроксимацию в виде нормального распределения (см. ИСО 11843-3 и ИСО 11843-4).

Сопоставление точности результатов при использовании распределения Пуассона и его аппроксимации нормальным распределением приведено в приложении С.

2    Нормативные ссылки

В настоящем стандарте использованы нормативные ссылки на следующие стандарты: Руководство ИСО 30 Термины и определения, используемые в области стандартных образцов (ISO Guide 30 Terms and definitions used in connection with reference materials)

ИСО 3534-1 Статистика. Словарь и условные обозначения. Часть 1. Общие статистические термины и термины, используемые в вероятностных задачах (ISO 3534-1 Statistics — Vocabulary and symbols — Part 1: General statistical terms and terms used in probability)

ИСО 11843-1 Способность обнаружения. Часть 1. Термины и определения (ISO 11843-1 Capability of detection — Part 1: Terms and definitions)

ИСО 11843-2 Способность обнаружения. Часть 2. Методология в случае линейной калибровки (ISO 11843-2 Capability of detection — Part 2: Methodology in the linear calibration case)

ИСО 11843-3 Способность обнаружения. Часть 3. Методология определения критического значения отклика без использования данных калибровки (ISO 11843-3 Capability of detection — Part 3: Methodology for determination of the critical value forthe response variable when no calibration data are used) ИСО 11843-4 Способность обнаружения. Часть 4. Метод сравнения заданного значения с минимальным обнаруживаемым значением (ISO 11843-4 Capability of detection — Part 4: Methodology for comparing the minimum detectable value with a given value)

Издание официальное

3    Термины и определения

В настоящем стандарте применены термины по ИСО 3534-1, ИСО 11843-1, ИСО 11843-2, ИСО 11843-3, ИСО 11843-4 и Руководству ИСО 30.

4    Измерительная система и обработка данных

Условия использования в расчетах распределения Пуассона обычно определяют экспериментально. Количество обнаруженных импульсов возрастает при увеличении времени наблюдений и ширины области наблюдаемого спектра. Эти два параметра должны быть установлены и сохраняться в процессе измерений.

Для достоверного определения минимального обнаруживаемого значения необходимо выполнение следующих требований:

a)    сигнал и фоновый шум подчиняются распределению Пуассона. Сигнал является средним наблюдаемых значений;

b)    исходные данные не обработаны и не сглажены;

c)    продолжительный непрерывный интервал времени измерений более предпочтителен, чем несколько коротких промежутков времени наблюдений. Измерения продолжительностью в одну секунду предпочтительнее, чем 10 измерений продолжительностью по 100 мс каждое. Аппроксимация распределения Пуассона нормальным распределением дает более точные результаты для более высоких значений среднего;

d)    количество измерений. Поскольку в приведенных методах использованы только средние, для их определения необходимо выполнение повторных измерений. Мощность теста возрастает с увеличением количества измерений;

e)    количество каналов, используемых датчиком. Не должно быть наложения соседних пиков. Количество используемых каналов для измерения фонового шума и выборочный спектр должны быть идентичными (приложение D, рисунок D.1);

f)    ширина пика. Для мониторинга единственного пика рекомендуется использовать полную ширину на половине максимума (FWHM). Это предпочтительно в случае измерения высоты вершины и/или ширины основания пика. Соответствующее значение FWHM должно быть определено заранее на основе результатов измерений стандартной выборки. Идентичное значение FWHM следует использовать и для фонового шума и для измерений элементов выборки.

Дополнительные факторы. Инструмент должен быть поверен (калиброван), датчик должен работать в пределах его линейного диапазона. Не должно быть сигналов, которые не могли бы быть четко идентифицированы, как не являющиеся шумом. Деградация объекта в процессе измерений должна быть незначительной. Должен быть обнаружен, по крайней мере, один сигнал или пик, относящийся к рассматриваемому объекту.

5    Аппроксимация распределения Пуассона

5.1 Критическое значение на основе нормального распределения

(1)

Решение о значимости измеренного сигнала (или нет) может быть принято на основе сопоставления среднего арифметического у_д результатов измерений с выбранным значением у_с. Значение у_с, рассматриваемое как критическое значение, удовлетворяет требованию

^р % ^> У_с I ^{х =} °) * «.

где вероятность вычислена при условии, что система находится в базовом состоянии (х = 0), а а — заданное значение вероятности.

Формула (1) задает вероятность того, что у_д >у_с при условии, что:

(2)

ГОСТ Р ИСО 11843-6-2015

где _Zl__a —

°ь —

% — J —

К —

квантиль нормированного нормального распределения уровня (1-а); стандартное отклонение при выполнении условий для отклика в базовом состоянии; среднее результатов измерений отклика в базовом состоянии;

количество повторных измерений чистой стандартной выборки. Оно представляет значение переменной в базовом состоянии;

количество повторных измерений контролируемой выборки. Оно дает значение фактической переменной состояния.

Знак «+» используют в формуле (2), если переменная отклика увеличивается при увеличении переменной состояния.

Знак «-» используют в противоположном случае.

X— отклик; У — базовое состояние; a — вероятность ошибки первого рода; |3 — вероятность ошибки второго рода Рисунок 1 — Принципиальное соотношение критического значения и результатов измерений отклика в фактическом и базовом состояниях

Определение критического значения соответствует ИСО 11843-1 и ИСО 11843-3. Его взаимосвязь с результатами измерений в фактическом и базовом состояниях показана на рисунке 1.

5.2 Определение критического значения переменной отклика

Если переменная отклика подчиняется распределению Пуассона с достаточно большим средним, оценка стандартного отклонения результатов повторных измерений переменной отклика в базовом состоянии (оценка о_ь) имеет вид Стандартное отклонение результатов повторных измерений переменной отклика в фактическом состоянии по выборке имеет вид и является оценкой a_g (см. приложение В).

з

Критическое значение у_с переменной отклика, которая подчиняется распределению Пуассона, аппроксимированному нормальным распределением, имеет вид:

Ус= Уь ^{+ 2}1-а⁰Ь ^ ⁺ X“^yb ^{+ Z}1-aV^^ ⁺ X’    ⁽³⁾

где у_ь — среднее арифметическое результатов измерений отклика в базовом состоянии.

5.3 Критерий достаточной способности обнаружения

Критерий достаточной способности обнаружения позволяет принимать решения об обнаружении сигнала на основе сопоставления критического значения вероятности с установленным значением уровня доверия (1-р). Если критерий выполнен, можно принимать решение о том, что минимальное обнаруживаемое значение x_d не превосходит значения переменной состояния х_д. Минимальное обнаруживаемое значение тогда определяет наименьшее значение переменной отклика п_д, для которого неправильное решение может быть принято с вероятностью (3. Это значение в ситуации, когда нет сигнала, а только фоновый шум, соответствует ошибке второго рода.

Если стандартное отклонение отклика для данного значения х_д равно а_д, критерий того, что вероятность больше или равна (1-р) устанавливает неравенство (4), из которого могут быть получены неравенства (5) и (6):

>

Ус + zi-p

+

2 д'

(4)

Если у_с заменить на у_с = ц_ь + z₁__a a_b

1/— + — (см. формулы (2) и (3)), то: f J К

Лд-ЛЬ^-а^Ь

Z1-P

+ 1 2 Т°д’

(5)

где

а — вероятность ошибки первого рода; р — вероятность ошибки второго рода; г)_ь — математическое ожидание отклика в базовом состоянии в фактических условиях; rig — математическое ожидание отклика по выборке с переменной состояния, равной х_д в фактических условиях.

Для р = а и К = J критерий можно упростить

Лд — ЛЬ ~ z1-a

(6)

Если о_ь заменить оценкой -^у^ (см. 5.2), а также а_д заменить оценкой yfjg (см. приложение В), критерий превращается в неравенство (7).

Лд — ЛЬ ~ z1-a

(7)

Примечание — При валидации метода в соответствии с ИСО 11843-4 обычно выбирают К = J = 1.

ГОСТ Р ИСО 11843-6-2015

5.4 Подтверждение критерия достаточной способности обнаружения

Стандартные отклонения и математические ожидания отклика обычно неизвестны, таким образом оценка при использовании неравенства критерия (6), должна быть получена по экспериментальным данным. Левая часть неравенства (6) неизвестна, а правая часть известна.

Доверительный интервал для ri_g -ц_ь определяют по N результатам повторных измерений в базовом состоянии и N повторных измерений выборки с переменной состояния равной х_д. Доверительный интервал с уровнем доверия 100(1-а/2) % для (ц_д - ц_ь):

(Уд-Уь)-²,

{b³

Ь ⁺

°о ⁵ Л_а

-Ль⁵ (Уд-Уь) ^{+ 2}|

(1 -а/2)

^°²ь ⁺ 1Т⁰1'

(8)

где z^__aJ2) — квантиль нормированного нормального распределения уровня 100 (1-а/2).

Для подтверждения критерия достаточной способности обнаружения используют односторонний критерий. Нижняя граница одностороннего доверительного интервала для (r)_g - г|_ь) с (5 = а и уровнем доверия 100(1-а) % имеет вид:

Лд Ль - (Уд Уь) ^{+ z}(1—а) У л/ ° Ь ⁺ л/ ° 9 ’    ^

где N — количество результатов повторных измерений каждого образца сравнения, используемого при определении оценки способности обнаружения; у_д — среднее арифметическое результатов измерений отклика по выборке с переменной состояния, равной х_д;

г)_ь — математическое ожидание отклика в базовом состоянии в фактических условиях; rig — математическое ожидание отклика по выборке с переменной состояния, равной х_д, в фактических условиях.

Нижнюю границу одностороннего доверительного интервала для (л_д - Ль)- полученную в соответствии с неравенством (9), сопоставляют с правой частью неравенства (6):

Лд - Ль * (Уд - Уь) ^{+ Z}(1-a) jjf °b ⁺ Jf °l ~ ^Z1-«    ⁺    M ⁺    ) •    (¹⁰)

Приближенную нижнюю границу T₀ доверительного интервала для (r)_g - г)_ь) с уровнем доверия 100(1—а) % получают, заменяя a_b и a_g на у/у^ и y/у^, соответственно (как в (3) и (7)):

^Г0 = (Уд " Уь) ^{+ Z}(1-a) УУЬ ⁺ Уд •    О¹¹)

Если Т₀ удовлетворяет неравенству (7), принимают решение о том, что минимальное среднее обнаруживаемое значение отклика у_д меньше или равно минимальному обнаруживаемому значению отклика у_& Значение x_d поэтому меньше или равно х и, для относительно больших значений N, нижняя доверительная граница в соответствии с формулой (11) является удовлетворительной.

5

ГОСТ Р ИСО 11843-6-2015

6 Отчет о полученной оценке способности обнаружения

Определение оценки способности обнаружения обычно является частью первичной валидации метода. Отчет должен включать:

a)    информацию об образцах сравнения, включая значения х_д, для эталонного состояния;

b)    количество повторных измерений N для каждого эталонного состояния;

c)    средние арифметические у_ь и у_д;

d)    выбранные значения а, (3, J и К;

e)    значения левых и правых частей неравенства (7) с использованием оценок, т. е. (у_д - у_ь) или, если это применимо, ((3 = а, К = J), значение (г|_д - г|_ь), соответствующий доверительный интервал и его приемлемую нижнюю границу

1) заключение относительно способности обнаружения;

д) при необходимости, минимальное обнаруживаемое значение для заданной величины фонового шума. Это значение получают, заменяя N и J бесконечностью и 1, соответственно, в формуле (10).

7 Отчет о результатах применения метода

Наблюдаемые значения должны быть зафиксированы, поскольку они представляют отклик переменной состояния. Тот факт, что наблюдаемые значения используют для проверки истинных значений, не является причиной для замены их верхней границей (равной критическому значению) или минимальным обнаруживаемым значением. В отчете также указывают критическое значение и минимальное обнаруживаемое значение.