місяць;
день;
рік;
ідентифікація оператора;
ідентифікація контрольного еталона;
ідентифікація міри;
д)середнє J повторень;
стандартне відхилення повторюваності від J повторень;
ступені свободи;
умови навколишнього середовища (за необхідності).
Рисунок 2 — Дворівневий ієрархічний експеримент
Для цього дворівневого ієрархічного експерименту таблицю ANOVA, таблицю 2, можна отримати з випадку з трирівневим експериментом.
Таблиця 2 — Таблиця ANOVA для дворівневого ієрархічного експерименту
Джерело |
Ступені свободи, V |
Сума квадратів, SS |
Середній квадрат, MS |
Очікуваний середній квадрат |
День |
К- 1 |
ssD |
MSD |
а2 + Уод |
Похибка |
K(J-1) |
SSe |
MSE |
а2 |
Стандартне відхилення похибки з К (J- 1) ступенями свободи обчислюють так:
, л к J _ 2
Де
Середній квадрат для складової дня, MSD, з К-1 ступенями свободи дорівнює:
К-1
де стандартне відхилення, що характеризує щоденну варіацію, дорівнює:
- _ MSD-MSE
gd ~ ~ і
якщо різниця під знаком квадратного кореня є додатною. Інакше, ”aD дорівнює нулю.
Наслідком використання класичної оцінки є те, що це спричиняє появу від'ємних оцінок дисперсії. Інші оцінки можуть не мати цієї властивості та можуть використовуватися, якщо це є доречним.
Конфігурація вимірювання
Інші джерела невизначеності
Вимірювання на об’єктах випробування зазвичай проводять в один день, одним оператором, на одному ЗВТ тощо. Якщо невизначеність використовують для характеризування всіх вимірювань, що проводять у лабораторії, вона має передбачати будь-які розбіжності через:
— ЗВТ,
операторів,
геометрію,
інше.
Складову неконтрольованих навколишніх умов в лабораторії часто можна оцінити з даних контрольного еталона за період часу. Методи обчислення відповідних складових невизначеності обговорюють в іншому місці цього стандарту. Невизначеності, що є результатом контрольованих факторів, таких як оператори чи вибрані для певних вимірювань ЗВТ, описано в цьому підпункті.
Треба враховувати, що оператори мають досліджуватися лише раз, або за видами експериментів, що залежать від часу, або за конфігурацією вимірювання. Прикладами причин для розбіжностей в добре обслуговуваній лабораторії є:
розбіжності між ЗВТ для вимірювань похідних одиниць, наприклад, поверхневий опір кремнію, де ЗВТ не можна безпосередньо відкалібрувати до вихідного еталона;
розбіжності між операторами для оптичних вимірювань, що не є автоматичними та суттєво залежать від відеорозпізнавання оператора;
розбіжності між геометричними чи електричними конфігураціями апаратури.
Відкалібровані ЗВТ зазвичай не відносять до цього класу, оскільки невизначеності, що пов'язані з калібруванням, часто визначають оцінюванням за типом В, а ЗВТ у лабораторії мають узгоджуватися у межах невизначеностей калібрування. Для ЗВТ, вхідні сигнали яких безпосередньо не калібрують до визначеної одиниці, застосовують оцінювання за типом А. При цьому розглядають ситуації, коли вимірювання проводять за процедурою випробовування або стандартною методикою з використанням певного типу ЗВТ.
Проте треба зауважити, що деякі систематичні складові не можна усунути калібруванням, наприклад, матричні складові в аналітичній хімії.
Важливість змісту невизначеності
Розбіжності, про які йшлося на початку 5.3.1, вважають або випадковими розбіжностями, або розбіжностями зміщення. Це визначають, головним чином, вимоги до заявленої невизначеності. Наприклад, якщо проблемою є складова ЗВТ, один підхід може розглядати, наприклад, ЗВТ у лабораторії як випадкову вибірку ЗВТ одного типу та оцінювати невизначеність щодо всіх результатів незалежно від певного ЗВТ, яким проводять вимірювання. У цьому випадку можна застосовувати дворівневий ієрархічний експеримент з 5.2.3, в якому другий рівень призначено для одного з джерел впливу, таких як джерело ЗВТ. Другий підхід — оцінювання невизначеності щодо результатів за допомогою конкретного ЗВТ, який розглядають як аналіз систематичної складової чи зміщення в 5.5.
Нижче наведено простий підхід з використанням дворівневого ієрархічного експерименту з випадковими складовими для оцінювання невизначеності для одного джерела впливу.
Збирання даних і обчислення складової дисперсії
Щоб оцінити невизначеність процесу вимірювання, обумовлену ЗВТ, слід зробити випадкову вибірку /(/> 1) ЗВТ з тих, що є в наявності. Необхідно провести вимірювання на Q(Q > 1) артефактах кожним ЗВТ. За умови наявності / • Q вимірювань, стандартне відхилення, що описує розбіжності між ЗВТ, обчислюють наведеним нижче чином із середнього значення для кожного ЗВТ та має І - 1 ступені свободи:
е _ < /=1
у І _ '
де для /-го ЗВТ:
1 Q лі
Уі.=-ЎУ:а, У..=-УУі.,
Iq
Приклад аналізу випадкових розбіжностей
У таблиці 3 наведено дані вимірювання питомого опору (Ом • см) з п’ятьма зондами (позначеними номерами 1, 281, 283, 2 062, 2 362) на Q = 5 пластинах (позначених номерами 138, 139, 140, 141, 142). Ці самі дані аналізують для зміщення значень для зонда № 2362 в 5.5 цього стандарту. Наведено середнє значення для кожного зонда на артефактах. Стандартне відхилення (середніх значень вимірювань питомого опору) для зондів складає 0,0219 Ом ■ см із чотирма ступенями свободи. Таким чином, S/nst= 0,0219.
Таблиця 3 — Вимірювання питомого опору п’яти пластин з використанням п’яти зондів
Значення в Ом ■ см
№ зонда |
Ідентифікаційний номер пластини |
Середнє значення |
||||
138 |
139 |
140 |
141 |
142 |
||
1 |
95,1548 |
99,3118 |
96,1018 |
101,1248 |
94,2593 |
97,1905 |
281 |
95,1408 |
99,3548 |
96,0805 |
101,0747 |
94,2907 |
97,1883 |
283 |
95,1493 |
99,3211 |
96,0417 |
101,1100 |
94,2487 |
97,1742 |
2 062 |
95,1125 |
99,2831 |
96,0492 |
101,0574 |
94,2520 |
97,1508 |
2 362 |
95,0928 |
99,3060 |
96,0357 |
101,0602 |
94,2148 |
97,1419 |
Для графічного аналізу розбіжності між виміряними значеннями та середнім значенням для кожного зонда можна зобразити залежно від пластини, для кожного зразка, до того ж зонди буде окремо позначено певним графічним символом. Графік розглядають, щоб визначити, чи завжди деякі ЗВТ показують верхні чи нижні значення відносно інших ЗВТ, та чи узгоджується ця характеристика із зондами. З графіка на рисунку 3, створеного на основі даних з таблиці 4 (див. 5.5.2.2), видно, що між зондами є розбіжності, наприклад, зонди № 2062 і № 2362 однаково показують нижні значення відносно інших зондів.
Неоднорідність матеріалів
Проблеми, спричинені неоднорідностями
Артефакти, електричні прилади, хімічні речовини тощо можуть бути неоднорідними відносно кількості, яку характеризують. Неоднорідність може стати фактором під час оцінювання невизначеності, коли:
артефакт характеризується одним значенням, проте є неоднорідним на поверхні тощо;
багатьом об’єктам задають однакове значення на підставі кількох вибірок з партії, а партія є неоднорідною від вибірки до вибірки.
Невдалим аспектом цієї ситуації є те, що неоднорідність може бути джерелом невизначеності, що переважає. Навіть якщо сам процес вимірювання є дуже точним і його статистично контролюють, сумарна невизначеність все ще може бути неприйнятною для практичних цілей через неоднорідності матеріалів. Дослідження однорідності для стандартних зразків детально обговорено в Guide ISO 35 [14].
Дії для випадкових неоднорідностей
Випадкові неоднорідності оцінюють за допомогою статистичних методів для кількісного визначення випадкових складових. Прикладом неоднорідності є хімічний стандартний зразок, що не може бути достатньо гомогенізований відносно ізотопів, що викликають інтерес. Ізотопні відношення може бути визначено вимірюваннями на кількох пробах, відібраних навмання з партії.
Збирання даних і обчислення складової для неоднорідності
Простою схемою для ідентифікування й кількісного визначення впливу неоднорідності на результат вимірювання є збалансований дворівневий ієрархічний експеримент. З партії навмання вибирають К(К > 1) об’єктів випробування, проводять J(J > 1) вимірювань на об’єкті випробування. Вимірювання позначають як:
Vl, Ук2’-- < ^KJ’
де індекс k - 1, ..., К означає об’єкти випробування для Рівня 2 та; = 1,...., J означає повторення для об’єкта випробування для Рівня 1.
Дисперсію неоднорідності (між об'єктами випробування), що визначають як складову дисперсії, обумовлену неоднорідністю об'єктів випробування, обчислюють як у 5.2.3 із використанням методики ANOVA, де:
2Msitem-MSE
"^inh і
d К 2 і К J 2
=—У(Ук.-У..) т—гУУ(Yk! -Yk.} ,
Де - 1 J
У«-=7Ж:
J /-і
MSitem— середній квадрат, обумовлений об’єктами випробування;
S'Y — оцінка складової дисперсії, обумовленого неоднорідністю матеріалу чи об’єктами випробування.
Якщо S?h— від’ємне, то вплив неоднорідності статистично вважають рівним нулю, і внеску у невизначеність немає. Тобто невизначеність, пов’язану з неоднорідністю, описують як:
uinh= max (Sinh. 0).
Оцінювання невизначеності, пов’язаної з неоднорідністю
Оцінювання невизначеності залежить від використання результату вимірювання. Зазвичай, неоднорідність є важливою, коли отримується середнє значення для низки об'єктів випробування з великої партії, і це середнє значення приписується кожному об’єкту випробування в партії. Для результату вимірювання, що обчислюють як середнє з результатів з К різних об'єктів випробування, стандартну невизначеність щпЛ, що виникає внаслідок неоднорідності та пов’язана з середнім результатом, обчислюють за 8/лЛ відповідно до:
Проте для результату вимірювання, що обчислюють як середнє з результатів з К різних об’єктів випробування та застосовують до кожного об’єкта в решті партії, стандартна невизначеність uinh, що виникає внаслідок неоднорідності та пов’язана з предикційним інтервалом для кожного об'єкта з решти партії, ґрунтується на цьому інтервалі (див. [5]) та її задають так:
ґ І Г “inh=max J1 + —-S/n/1,0 І у Л J
Дії для систематичних неоднорідностей
Систематичні неоднорідності потребують дещо іншого підходу. Наприклад, шорсткість може змінюватися систематично на поверхні 50 мм квадратної металевої деталі, що має певний профіль шорсткості. Лабораторія з сертифікації може вимірювати цю деталь у кількох місцях, проте, за винятком, коли можливо охарактеризувати шорсткість як функцію координат на деталі, необхідно оцінювати неоднорідність як джерело невизначеності.
У цьому разі за однією дією треба обчислювати заявлене значення як середнє з вимірювань на поверхні деталі та оцінювати невизначеність для відхилень від середнього значення. Складові невизначеності можна оцінити одним із кількох методів для оцінювання невизначеності за типом В, наведених у GUM.5.5 Зміщення, зумовлене конфігураціями вимірювання
Загальні положення
У статистиці для того, щоб оцінити параметр 0, визначають зміщення оцінки 0 як різницю між математичним сподіванням 0 та істинним значенням 0. А саме b = Е0- 0. Оскільки істинне значення 0 невідоме, то b — невідоме. Коли існує оцінка Ь, її називають поправкою та позначають Ь. Поправка може бути по відношенню до опорного значення або до певного виду середнього значення. За заданого набору поправок, bv...,bn, зміщення оцінки можна визначити середнім із поправок, тобто:
п
л
Якщо поправки вважають випадковими, то можна прийняти нормальний розподіл для її оцінювання.
Якщо поправки зосереджуються біля нуля, часто припускають, що розподілення ймовірності для поправок має середнє значення, що дорівнює нулю, і цей випадок часто називають «нульовою» поправкою. У цьому випадку 0 є незміщеною оцінкою 0. Якщо поправки нормально розподілено, або існує велика низка поправок, то невизначеність зміщення можна оцінити стандартним відхиленням вибіркового середнього з поправок. Якщо про розподіл існує недостатньо інформації, то можна припустити, що поправки {Ь(, і = 1, ..., п} рівномірно розподілено між мінус а та а. Зміщення, таким чином, оцінюють за нулем. Величину а можна оцінити так: