Umk_TV_2009__MATAN_2
.pdf
3.4.Интервальные оценки
Встатистике имеются два подхода к оцениванию неизвестных параметров распределений: точечный и интервальный. В соответствии с точечным оцениванием, которое рассмотрено в предыдущем разделе, указывается лишь точка, около которой находится оцениваемый параметр. Желательно, однако,
знать, как далеко может отстоять в действительности этот параметр от возможных реализаций оценок в разных сериях наблюдений.
Ответ на этот вопрос – тоже приближенный – дает другой способ оценивания параметров – интервальный. В соответствии с этим способом оценивания находят интервал, который с вероятностью, близкой к единице,
накрывает неизвестное числовое значение параметра.
Понятие интервальной оценки
Точечная оценка a* a* ( X1, , X n ) является случайной величиной и для возможных реализаций выборки принимает значения лишь приближенно равные истинному значению параметра a . Чем меньше разность a a* , тем
точнее |
|
оценка. Таким образом, положительное число |
, |
для которого |
||||||||||
|
a a* |
|
, характеризует точность оценки и называется ошибкой оценки (или |
|||||||||||
|
|
|||||||||||||
предельной ошибкой). |
|
|
|
|
|
|
|
|
|
|
||||
|
Доверительной вероятностью (или надежностью) |
называется |
||||||||||||
|
|
|
|
|
|
|
|
, т. е. |
||||||
вероятность β, с которой осуществляется неравенство |
a a* |
|||||||||||||
|
|
|
|
|
P |
|
a a* |
|
. |
|
(3.20) |
|||
|
|
|
|
|
|
|
|
|||||||
|
a a* |
|
равносильным ему двойным неравенством |
|||||||||||
Заменив неравенство |
|
|
||||||||||||
a a* , или |
a* a a* , получим |
|
|
|||||||||||
81
P a* a a* . |
(3.21) |
Интервал I a* , a* , накрывающий |
с вероятностью β, |
0 1, неизвестный параметр a , называется доверительным интервалом
(или интервальной оценкой), соответствующим доверительной вероятности β.
Случайной величиной является не только оценка a* , но и ошибка
значение зависит от вероятности β и, как правило, от выборки. Поэтому
доверительный интервал случаен и выражение (3.21) следует читать так:
“Интервал I накроет параметр a с вероятностью β ”, а не так: “Параметр a
попадет в интервал I с вероятностью β ”.
Смысл доверительного интервала состоит в том, что при многократном повторении выборки объема n в относительной доле случаев, равной β,
доверительный интервал, соответствующий доверительной вероятности β,
накрывает истинное значение оцениваемого параметра. Таким образом,
доверительная вероятность β характеризует надежность доверительного оценивания: чем больше β, тем вероятнее, что реализация доверительного интервала содержит неизвестный параметр.
Следует, однако, иметь в виду, что с ростом доверительной вероятности β в
среднем растет длина доверительного интервала, то есть уменьшается точность доверительного оценивания. Выбор доверительной вероятности определяется
конкретными условиями; |
обычно используются значения β, равные 0,90; 0,95; |
|
0,99. |
|
|
Вероятность |
1 |
(3.22) |
82
называется уровнем значимости и характеризует относительное число ошибочных заключений в общем числе заключений.
В формуле (3.21) границы доверительного интервала симметричны относительно точечной оценки. Однако не всегда удается построить интервал,
обладающий таким свойством. Более общим является следующее определение.
Доверительным интервалом (или интервальной оценкой) параметра a с
доверительной вероятностью β, 0< β <1, называется интервал со случайными
границами a * a * X |
1 |
, , X |
n |
, a * a |
* X |
1 |
, , X |
n |
, накрывающий с |
|
1 |
1 |
|
2 |
2 |
|
|
||||
вероятностью β неизвестный параметр a , т. е. |
|
|
|
|
|
|||||
|
P a* |
a a* |
. |
|
|
|
|
(3.23) |
||
|
1 |
|
|
2 |
|
|
|
|
|
|
Иногда вместо двусторонних доверительных интервалов рассматривают односторонние доверительные интервалы, полагая a1* или a2* .
Построение интервальных оценок
Доверительный интервал задается своими концами a1* и a2* . Однако найти функции a1* a1* X1, , X n и a2* a2* X1, , X n из условия (3.23) невозможно, поскольку закон распределения этих функций зависит от закона распределения ξ и, следовательно, зависит от неизвестного параметра a .
Используют следующий прием, позволяющий в ряде случаев построить
доверительный интервал. Подбирается такая функция , чтобы:
G X , a
- ее закон распределения был известен и не зависел от неизвестного параметра a ;
- функция была непрерывной и строго монотонной по .
G X , a a
83
Тогда для любого β можно выбрать два числа g1 и g2 так, чтобы выполнялось равенство
g 2
P g1 G g2 F g2 F g1 fG t dt . (3.24)
g1
Отсюда находят g1 и g2 как квантили функции распределения .
G X , a
Границы искомого доверительного интервала выражают через найденные квантили и выборочные данные, используя для этого соотношения,
связывающие новую и старую случайные величины.
Если плотность распределения случайной величины
симметрична, то доверительный интервал симметричен относительно
G X , a
точечной оценки a* , и для нахождения границ доверительного интервала вместо условия (3.23) можно использовать соотношение (3.21).
Основные статистические распределения
Построение разного рода оценок и статистических критериев часто основывается на использовании ряда специальных распределений случайных величин.
Нормальное распределение. Случайная величина X имеет нормальное
распределение с параметрами m и 2 |
0, что обозначается как |
X ~ N m; 2 , если плотность вероятности этой случайной величины имеет
вид
|
|
1 |
|
|
|
2 |
|
|
f x |
|
|
|
x m |
|
|
||
|
|
|
exp |
|
|
. |
(3 .25) |
|
|
|
|
|
|
||||
2 |
|
2 |
||||||
|
|
|
2 |
|
|
|||
|
|
|
|
|
|
|
|
|
84
График плотности вероятности случайной величины, имеющей нормальное распределение, представлен на рисунке 3.5, на котором видно, что максимум функции находится в точке
Поскольку нормальное распределение подробно изучается в курсе теории вероятностей, напомним свойства нормальной случайной величины, которые будут использоваться в дальнейшем.
Рис. 3.5
1)M X m, D X 2 .
2)Случайная величина называется центрированной, если ее математическое
ожидание равно нулю. Для того чтобы центрировать случайную величину, надо
вычесть из нее математическое ожидание:
M X M X M X M X 0.
3)Случайная величина называется нормированной, если ее дисперсия равна единице, а математическое ожидание равно нулю.
85
Для того чтобы нормировать случайную величину, надо ее поделить на среднее квадратическое отклонение:
D |
X |
|
|
1 |
D X 2 1. |
|
|
||||
|
|
|
2 |
2 |
|
|
|
||||
Центрированная и нормированная нормальная случайная величина называется стандартной. Таким образом, стандартной будет случайная величина
|
|
|
|
|
|
Z |
X m |
~ N 0; |
1 . |
|
(3.26) |
|||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вероятность |
попадания |
случайной величины |
X в интервал |
(α,β) |
|||||||||
вычисляется по формуле |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
m |
|
m |
|
|
||
|
P X |
|
|
|
|
, |
(3.27) |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
z |
|
1 |
|
z |
t 2 |
|
|
|
|
|
|
|
|
где |
|
|
e 2 dt - интеграл |
вероятности, |
представляющий |
собой |
||||||||
|
|
|
||||||||||||
|
||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
функцию распределения стандартной нормально распределенной случайной величины. Интеграл вероятности табулирован. Его значения приведены в таблице В Приложения.
Для стандартной нормальной случайной величины и симметричного промежутка формула (3.27) принимает следующий вид:
P X P |
|
X |
|
2 1. |
(3.28) |
|
|
||||
|
|
|
|
|
|
86
Распределение 2 |
|
X i ~ N 0; |
|
i |
|
(хи-квадрат). Если |
, |
1, k |
|||
|
|
1 |
|
|
независимые стандартные нормальные случайные величины, то говорят, что случайная величина
k |
2 |
|
Zk X i |
(3.29) |
|
i 1 |
|
|
имеет распределение хи-квадрат с k степенями свободы, что обозначается как
Zk ~ 2 k . Графики плотности вероятности для двух значений степени свободы приведены на рис.3.6.
|
|
|
|
|
Рис. 3.6 |
|
|
|
|
С увеличением числа степеней свободы k |
плотность вероятности стремится к |
||||||||
нормальной. При k 2 плотность вероятности постоянно убывает, |
а при k 2 |
||||||||
имеет единственный максимум x k 2, M Zk k , |
D Zk 2k . |
||||||||
Распределение Стьюдента. |
Пусть X 0 , |
X1, , |
X k - |
независимые |
|||||
стандартные нормальные случайные величины. Тогда случайная величина |
|||||||||
|
|
|
|
|
X 0 |
|
|
(3.30) |
|
Tk |
|
|
|
|
|
|
|||
|
k |
|
|
|
|||||
1 |
|
|
|
|
|||||
|
X i 2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
||
|
|
|
k i 1 |
|
|
|
|||
87
имеет распределение Стьюдента с k степенями свободы, что обозначается как Tk ~ S k , при этом
|
|
M Tk 0 , |
D Tk |
k |
. |
|
|
|
|
|
|||
|
|
k 2 |
|
|||
|
|
|
|
|
|
|
На |
рис.3.7 |
приведены кривые |
стандартного |
нормального |
||
распределения |
(кривая 1) и плотности распределения Стьюдента (кривая 2). |
|||||
Рис. 3.7
При n плотность распределения Стьюдента стремится к плотности стандартной нормальной случайной величины.
На практике, как правило, используется не плотность вероятности, а
квантиль распределения. Напомним, что квантилью порядка (или уровня) p
непрерывной случайной величины X называется такое ее значение x p , которое удовлетворяет равенству F x p P X x p p ,
где F (x) - функция распределения, а p - заданное значение вероятности.
Рис.3.8 поясняет понятие квантили порядка p .
88
Рис. 3.8
Следующая теорема устанавливает свойства основных выборочных характеристик, вычисленных по выборке, соответствующих нормальному
распределению.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Теорема Фишера. Пусть |
X X1, , X n - |
случайная выборка из |
|||||||||||
генеральной |
совокупности ~ N m; 2 , тогда |
выборочное |
среднее |
|||||||||||||
|
|
|
1 |
n |
|
|
|
|
|
|
|
1 |
n |
|
|
|
|
|
|
X i |
и несмещенная выборочная дисперсия S 2 |
|
X i |
|
2 |
||||||||
|
X |
X |
||||||||||||||
|
|
|
||||||||||||||
|
|
|
n i 1 |
|
|
|
|
|
|
|
n 1i 1 |
|
|
|
||
независимы, и при этом |
|
|
|
|
|
|
|
|
|
|||||||
|
1) случайная величина |
|
m |
|
|
/ имеет распределение N 0; 1 |
|
|
|
|||||||
|
|
|
n |
|
|
|
||||||||||
|
X |
; |
|
|
||||||||||||
2)случайная величина n 1 S 2 / 2 имеет распределение 2 n 1 ;
3)случайная величина X m 
n / S имеет распределение S n 1 .
Доказательство теоремы приведено в [2].
89
Интервальные оценки математического ожидания нормального
распределения
Интервальная оценка математического ожидания при известной дисперсии.
Построим доверительный интервал для математического ожидания наблюдаемой |
||||||||||||||
случайной величины ~ N m; |
2 при известной дисперсии |
D 2 |
по |
|||||||||||
|
|
|
|
X1, , X n . |
|
|
|
|
|
|
|
|
||
выборке X |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
Z |
m* m |
|
||||
Образуем |
|
вспомогательную |
случайную величину |
|
|
|
|
, |
где |
|||||
|
|
|
|
|||||||||||
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
n |
|
|||
|
|
|
1 |
|
n |
|
|
|
|
|
|
|
|
|
m* X |
|
|
|
|
|
|
|
|
||||||
|
|
|
X i - точечная оценка математического ожидания m . Согласно |
|||||||||||
|
|
|
||||||||||||
|
|
|
|
n i 1 |
|
|
|
|
|
|
|
|
||
утверждению 1 теоремы Фишера, случайная величина |
Z имеет нормальное |
|||||||||||||
распределение |
N 0; 1 и ее функция распределения |
z P Z z не |
||||||||||||
зависит от неизвестного параметра.
Доверительный интервал, соответствующий надежности β, определяется из условия (3.20), которое в нашем случае имеет вид
|
|
P |
|
|
|
m* m |
|
|
. |
|
(3.31) |
|||||||||||||||
|
|
|
|
|
||||||||||||||||||||||
Неравенства |
|
m* m |
|
и |
|
|
m* m |
|
|
|
|
|
z являются равносильными, |
|||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
/ |
|
|
|
|
/ |
|
|
||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
n |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
то есть для |
любой |
|
|
|
выборки x |
|
они выполняются или не выполняются |
|||||||||||||||||||
одновременно, поэтому соотношение (3.31) можно записать в виде |
||||||||||||||||||||||||||
|
|
|
|
m |
* |
m |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
P |
|
|
|
|
|
|
|
|
|
|
|
|
z |
. |
(3.32) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
90
