ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 08.02.2025
Просмотров: 1280
Скачиваний: 2
^ |
||||||||
Если все выборочные значения x1;:::;xn различны, то функцию F (x;~) можно записать в виде |
||||||||
0; |
x x(1); |
¡ |
||||||
> n |
||||||||
Fn(x) = |
8 |
i |
; |
x(6i) < x 6 x(i+1); i = |
1;n |
1; |
||
> |
||||||||
: |
1; |
x > x(n); |
||||||
< |
||||||||
т.е. в каждой точке x(i) функция Fn(x) имеет скачок величиной 1=n. График функции Fn(x) изображен на рис. 14.1.
При больших объемах выборки n (n > 50)обычно |
|
производят группирование исходных данных следу- |
|
ющим образом. Промежуток J = [x(1);x(n)], содержа- |
|
щий все выборочные значения, разбивают на m полу- |
|
интервалов J1;:::;Jm, как правило, одинаковой дли- |
|
ны ¢ и таких, что каждый из них, кроме послед- |
|
него, содержит левую границу, а последний содер- |
|
жит обе границы, и подсчитывают число ni элемен- |
|
тов выборки, попавших в i-ый промежуток Ji, i = 1;m, |
|
n = n1 + ¢¢¢ + nm, а результаты представляют в виде |
Рис. 14.1. |
следующей таблицы 14.2, которую называют интер- |
вальным статистическим рядом. |
||||||||
Иногда в верхней строке таблицы 14.3 указывают не интервал, а |
||||||||
J1 |
J2 |
. . . |
Jm |
|||||
его середину xi, а в нижней строке вместо частоты ni записывают от- |
||||||||
. . . |
m |
|||||||
носительную частоту ni=n. |
n1 |
n2 |
nm |
ni |
||||
i=1 |
||||||||
промежутков m, на которые разбивают промежуток J = |
||||||||
Число |
e |
Таблица 14.3. |
P |
|||||
[x(1);x(n)], содержащий все выборочные значения, выбирают в зави- |
||||||||
симости от объема выборки n. Для ориентировочной оценки величины m можно пользоваться следующей приближенной формулой.
m ¼ log2 n + 1; |
|||||
Определение 14.8 График функции |
8 n¢i |
x 2 Ji; |
(14.3) |
||
pn(x) = |
; |
||||
n |
|||||
: |
0; |
2 |
|||
< |
x = J; |
||||
представляющий собой кусочно постоянную функцию называют гистограммой (см. рис. 14.2).
Часто гистограммой называют диаграмму, составленную из прямоугольников с основанием ¢ и высотами ni=(n¢), i = 1;m. Нетрудно увидеть, что суммарная площадь всех прямоугольников, образующих такую диа-
грамму, равна 1, так как Pm nn¢i ¢ = n1 Pm
i=1 i=1
каждого прямоугольника ni=n есть частота попадания элементов выборки в соответствующий интервал Ji статистического ряда.
~
Рассмотрим случайную величину ni(Xn)=n, которая для
~
каждой реализации ~xn случайной выборки Xn ~ равна частоте ni=n. В соответствии с законом больших чисел в форме Бернулли ni(Xn)=n при n ! 1 будет сходиться по вероятности к вероятности попадания случайной величины X в промежуток Ji, i = 1;m, т.е.
~
ni(Xn)
n
© |
ª |
Ji |
|
P |
= Z p(x)dx; |
||
n¡!!1 P X 2 Ji |
|||
где p(x) плотность распределения генеральной совокупности X. Если длина ¢ промежутков достаточно мала и объем выборки n велик, то с вероятностью, близкой к 1, можно утверждать, что
nni ¼ p(xei)¢;
или |
ni |
¼ p(xei); |
n¢ |
67
где xei середина промежутка Ji, i = 1;m. Таким образом, при большом объеме выборки n и достаточно малом ¢ с вероятностью, близкой к 1, можно считать, что pn(x) ¼ p(x). Иными словами, функция pn(x) является статистическим аналогом плотности распределения p(x), наблюдаемой в эксперименте случайной величины X, а гистограмма выглядит приблизительно как график плотности случайной величины X.
~
Пусть Xn случайная выборка из генеральной совокупности X с функцией распределения F (x) (и плотностью распределения p(x) в случае непрерывной статистической модели).
Случайную величину
1 |
n |
||||
~ |
Xi |
||||
k |
|||||
¹k(Xn) = |
n |
Xi |
(14.4) |
||
=1 |
|||||
называют выборочным начальным моментомb |
k-го порядка. В частности, выборочный на- |
||||
~ |
||||||||||||||||||||||||||||||||
чальный момент первого порядка X = ¹ (X ) называют выборочным средним. |
||||||||||||||||||||||||||||||||
Случайную величину |
b1 |
n |
1 n |
¡ |
¢ |
k |
||||||||||||||||||||||||||
~ |
Xi |
|||||||||||||||||||||||||||||||
ºk(Xn) = |
n |
=1 |
Xi ¡ X |
(14.5) |
||||||||||||||||||||||||||||
называют выборочным центральнымb |
моментом k-го порядка. В частности, выборочный |
|||||||||||||||||||||||||||||||
центральный момент 2-го порядка ¾ |
2 |
~ |
~ |
|||||||||||||||||||||||||||||
(Xn) = º2(Xn) называют выборочной дисперсией. |
||||||||||||||||||||||||||||||||
~ |
) = |
2 |
~ |
|||||||||||||||||||||||||||||
X |
¾b(Xn) называют выборочным средним квадратич- |
|||||||||||||||||||||||||||||||
Выборочную характеристику ¾(b n |
||||||||||||||||||||||||||||||||
ным отклонением. |
b |
q |
||||||||||||||||||||||||||||||
Случайную величину |
b |
1 |
n |
|||||||||||||||||||||||||||||
~ |
~ |
|||||||||||||||||||||||||||||||
b |
Xi |
¡ |
¢¡ |
¢ |
(14.6) |
|||||||||||||||||||||||||||
K(Xn |
;Yn) = n =1 |
Xi ¡ X |
Yi ¡ Y |
|||||||||||||||||||||||||||||
называют выборочным корреляционным моментом. |
||||||||||||||||||||||||||||||||
Выборочную характеристику |
||||||||||||||||||||||||||||||||
~ |
~ |
|||||||||||||||||||||||||||||||
K(X ;Y ) |
||||||||||||||||||||||||||||||||
½(X~n;Y~n) = |
¾x(bX~n)n¾y(nY~n) |
; |
(14.7) |
|||||||||||||||||||||||||||||
где |
nb |
¢ |
b |
b |
n |
¡ |
¢ |
|||||||||||||||||||||||||
1 |
X¡ |
1 |
Xi |
|||||||||||||||||||||||||||||
2 |
~ |
2 |
2 |
~ |
2 |
|||||||||||||||||||||||||||
¾x |
(Xn) = |
n |
i=1 |
Xi ¡ X |
; |
¾y |
(Yn) = |
n |
=1 |
Yi ¡ Y |
; |
|||||||||||||||||||||
называют выборочнымb |
коэффициентом корреляцииb . |
|||||||||||||||||||||||||||||||
Значения выборочных моментов (неслучайные числа) будем для краткости называть теми же терминами, что и сами моменты (случайные величины).
Основное свойство выборочных моментов состоит в том, что при увеличении объема выборки n
они |
сходятся по вероятности к соответствующим теоретическим моментам. В частности, при n |
! 1 |
|||||
P |
P |
||||||
2 |
~ |
||||||
имеем X n¡!!1 MX, а ¾b |
(Xn) n¡!!1 DX. |
||||||
68
Лекция 15
Точечные оценки
Одной из задач математической статистики является оценка неизвестных параметров выбранной
параметрической модели.
Предположим, что закон распределения генеральной совокупности принадлежит множеству
~ ~ ~
fF (x;µ) : µ 2 £g, где вид функции распределения задан, а вектор параметров µ = (µ1;:::;µr) неизве-
~
стен. Требуется найти оценку для µ или некоторой функции от него (например, математического ожидания, дисперсии) по случайной выборке (X1;:::;Xn) из генеральной совокупности X.
Например, предположим, что масса X детали имеет нормальный закон распределения, но его параметры µ1 = MX и µ2 = DX неизвестны. Нужно найти приближенное значение параметров по результатам наблюдений x1, . . . , xn, полученным в эксперименте (по реализации случайной выборки).
~ |
||||||
Пусть Xn = (X1;:::;Xn) случайная выборка из генеральной совокупности X, функция распре- |
||||||
деления F |
x;µ которой известна, а µ неизвестный параметр, т.е. рассматривается параметриче- |
|||||
~ ~ |
£ |
(для простоты изложения будем считать пока, что µ скаляр). |
||||
ская |
модель F (x;µ) : µ |
|||||
¡ f ¢ |
2 g |
~ |
||||
Требуется построить статистику µ(X ), которую можно было бы принять в качестве точечной |
||||||
оценки параметра µ. |
b |
n |
||||
Определение 15.1 Точечной оценкой параметра µ 2 £ назовем любую функцию от наблюдений
^ ~
(т.е. любую статистику) µ(Xn).
b ~
Определение 15.2 Статистику µ(Xn) называют состоятельной оценкой параметра µ 2 £, если с ростом объема выборки n она сходится по вероятности к оцениваемому параметру µ, т.е.
~ |
) |
P |
µ: |
||
µ(X |
n¡!!1 |
||||
b |
n |
||||
b ~
Определение 15.3 Статистику µ(Xn) называют несмещенной оценкой параметра µ, если ее
b ~
математическое ожидание совпадает с µ, т.е. Mµ(Xn) = µ для любого фиксированного n.
Замечание 15.1 Можно показать, что статистика
n |
|||||||||
¡ |
Xi |
¡ |
¢ |
2 |
|||||
S2 |
= |
1 |
Xi ¡ X |
||||||
n |
1 |
=1 |
|||||||
является несмещенной и состоятельной оценкой дисперсии DX генеральной совокупности. Ее называют исправленной выборочной дисперсией.
Замечание 15.2 Можно доказать, что выборочные начальные и центральные моменты являются состоятельными оценками соответствующих теоретических моментов, если только они существуют. Однако эти оценки, кроме X, являются смещенными.
Пример 15.1 Пусть X1, . . . , Xn случайная выборка из генеральной совокупности X, имеющей нормальное распределение с неизвестным средним значением µ и известной дисперсией ¾2.
b b
Оценка µ = µ(X1;:::;Xn) = X1 является несмещенной для µ, ибо MX1 = MX = µ, но не является состоятельной, так как, во-первых, X1 не зависит от объема выборки и, следовательно, ее распределение не меняется с ростом n, а во-вторых,
P jX1 |
¡ µj < " = ¾p2¼ |
Z |
e¡2¾2 |
dt 6= 1: |
|
© |
ª |
2 |
" |
t2 |
|
0 |
|||||
69
Метод моментов
Метод моментов был предложен английским статистиком К. Пирсоном и является одним из
первых общих методов оценивания. Он состоит в следующем. |
|||||||||
~ |
;:::;Xn) из генеральной совокупности X, распреде- |
||||||||
Пусть имеется случайная выборка Xn = (X1 |
|||||||||
~ |
~ |
;:::;µr). Требуется найти |
|||||||
ление которой p x;µ |
известно с точностью до вектора параметров µ = (µ1 |
||||||||
~ |
~ |
. |
|||||||
оценку |
параметра µ по случайной выборке X |
n |
|||||||
¡ ¢ |
|||||||||
Обозначим |
~ |
k |
~ |
k |
|||||
¹k(µ) = M(X) ; ºk(µ) = M(X ¡ MX) |
|||||||||
начальный и центральный моменты порядка k, k = 1;2;:::. В методе моментов в качестве точечной
~ |
~ |
~ |
~ |
~ |
r |
|||||||||||||
µ(X ) = (µ |
(X |
n |
);:::;µ (X )) |
вектора~параметров |
µ |
берут решение системы |
уравнений |
|||||||||||
оценки b |
n |
1 |
r |
n ~ |
||||||||||||||
b |
b |
¹i® (Xn) =¹i® (µ); ® = |
||||||||||||||||
1;k; |
||||||||||||||||||
( |
º (X~ |
) =º |
(µ~); ¯ = |
1;l; |
; |
k + l = r; |
(15.1) |
|||||||||||
b |
n |
j¯ |
||||||||||||||||
bj¯ |
||||||||||||||||||
относительно неизвестных µ1;:::;µr. Индексы i® и j¯ выбирают так,чтобы система уравнений решалась как можно проще. Можно показать,что при условии непрерывной зависимости решения этой системы от ¹bi® и ºbj¯ , оценка, полученная методом моментов, является состоятельной и имеет асимптотически нормальное распределение, т.е. ее распределение при n ! 1 стремится к нормальному. При этом уравнения (15.1) во многих случаях просты и их решение не вызывает больших вычислительных сложностей.
Понятно, что метод моментов не примен´им, когда моменты генеральной совокупности нужного порядка не существуют (например, для распределения Коши, у которого не существует даже начальный момент первого порядка математическое ожидание).
Пример 15.2 Пусть дана случайная выборка (X1;:::;Xn) объема n из генеральной совокупности X, имеющей равномерный закон распределения
p(x) = |
1 |
; x 2 (a; b); |
|
b |
a |
||
( |
¡0; x = (a; b); |
||
2 |
|||
с неизвестными параметрами a и b. Найдем методом моментов точечные оценки этих параметров. Известно, что для равномерно распределенной случайной величины X
MX = |
a + b |
; |
DX = |
(b ¡ a)2 |
: |
||||||||||||||||||
2 |
12 |
||||||||||||||||||||||
выборочная дисперсия ¾ |
2 |
~ |
|||||||||||||||||||||
X |
|||||||||||||||||||||||
Выборочное среднее |
и |
1 n |
2 |
b~ |
n |
1 |
n |
2 |
|||||||||||||||
X |
Xi |
¡ X) : |
|||||||||||||||||||||
X = n i=1 Xi; |
¾ (Xn) = n =1(Xi |
||||||||||||||||||||||
b |
|||||||||||||||||||||||
Составляем систему двух уравнений |
|||||||||||||||||||||||
8 |
2 |
2 |
|||||||||||||||||||||
> |
a + b |
= x; |
|||||||||||||||||||||
(b ¡ a) |
= ¾2(~x ): |
||||||||||||||||||||||
< |
12 |
n |
|||||||||||||||||||||
> |
b |
||||||||||||||||||||||
: |
|||||||||||||||||||||||
Решая систему, получаем |
p |
||||||
Окончательно имеем |
b= x + |
3¾b(~xn); |
|||||
~ |
p |
~ |
|||||
b(Xn) = X + 3¾b(Xn);
p
ba = x ¡ 3¾b(~xn):
p
~ ~
ba(Xn) = X ¡ 3¾b(Xn):
Метод максимального правдоподобия
Одним из наиболее универсальных методов оценивания параметров является метод максимального правдоподобия (предложенный Р. Фишером), суть которого состоит в следующем.
~ |
из генеральной совокупности X, |
|||||
Рассмотрим функцию правдоподобия случайной выборки Xn |
||||||
~ |
~ |
£: |
||||
p x;µ |
точностью до параметра µ |
2 |
||||
плотность распределения |
¡ ¢ которой известна с |
n |
¢ |
|||
~ |
iY ¡ |
|||||
L(X1;:::;Xn;µ) = |
p Xi;µ : |
|||||
=1 |
||||||
70