Umk_TV_2009__MATAN_2
.pdf
Полигон относительных частот этого распределения изображен на рис.3.2.
Рис. 3.2
Для полученного статистического ряда вычислим значения эмпирической функции распределения, используя формулу (3.3)
|
|
0, |
при x -4, |
|
|
|
|
1 |
|
|
|
|
, при - 4 x -2, |
|
|
|
|
||
|
|
|
20 |
|
|
|
|
4 |
, при - 2 x -1, |
|
|
|
|
|
|
|
20 |
||
|
|
|
|
|
F * |
x |
|
9 |
, при -1 x 0, |
|
|
|||
20 |
|
20 |
|
|
|
|
|
|
|
|
|
|
15 |
, при 0 x 2, |
|
|
|
|
|
|
|
20 |
|
|
|
|
|
|
|
|
|
|
18 |
, при 2 x 3, |
|
|
|
20 |
|
|
|
|
|
|
|
|
1, при x 3. |
||
График этой функции представлен на рис. 3.3.
71
Рис. 3.3
Гистограмма
Если наблюдаемая случайная величина ξ непрерывна или объем выборки n большой, то вариационный и статистический ряды будут трудно обозримыми множествами, практически не будет равных элементов выборки. В этом случае
используется процедура группировки выборки, которую рассмотрим для
|
x1, , xn . Интервал возможных значений |
|
реализации выборки x |
ξ делят |
точками z0 z1 zk на k непересекающихся полуинтервалов (разрядов)
i |
zi 1, zi , |
i |
|
. Для каждого разряда i подсчитывают частоту mi |
|
|
1, k |
- |
|||||
|
|
|
|
k |
|
|
число элементов выборки, попавших в этот разряд. При этом mi |
n . |
В |
||||
|
|
|
|
i 1 |
|
|
интервал включают значения, большие или равные нижней границе и меньше верхней границы. Далее находят относительные частоты (статистические
вероятности) pi* mni . Группированные данные удобно представить в виде
интервального статистического ряда – последовательности пар i , pi* , или в
виде таблицы (табл. 3.3). Часто группу элементов выборки, входящих в интервал
72
i , заменяют средней точкой zi zi 12 zi .
Таблица 3.3
|
i |
|
1 |
2 |
… |
k |
|
|
|
||
|
p |
* |
|
|
* |
* |
… |
* |
|
|
|
|
|
i |
|
p1 |
p2 |
|
pk |
|
|
|
|
|
Обычно длина разрядов выбирается одинаковой, т.е. равной h |
zk z0 |
. |
||||||||
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
k |
|
Число разрядов k |
выбирается в зависимости от объема выборки n так, чтобы |
||||||||||
построенный ряд не был громоздким и в то же время позволял выявить
характерные особенности изменения случайной величины. |
Для определения k |
можно рекомендовать формулу Стерджеса |
|
k 1 1,44ln n, |
(3.4) |
которая дает нижнюю оценку величины k . В качестве значения k следует брать ближайшее целое число.
Группированный статистический ряд наглядно можно изобразить в виде гистограммы. Для ее построения на оси абсцисс откладывают разряды i
длиною h , и на каждом из них, как на основании, строят прямоугольник. В
результате получают ступенчатую фигуру, которую называют гистограммой.
Высота i–го частичного прямоугольника при построении гистограммы |
|
||
частот равна отношению mi h (плотность частоты). |
|
|
|
Площадь i–го частичного прямоугольника численно равна |
mi |
h m |
, |
|
|||
|
h |
i |
|
|
|
|
|
а площадь гистограммы частот численно равна объему выборки, т.е.
73
k
mi n . (3.5) i 1
При построении гистограммы относительных частот:
высота i–го частичного прямоугольника равна отношению относительной
частоты к длине интервала |
mi / n |
(плотность относительной частоты); площадь |
||||||||
h |
||||||||||
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
h |
n |
h |
n |
||
i–го частичного прямоугольника |
численно равна |
mi |
|
mi |
; площадь |
|||||
|
|
|
|
|||||||
гистограммы относительных частот численно равна |
|
|
|
|
|
|||||
|
k m |
|
|
|
|
|
||||
|
|
|
i |
1. |
|
|
|
(3.6) |
||
|
|
|
|
|
|
|||||
|
i 1 n |
|
|
|
|
|
||||
Гистограмма относительных частот является статистическим аналогом плотности распределения наблюдаемой случайной величины ξ.
Гистограмма изображена на рис. 3.4.
Рис. 3.4
74
3.3. Точечные оценки параметров распределения
На практике часто удается предсказать или оценить с помощью гистограммы вид распределения наблюдаемой случайной величины ξ с точностью до неизвестного параметра a (или нескольких параметров). Одной из основных
задач математической статистики является нахождение оценки (приближенного значения) неизвестного параметра по имеющейся выборке.
|
Основные понятия |
|
||
Пусть |
наблюдается случайная величина ξ с функцией |
распределения |
||
F x и плотностью распределения |
f x . |
Случайная выборка |
представлена |
|
|
X1, , X n с реализацией |
|
|
|
вектором X |
x x1, , xn . |
(3.7) |
||
Параметром распределения a случайной величины называется любая
числовая характеристика этой случайной величины (математическое ожидание,
дисперсия и т.п.) или любая константа, явно входящая в выражение для функции
или плотности распределения.
Если параметр a неизвестен, то его точечной оценкой называется произвольная функция элементов выборки
|
. |
|
an * g( X1, X 2 , , X n ) g X |
(3.8) |
Реализацию оценки, т.е. значение оценки для наблюдавшейся в эксперименте
реализации выборки, принимают за приближенное значение неизвестного параметра a an* (x1, x2 , , xn ).
Из соотношения (3.8) видно, что an * как функция случайных величин
сама также является случайной величиной. Закон распределения оценки an*
зависит от вида функции g , числа наблюдений и значения оцениваемого параметра.
75
Ясно, что существует много разных способов построения точечной оценки,
ине всякая зависимость g X может давать удовлетворительную оценку
неизвестного параметра a . Рассмотрим некоторые свойства, которыми должна обладать оценка, чтобы ее можно было считать хорошим приближением к неизвестному параметру.
Оценка |
an* |
параметра |
|
a называется |
несмещенной, |
если ее |
||
математическое ожидание равно оцениваемому параметру, то есть |
|
|||||||
|
|
|
M |
a |
* a . |
|
|
(3.9) |
|
|
|
|
|
n |
|
|
|
Если свойство (2.2) не выполняется, то есть |
|
|
|
|||||
|
|
|
M a |
* a b 0 |
, |
|
(3.10) |
|
|
|
|
|
n |
|
|
|
|
то оценку an* |
называют смещенной, при этом величину b называют |
|||||||
систематической ошибкой оценки an*. |
|
|
|
|||||
Требование несмещенности означает, что выборочные значения оценок,
полученных в результате повторения выборок, группируются около оцениваемого параметра.
Оценка |
a |
* |
параметра a называется состоятельной, если при n |
||||
|
|
n |
|
|
|
|
|
она сходится по вероятности к оцениваемому параметру a , т.е. |
для любого ε > 0 |
||||||
выполняется равенство |
|
||||||
|
|
|
lim P |
an * a |
|
1 . |
(3.11) |
|
|
|
|
||||
|
|
|
n |
|
|
|
|
Следующая теорема устанавливает достаточные условия состоятельности оценки an* параметра a .
76
Теорема. Если при n M an* a и D an* 0, то оценка an* параметра a является состоятельной.
Состоятельность оценки означает, что, при достаточно большом объеме
выборки с вероятностью близкой к единице, отклонение оценки от истинного значения параметра меньше ранее заданной величины.
Обычно в качестве меры точности оценки an* используется среднеквадратическая ошибка (среднее значение квадрата ошибки)
|
* |
2 |
|
. Очевидно, чем меньше эта ошибка, тем теснее сгруппированы |
M an |
a |
|
||
|
|
|
|
|
значения оценки около оцениваемого параметра. Поэтому всегда желательно,
чтобы ошибка оценки была по возможности малой. Используя свойства математического ожидания, нетрудно получить
|
2 |
|
D an* b2 . |
|
M an* a |
|
(3.12) |
||
|
|
|
|
|
Для несмещенных оценок
|
2 |
|
D an* , |
|
M an |
* a |
|
(3.13) |
|
|
|
|
|
|
то есть их мерой точности является дисперсия.
Несмещенная оценка параметра a называется его эффективной оценкой,
если ее |
дисперсия |
D an* |
является наименьшей |
среди дисперсий всех |
|
возможных оценок |
параметра |
a , вычисленных по одному и тому же объему |
|||
выборки. |
|
|
|
|
|
|
Точечные оценки математического ожидания и дисперсии |
||||
|
|
|
|
|
|
Пусть случайная |
выборка X X1, , X n |
порождена наблюдаемой |
|||
случайной |
величиной |
ξ, математическое ожидание |
M m и дисперсия |
||
77
D D 2 которой неизвестны. В качестве оценок для этих характеристик было предложено использовать выборочное среднее
|
|
|
|
1 |
n |
|
|
|
|
|
|
|||
m* X |
|
|
|
|
|
|
||||||||
|
|
X |
i |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
||||||||
n |
|
|
|
|
n i 1 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||
и выборочную дисперсию |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
n |
|
|
|
|
|
|
|
D* |
*2 |
|
X |
|
|
|
2 . |
|
||||||
i |
X |
(3.14) |
||||||||||||
|
||||||||||||||
n |
|
n |
|
|
n i 1 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
Рассмотрим некоторые свойства оценок математического ожидания и дисперсии.
1. Вычислим математическое ожидание выборочного среднего:
|
|
|
1 |
n |
|
|
|
|
|
||||
M X M |
|
X i |
|
|||
|
||||||
|
|
n i 1 |
|
|
||
1 |
n |
1 |
n |
|
M X i |
m m . (3.15) |
|||
|
|
|||
n i 1 |
n i 1 |
|||
Следовательно, выборочное среднее является несмещенной оценкой для m .
2. Напомним, что результаты наблюдений – независимые
случайные величины, каждая из которых имеет такой же закон распределения, как
|
, |
|
|
|
M X |
|
m, D X |
|
2 |
|
|
|
|
|||||||||
и величина |
а значит, |
|
i |
i |
, |
i 1,2 ,...n . Будем |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
предполагать, что дисперсия D конечна. Тогда, согласно теореме Чебышева о |
||||||||||||||||||||||
законе больших чисел, для любого ε > 0 имеет место равенство |
|
|
|
|||||||||||||||||||
|
|
|
|
1 |
n |
1 |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
lim |
P |
|
|
X i |
M X i |
|
|
1, |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|||||||||||||||
n |
|
|
n i 1 |
n i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
P |
|
|
m |
|
1. |
|
|
|
||||||||||||
которое можно записать так: |
lim |
X |
|
|
(3.16) |
|
||||||||||||||||
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
||||||
Сравнивая (3.16) с определением свойства состоятельности (3.11), видим, что оценка X является состоятельной оценкой математического ожидания M .
78
3. Найдем дисперсию выборочного среднего:
|
|
1 |
|
n |
|
1 |
|
|
n |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
||||||||||||||
D X D |
|
|
|
X i |
|
|
|
|
D X i |
|
|
|||||||||
|
|
|
|
|
|
|
||||||||||||||
|
|
n i 1 |
|
n2 i 1 |
|
|
|
|
||||||||||||
|
|
|
|
|
1 |
|
|
n |
|
|
1 |
|
|
|
2 |
|
||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
n 2 |
|
|
. |
(3.17) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
n2 i 1 |
|
n2 |
|
|
|
n |
|
||||||||
Таким образом, дисперсия оценки математического ожидания уменьшается |
||||||||||||||||||||
обратно пропорционально объему выборки. |
|
|
|
|
|
|||||||||||||||
Можно доказать, что если случайная величина ξ распределена нормально, |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
то выборочное среднее |
X |
является |
эффективной |
оценкой |
математического |
|||||||||||||||
ожидания M , то есть дисперсия |
D |
|
принимает наименьшее значение по |
|||||||||||||||||
X |
||||||||||||||||||||
сравнению с любой другой оценкой математического ожидания. Для других законов распределения ξ это может быть и не так.
Выборочная дисперсия *2 |
является смещенной оценкой дисперсии 2 |
, так как |
|||||
|
|
n |
|
|
|
|
|
M *2 |
|
n 1 |
2 |
2 . |
(3.18) |
|
|
|
|
||||||
n |
|
n |
|
|
|
|
|
|
|
|
|
|
|
||
Действительно, используя свойства математического ожидания и формулу (3.17),
найдем
|
1 |
n |
2 |
|
|
1 |
n |
|
|
|
|
|
M n*2 M |
|
|
|
|
|
|
||||||
|
X i X |
|
|
|
M X i m X m 2 |
|
|
|||||
|
|
|||||||||||
n i 1 |
|
|
|
|
n |
i 1 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
n |
|
|
|
|
|
n |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
M X i m 2 |
2 X m X i m X m 2 |
|
|
|||||||||||||
n |
||||||||||||||||||
|
|
|
i 1 |
|
|
|
|
|
i 1 |
|
i 1 |
|
|
|||||
|
|
|
1 |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 X m nX nm n X m 2 |
|
||||||||||||||
|
|
|
M |
X i m 2 |
|
|
||||||||||||
|
|
|||||||||||||||||
|
|
|
n |
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
||
79
|
1 |
n |
|
|
|
|
|
|
|
n X m 2 |
|
||||||
|
|
M X i m 2 |
|
|
||||
n |
||||||||
|
i 1 |
|
|
|
|
|
||
|
1 |
n |
|
|
|
|
|
|
|
|
|
|
|||
|
|
M X i m 2 |
nM X m 2 |
|
|
||
|
|||||||
|
n i 1 |
|
|
|
|
|
|
|
1 |
|
2 |
|
2 |
|
|
n 1 |
|
|
|
|
|
|
n |
|
|
|
|
|
|
2 . |
|
|
|
|
|
|
|||||||
|
n |
|
n |
|
n |
|
|||||
|
n |
|
|
|
|
|
|
||||
Чтобы получить несмещенную оценку дисперсии, оценку (3.14) нужно исправить,
то есть домножить на |
|
n |
. |
Тогда |
получим несмещенную |
выборочную |
|||||||
|
|
|
|||||||||||
|
n 1 |
||||||||||||
дисперсию |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
1 |
n |
|
|
|
|
|
||
S 2 |
D* |
X |
|
|
|
2 . |
|
||||||
i |
X |
(3.19) |
|||||||||||
|
|
||||||||||||
|
n 1 |
n |
|
n 1i 1 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
||||||
Отметим, что формулы (3.14) и (3.19) отличаются лишь знаменателем, и при больших значениях n выборочная и несмещенная дисперсии отличаются мало.
Однако при малом объеме выборки n 30 следует пользоваться соотношением
(3.19).
Для оценки среднего квадратического отклонения случайной величины используют так называемое “исправленное” среднее квадратическое отклонение,
которое равно квадратному корню из несмещенной дисперсии: S 
S 2 .
80
