ГОСТ Р ИСО 11843-6-2015 Статистические методы. Способность обнаружения. Часть 6. Методология определения критического значения и минимального обнаруживаемого значения с применением аппроксимации распределения Пуассона нормальным распределением

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

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

Статистические методы

СПОСОБНОСТЬ ОБНАРУЖЕНИЯ

Часть 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

ОКС 03.120.30

17.020*

_____________________

* По данным официального сайта Росстандарта ОКС 03.120.30,

здесь и далее. - .

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

Предисловие

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

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

3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 6 октября 2015 г. N 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-2012 (пункт 3.5).

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

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

6 ПЕРЕИЗДАНИЕ. Ноябрь 2019 г.

Правила применения настоящего стандарта установлены в статье 26 Федерального закона от 29 июня 2015 г. N 162-ФЗ "О стандартизации в Российской Федерации". Информация об изменениях к настоящему стандарту публикуется в ежегодном (по состоянию на 1 января текущего года) информационном указателе "Национальные стандарты", а официальный текст изменений и поправок - в ежемесячном информационном указателе "Национальные стандарты". В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ближайшем выпуске ежемесячного информационного указателя "Национальные стандарты". Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования - на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет (www.gost.ru)

Введение

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

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

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

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

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

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

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

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

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

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

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

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

В настоящем стандарте использованы нормативные ссылки на следующие стандарты:

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

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

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

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

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

ISO 11843-4, Capability of detection - Part 4: Methodology for comparing the minimum detectable value with a given value (Способность обнаружения. Часть 4. Метод сравнения заданного значения с минимальным обнаруживаемым значением )

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)

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

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

, (2)

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

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

- среднее результатов измерений отклика в базовом состоянии;

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

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

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

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

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

X - отклик; Y - базовое состояние; - вероятность ошибки первого рода; - вероятность ошибки второго рода

Рисунок 1 - Принципиальное соотношение критического значения и результатов измерений отклика в фактическом и базовом состояниях

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

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

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

, (3)

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

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

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

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

. (4)

Если заменить на [см. формулы (2) и (3)], то:

, (5)

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

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

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

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

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

. (6)

Если заменить оценкой (см. 5.2), а также заменить оценкой (см. приложение B), критерий превращается в неравенство (7).

. (7)

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

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

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

Доверительный интервал для определяют по N результатам повторных измерений в базовом состоянии и N повторных измерений выборки с переменной состояния равной . Доверительный интервал с уровнем доверия 100()% для ()

, (8)

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

Для подтверждения критерия достаточной способности обнаружения используют односторонний критерий. Нижняя граница одностороннего доверительного интервала для () с и уровнем доверия 100()% имеет вид

, (9)

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

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

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

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

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

. (10)

Приближенную нижнюю границу доверительного интервала для () с уровнем доверия 100()% получают, заменяя и на и , соответственно [как в (3) и (7)]:

. (11)

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

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

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

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

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

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

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

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

;

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

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

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

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

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

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

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

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

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

, (B.1)

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

у - фактическое количество событий, зарегистрированных за то же время.

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

. (B.2)

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

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

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

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

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

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

Это означает, что распределение можно описать формулой (C.1), где .

. (C.1)

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

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

. (C.2)

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

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

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

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

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

Продолжение таблицы С.1

Окончание таблицы С.1

X - значение ; Y - разность значений , полученных точным и приближенным методом в %

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

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

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

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

X - каналы; Y - относительная интенсивность (учитываемых импульсов); S - область сигнала; B1 - левая фоновая область; B2 - правая фоновая область

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

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

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

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

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

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

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

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

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

- ширина канала: 0,02 градуса Брэгговского угла 2;

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

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

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

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

X - Брэгговский угол 2; Y - относительная интенсивность сигналов

Рисунок E.1 - Образец дифракции рентгеновских лучей асбеста хризотилового

E.1.1 Статистический анализ

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

.

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

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

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

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

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

,

,

.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

X - энергия связи, эВ; Y - относительная интенсивность сигналов

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

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

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

.

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

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

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

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

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

Сведения о соответствии ссылочных международных стандартов национальным и межгосударственным стандартам

Таблица ДА.1

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

Электронный текст документа

и сверен по:

, 2019