Товары в корзине: 0 шт Оформить заказ
Стр. 1
 

24 страницы

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

 Скачать PDF

Идентичен ISO 11843-6:2013

Оглавление

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

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

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

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

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

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

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

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

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

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

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

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

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

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

 

24 страницы

Дата введения01.12.2016
Добавлен в базу01.02.2017
Актуализация01.01.2019

Этот ГОСТ находится в:

Организации:

06.10.2015УтвержденФедеральное агентство по техническому регулированию и метрологии1471-ст
ИзданСтандартинформ2016 г.
РазработанОАО НИЦ КД

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

Стр. 1
стр. 1
Стр. 2
стр. 2
Стр. 3
стр. 3
Стр. 4
стр. 4
Стр. 5
стр. 5
Стр. 6
стр. 6
Стр. 7
стр. 7
Стр. 8
стр. 8
Стр. 9
стр. 9
Стр. 10
стр. 10
Стр. 11
стр. 11
Стр. 12
стр. 12
Стр. 13
стр. 13
Стр. 14
стр. 14
Стр. 15
стр. 15
Стр. 16
стр. 16
Стр. 17
стр. 17
Стр. 18
стр. 18
Стр. 19
стр. 19
Стр. 20
стр. 20
Стр. 21
стр. 21
Стр. 22
стр. 22
Стр. 23
стр. 23
Стр. 24
стр. 24

ГОСТ Р исо

11843-6—

2015

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

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

СТАНДАРТ

РОССИЙСКОЙ

ФЕДЕРАЦИИ

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

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

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

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

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

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

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

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

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

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

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

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

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

7

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

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

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

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

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

У!

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

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

j

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

8

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

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

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

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

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

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

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

оо

(С.1)

Рг[у] = е- Е %Mhy) [j!(/- у)!]"1 = е-20 /у( 20).

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

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

/ 01 \ у/2

Рг[у] = е-(01 + 02) L 0/ 02J-y [7! (J - у)! Г1 = е-<01+ е2> (-^-) 1у ( 2    (С.2)

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

В таблице С.1 приведено минимальное обнаруживаемое значение, когда параметр уь соответствует значению в базовом состоянии от 1 до 200, а также разности значений yd, вычисленных с помощью точного и приближенного метода [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

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


х

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

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

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

У

X— каналы; У — относительная интенсивность (учитываемых импульсов); S — область сигнала; В1 —левая фоновая область;

В2 — правая фоновая область


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

Рисунок 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 % имеет вид;

— -^174 + 261 =71,7.

<У,~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“3 % / (количество сигналов).

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

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 — вероятность ошибочного не обнаружения того, что система не находится в базовом состоянии, если значение переменной состояния равно минимальному обнаруживаемому значению (xd).

Настоящий стандарт полностью согласован с ИСО 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 Определение критического значения переменной отклика

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

з

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

Ус= Уь + 21-а0Ь ^ + X“yb + Z1-aV^^ + X’    (3)

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

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

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

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


>


Ус + zi-p

+


2

д'


(4)


Если ус заменить на ус = ць + z1_a ab


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.



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

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

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


(Уд-Уь)-2,


{b3


Ь +


°о 5 Ла


-Ль5 (Уд-Уь) + 2|


(1 -а/2)


2ь +01'


(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 +    ) •    (10)

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

Г0 = (Уд " Уь) + Z(1-a) УУЬ + Уд •    О11)

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


5


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

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

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

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

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

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

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

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

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

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

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

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