Большинство случайных величин, с которыми
сталкиваются на практике, являются непрерывными. Для того, чтобы математически
описать случайную величину, необходимо указать множество ее значений и
соответствующее случайно величине распределение вероятностей для этого
множества (таблично, аналитически или графически).
2.1 Оценивание генеральных параметров
Математическое ожидание
Оценками математического ожидания μобычно
служит среднее арифметическое
, а иногда и
среднее геометрическое хr, медиана (серединное значение)
или
мода. Если провели n измерений случайной величины Х: х1, х2, … хn, то их
среднее арифметическое или выборочное среднее можно найти по формуле:
Если результат изменения выдается в виде
логарифма искомой величины, то конкретной оценкой математического ожидания
является среднее геометрическое xr:
Значение xrв сегда немного меньше
,
найденного по той же выборке. Если случайная погрешность измерений не столь
велика (менее 20%), то различие между xrи
несущественно,
поэтому на практике обычно для характеристики μ пользуются
средним арифметическим.
Чтобы найти медиану, результаты измерения
располагают в ряд (ранжированный ряд) в порядке возрастания их значений
Если
n- нечетное число, то медиана
равна среднему
члену ряда. Если n - четное число, то значение
равно
полусумме двух средних членов ряда. Медиана в отличие от среднего
арифметического нечувствительна к выделяющимся измерениям, поэтому при низком
качестве измерений и малом числе n она является лучшей оценкой математического
ожидания.
Мода - наиболее часто встречающееся значение величины при большой серии измерений.
Генеральная дисперсия
Оценка генеральной дисперсии σ2
служит выборочная дисперсия, которую рассчитывают по формуле:
Величину f=n-1 называют числом степеней свободы выборочной характеристики и определяют как число измерений, используемое для расчета данной характеристики, за вычетом числа связей, которые наложили на выборку при вычислении этой характеристики. Значение fуказывает на надежность определения генерального параметра по данной выборке.
Если известно математическое ожидание данной
совокупности измерений, то дисперсию S2вычисляют по формуле:
При использовании данной формулы на выборку не наложили ни одной связи, поэтому f=n.
Положительное значение корня квадратного из
дисперсии называют стандартным отклонением:
Относительное стандартное отклонение (коэффициент вариации V) равно:
в частях
;
в процентах
.
2.2 Распределение случайных величин
Основной характеристикой случайной величины является закон ее распределения, который можно задать в дифференцированной и интегрированной формах.
Если случайная величина Х принимает значения х1, х2, … хn и нам известна вероятностьР1, Р2, … Рn появления каждого из этих значений, то графическое или аналитическое представление зависимости P=φ(х)называется дифференцированным законом распределения или плотностью вероятности случайной величины.
Интегрированный закон распределения F(x) случайной величины указывает вероятность события, когда случайные величина Х принимает i-е значения, для которых xi<xa, где ха - некоторое определенное значение этой случайной величины. Интегральная функция представляет монотонно возрастающую кривую, которая стремится к единице (или к 100%). Ордината кривой соответствует точке ха, определяет вероятность события, когда случайная величина Х окажется меньше ха.разность двух ординат, соответствующих точкам ха и хb, определяет вероятность того, что случайная величина Х принимает значения из интервала ха ÷ хb: P(ха >xi> хb)=F(хb)-F(ха).
Функция распределения случайной величины является некоторой абстрактной математической моделью, при помощи которой описываются экспериментально наблюдаемые величины. Одна из задач статистической обработки цифрового материала заключается в нахождении этой функции, чтобы в дальнейшем использовать ее для статистической обработки данных.
Ниже рассмотрим некоторые распределение случайных величин, которые чаще всего применяются в аналитической химии.
Нормальное распределение случайной величины
Наиболее применяемым в практике измерения является нормальный закон распределения Гаусса. Этот закон является основой классической теории погрешностей измерений и современной статической обработки результатов эксперимента. Кроме того, это предельный закон для некоторых других законов распределения.
Функция плотности вероятностей нормального
закона распределения определяется формулой:
Табулировать его относительно хiне
представляется возможным, т.к. параметры
и
зависят
от абсолютных значений измеренной случайной величины, поэтому в практике обычно
пользуются нормальным законом распределения, где вместо переменной xi
используется переменная ui:
Рис. 1 График
дифференциальной кривой нормированного закона распределения
В этом случае независимо от конкретного значения
случайной величины математическое ожидание u=0, а
(u)=1.
Функция плотности вероятности приобретает вид:
Для кривой
центр
симметрии совпадает с началом координат, а на ось абсцисс наносится величина
ui- отклонение от среднего выражения в долях
.
Пользуясь табулированными значениями функции
можно подсчитать, что в области
(u=1) лежит 68,2%
всех результатов измерения, то есть около 16% результатов измерения отклоняются
от среднего менее чем
, и около 16% -
больше чем
. Вероятность
появления результатов за пределами
равно
соответственно 5% и 0,3% (Рис 1).
Распределение Пуассона
Результаты измерения целочисленных величин часто
подчиняются закону распределения Пуассона, например, подсчет числа квантов в
радиохимии или структурных единиц в минералогическом анализе. Функция плотностей
вероятности в этом случае имеет вид:

Рис 2. График
плотности вероятности распределения Пуассона
Из формулы видно, что плотность вероятностей
закона Пуассона характеризуется одним параметром
.
Генеральная дисперсия связана с
соотношением
Для малых значений
функция
существенно
ассиметрична, но асимметрия снижается с ростом
,
и форма кривой распределения приближается к нормальному сопротивлению со
значением генеральных параметров
и
.
Для практических целей при
можно считать, что
результаты измерения подчиняются нормальному закону распределения.
Распределение Стьюдента
Рис. 3 График
плотности вероятности распределения Стьюдента
Нормальный закон распределения случайных величин показывает, что вероятность появления малых отклонений от среднего значительно больше вероятности появления больших отклонений. Поэтому естественно ожидать того, что если исследователь выполнил небольшое число измерений одной величины, то среди них результатов с большим отклонением от среднего не будет. Если по этой выборке оценить дисперсию S2, то ее значение будет меньше соответствующей генеральной дисперсии. В силу этих обстоятельств классическая теория ошибок, основанная на нормальном законе распределения, не применима для обработки малого числа измерений. Статистика малых выборок стала возможна с появлением распределения Сттьюдента (t-распределение).
Если обработке подвергается небольшое количество
измерений, то вместо генеральной дисперсии
вынуждены
пользоваться выборочной дисперсией S2, близость значения которой к генеральной
дисперсии
зависит от числа
степеней свободы f, по которой рассчитывали S2. В этом случае переменную
заменяют
на новую переменную
, не содержащую
неизвестной величины
Функция
потности
вероятностей t-распределения симметрична относительно
,
то есть максимум t-распределения и нормального распределения соответствует
одному значению абсциссы. Однако в случае с t-распределением высота и ширина
кривой зависит не только от S, но и от f, по которым рассчитывается S.
В Таблице 1 указаны значения t в зависимости от
степеней свободыf.
Таблица 1. Коэффициенты t распределения Стьюдента
|
Число степеней свободы f |
Уровень значимости α |
|||
|
|
0,10 |
0,05 |
0,02 |
0,01 |
|
1 |
6,31 |
12,71 |
31,82 |
63,66 |
|
2 |
2,92 |
4,30 |
6,97 |
9,93 |
|
3 |
2,35 |
3,18 |
4,54 |
5,84 |
|
4 |
2,13 |
2,78 |
3,75 |
4,60 |
|
5 |
2,02 |
2,57 |
3,37 |
4,03 |
|
6 |
1,94 |
2,45 |
3,14 |
3,71 |
|
7 |
1,90 |
2,37 |
3,00 |
3,50 |
|
8 |
1,86 |
2,31 |
2,90 |
3,36 |
|
9 |
1,83 |
2,26 |
2,82 |
3,25 |
|
10 |
1,81 |
2,23 |
2,76 |
3,17 |
|
20 |
1,73 |
2,09 |
2,53 |
2,75 |
|
|
1,66 |
1,96 |
2,33 |
2,58 |
Рис. 4 График плотности вероятности распределения Фишера
Распределение Фишера
Отношение двух выборочных дисперсий
,
пренадлезащих различным совокупностям, подчиняется распределению Фишера или
F-распределению. Функция
плотности
вероятностей определяется сложным выражением, зависящим от значений Fи степеней
свободы f1 иf2, по которым подсчитываются выборочные дисперсии. Величина F
изменяется в пределах
. Кривая плотности
вероятностей ассиметрична, поэтому
,
но асимметрия уменьшается с увеличением значений f1 и f2.
Таблица 2. Значение F для уровня значимости α=0,10
|
f2 |
f1 |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
1 |
39,86 |
49,50 |
53,83 |
55,83 |
57,24 |
58,20 |
58,91 |
59,44 |
59,86 |
60,20 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
2 |
8,36 |
9,00 |
9,16 |
9,24 |
9,29 |
9,33 |
9,35 |
9,37 |
9,38 |
9,39 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
3 |
5,54 |
5,46 |
5,39 |
5,34 |
5,31 |
5,28 |
5,27 |
5,25 |
5,24 |
5,23 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
4 |
4,54 |
4,32 |
4,19 |
4,11 |
4,05 |
4,01 |
3,98 |
3,95 |
3,94 |
3,92 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
5 |
4,06 |
3,78 |
3,61 |
3,52 |
3,45 |
3,40 |
3,37 |
3,34 |
3,32 |
3,30 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
6 |
3,78 |
3,46 |
3,29 |
3,18 |
3,11 |
3,05 |
3,01 |
2,98 |
2,96 |
2,94 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
7 |
3,59 |
3,26 |
3,07 |
2,96 |
2,88 |
2,83 |
2,78 |
2,75 |
2,72 |
2,70 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
8 |
3,46 |
3,11 |
2,92 |
2,81 |
2,73 |
2,67 |
2,62 |
2,59 |
2,56 |
2,54 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
9 |
3,36 |
3,01 |
2,81 |
2,69 |
2,61 |
2,55 |
2,51 |
2,47 |
2,44 |
2,42 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
10 |
3,28 |
2,92 |
2,73 |
2,61 |
2,52 |
2,46 |
2,41 |
2,38 |
2,35 |
2,32
Таблица 4. Значение F для уровня значимости α=0,025
Χ2- распределение
Рис.
5 График плотности вероятности Χ2-
распределение
Пусть имеем выборочную совокупность х1, х2, …
хn, которая подчиняется нормальному закону распределения с генеральными
параметрами μ и σ2.
Если в нормальном законе распределения переменную Таблица 5. Значения
r-распределение
Ели в переменной
то получим распределение относительно отклонения или r-распределение. Функция φ(r) плотности вероятностей зависит от значения riи числа степеней свободы f=n-2. При больших значениях nr-распределение переходит в нормальное распределение. Таблица 6. Значение r для различных уровней значимости
Глава 3. Статистическое оценивание результатов
измерений
3.1 Методы исключения грубых ошибок
При проведении любого измерения встречаются результаты, которые резко отличаются от остальных. Вполне естественно, что в подобных случаях возникают подозрения о наличии грубых ошибок (промахов), допущенных исследователем при получении результатов. Единственным надежным методом выявления грубых ошибок является детальный анализ условий эксперимента. Если было установлено, что были нарушены стандартные условия, то сомнительный результат следует отбросить вне зависимости от его величины. Если в ходе проведения эксперимента был получен результат, который резко отличается от предыдущих, и условия оставались прежними, то этот сомнительный результат надо исключить, но заменить тремя новыми, подряд проведенными измерениями. На практике далеко не всегда удается
воспользоваться приведенными выше рекомендациями для исключения выбросов,
поэтому для их обнаружения обычно обращаются статистическим критериям. Стоит
отметить, что данные критерии являются условными, поскольку они базируются на
нормальном законе распределения случайных величин, допускающем изменение xот Метод исключения выбросов при известном σ Если оценка S стандартного отклонения σ
надежно
установлена (оценка при
где V - коэффициент вариации. Метод исключения выбросов при неизвестном σ Пусть имеется n изменений x1, x2, x3, … xb … xnи
среди них значение xi вызывает сомнение, так как существенно отличается от
остальных значений. Выдвигается нуль-гипотеза, которая заключается в том, что
сомнительное значение xi принадлежит данной совокупности, то есть не является
выбросом. Для принятия или отклонения ее можно используют r-критерий, значение
которого рассчитывается по формуле:
где Полученное значение ri подчиняется
r-распределению с числом степеней свободы f=n-2; его сопоставляют с табличным
значением rmax(min), значения которых приведены в таблице 7. Если рассчитанное
значение Таблица 7. Значения
Можно дать следующие рекомендации для выявления выбросов: . Если среди n измерений имеются два результата,
из которых один выказывает сомнения в силу того, что он значительно больше
остальных, а другой из-за того, что он значительно меньше, то следует по всем n
результатам рассчитать . Если среди измерений имеются два результата,
которые высказывают сомнения, потому что они значительно меньше (или больше)
остальных, то сначала произвольно отбрасывают худший из значительных
результатов. По оставшимся (n-1) измерениям подсчитывают 2.2 Сравнение двух дисперсий
В любой экспериментальной работе, включая решение задач по аналитической химии, нередко требуется сравнить точность результата измерений. Если мы имеем две серии измерений: x1, x2, x3 …,
xb … xn1иy1, y2, y3 … yb… yn2. По результатам измерений рассчитаны выборочные
дисперсии Принятие или подтверждение выдвинутой
нуль-гипотезы проводят с помощью F-критерия. Для этого находят отношение В первом случае заранее известно, что
Если Во втором случае сравнивают дисперсии 2.3 Сравнение нескольких дисперсий
При создании стандартных образцов, проверке
качества работа нескольких однотипных приборов и т.п. возникает задача оценить
однородность нескольких (m) дисперсий При сравнении дисперсий с помощью критерия
Бартлетта рассчитывают величину В:
где Бартлетт показал, что величина В распределена
приближённо как критерий χ2 с числом
степеней свободы Если Если число степеней свободы для всех выборок
одинаково, то есть
Для установления однородности/неоднородности
дисперсий используют табличные значения 2.3 Сравнение двух средних результатов
В экспериментальной работе часто возникает ситуация, когда требуется сравнить между собой два результата измерения: сопоставление результатов одной пробы, полученные разными методиками; сравниваются измерения, сделанные разными приборами и т.д. С позиции статистически данные задачи решаются одинаково. Пусть в разных условиях получены две выборки: х1, х2, х3, …, хn1 (проведено n1измерений), у1, у2, у3, …, уn2 (проведено n2измерений), и требуется сравнить их средние результаты Однако прежде чем оценивать значимость различия
между х и у, следует проверить однородность дисперсий Сравнение равноточных результатов При
Рассчитанное значение t сравнивают с табличным
значением t-критерием, установленным для уровня значимости α
и
числа степеней свободы
Сравнение неравноточных результатов Если при оценке однородности дисперсий
Величину fв общем случае рассчитывают по
формуле:
Следует отметить, что в зависимости от того,
насколько различаются по величине дисперсии После определения величины f находят табличное
значение t(α, f) для принятия
или опровержения нуль-гипотезы применяют алгоритм, который был описан выше для
равноточных результатов.
Заключение
Исходя из вышеперечисленного, можно сделать вывод, что методы математической статистики обработки результатов находят широкое применение на практике. Они не только позволяют оценить вклад случайной погрешности, но и помогают определить промахи (грубые погрешности) и исключить их без проведения дополнительных измерений; сравнить две или несколько дисперсий, что, например, применимо в случаях, когда нужно сравнить воспроизводилось двух методик анализа или при установке нового оборудования возник вопрос: не изменилось ли воспроизводилось результатов измерения. Сравнение нескольких дисперсий находит свое применение при изготовлении стандартных образцов, проверке качества однотипных приборов, установленных в лаборатории и т.д. Таки образом, математическая статистика является
одним из ключевых компонентов при проведении анализа.
Список литературы
1. ГармашA.B. Метрологические основы аналитической химии / А.В. Гармаш, А.М.Сорокина - Москва 2012 г. . Дёрффель К. Статистика в аналитической химии / К. Дёрффель - М.: Мир 1994 г. . Классификация погрешностей измерений [Электронный ресурс] Режим доступа:http://metrob.ru/HTML/pogreshnost/klassifikacia-pogreshnosti.htmlДата доступа: 27.11.2015 погрешность химия случайный дисперсия . Краткие сведения о статической обработке эксперементальных данных [Электронный ресурс] Режим доступа:http://alnam.ru/book_a_chem2.php?id=26Дата доступа: 1.12.2015 . Сергеев А.Г. Метрология, стандартизация, сертефикация / А.Г. Сергеев, М.В. Латышев, В.В. Терегеря - Москва 2005 г. . Смагунова А.Н. Методы математической статистики в аналитической химии / А.Н. Смагунова, О.Р. Карпукова - Феникс, 2012. . Теория вероятностей и статистика в аналитической химии [Электронный ресурс] Режим доступа:http://chemstat.com.ru/lections Дата доступа: 1.12.2015 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||